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ABSTRACT 

The starting point of this investigation is the partially transparent 
echo from midlatitude sporadic-# layers as recorded by ionosondes; this is 
the range of frequencies between the blanketing frequency and the maximum 
frequency. One interpretation is that it is a partial reflection at the 
steep gradients of electron density which characterize the vertical profile 
of the layers as recorded in rocket flights. The theory of reflection at 
gradients is developed and evaluated in the case of five layers. It is 
found that this mechanism is not sufficient to explain the phenomenon. An 
alternative explanation is that the midlatitude sporadic-# layers are not 
uniform in the horizontal plane but contain localized regions of high electron 
density. This is investigated using data obtained by incoherent-scatter 
radar and is found to provide a satisfactory explanation. 

The main features of midlatitude sporadic-# layers are consistent with 
the convergence of metallic ions as describ-d by the wind-shear theory 
applied to gravity waves and tides. The interference of gravity waves with 
other gravity waves and tides can be recognized in the altitudes of 
occurrence and the structure of the layers. The small-scale horizontal 
irregularities are attributed in some cases to critical level effects and 
i.i others to fluid instabilities. The convergence of a meteor trail can, 
under some circumstances, account for localized enhancement of the electron 
density in the layer. 
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1 . INTRODUCTION 

1.1 General Introduction 

The ionosphere, a region of partially-ionized gases in the earth's upper 
atmosphere, has been the subject of much research in the past years. Among 
the most striking features encountered in nearly every technique by which 
the ionosphere is studied are sporadic-E layers. 

Sporadic-E layers are common features in the region of the ionosphere 
between the altitudes of 95 and 120 km. This is the lower part of the E region, 
an ionospheric region extending from 90 to 160 km in altitude. An electron- 
density profile of the E region containing a sporadic-E layer is shown in 
Figure 1.1. Detailed vertical electron -density profiles measured by sounding 
rockets show midlatitude sporadic E to be thin layers of enhanced ionization 
[Smith and Mechtly, 1972]. The layers are usually only one or two kilometers 
thick, but may extend horizontally over hundreds of kilometers. 

The ionosphere at altitudes where sporadic E occurs contains high 
concentrations of metallic ions with complementary electrons to preserve 
charge neutrality [ Narcisi , 1973; Zbinden et al. t 1975]. The metallic ions 
are presumed to be of meteoric origin. Most are ionized either by collisional 
processes as the meteor enters the atmosphere, or by charge exchange with 
ambient molecular ions. The lifetime with respect to recombination of a 
metallic ion in the E region is several months. 

The most widely held theory of the formation of midlatitude sporadic-E 
layers is that of vertical convergence of the long-lived metallic ions by the 
combined action of the horizontal neutral wind and the geomagnetic field. The 
ions collect in thin layers determined by vertical shears in the horizontal wind. 

Layers of metallic ions are apparently a constant feature of the E region. 
However, the concentration in the layers is not always enough in daytime to 
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produce layering in the electron density. LaLonde [1966] observed layers of 
ionization being exposed at sunset as the E region ionization density decayed. 
This was also observed during a solar eclipse by Narcisi [1966], Layering of 
• the electron concentration is nearly always in evidence in the nighttime 

ionosphere [Smith and Mechtly , 1972]. 

Sporadic E was first studied by reflecting radio waves from the ionosphere 
and examining the characteristics of the reflected wave. It is seen as a strong 
reflection from the E region at a virtual height that is constant with frequency. 
The echoes from sporadic E are detected at frequencies much greater than the 
plasma frequency of the normal E region. As these echoes are not always present, 
appearing in fact to occur randomly, they were given the name, "sporadic £". 

Most observations have been made using an ionosonde, a type of ionosphere 
sounder which measures the time it takes a transmitted pulse of electromagnetic 
radiation to be reflected by the ionosphere and return. The time delay gives a 
measure of the virtual height of the reflecting ionization, which can be related 
to the true height of reflection by considering the retardation of the group 
velocity of the radiation passing through the ionized medium. The height of 
reflection is the height at which the local plasma frequency of the ionosphere 

equals the frequency of the incident radiation. By varying the frequency of 

i 

the pulse, the height and electron density of the underside of the ionosphere 
can be measured. The ionosonde typically makes a record of the virtual height 

l 

i of reflection versus frequency. Examples of this type of record, called 

( 

I 

iono grams, are shown in Figure 1.2. 

I • 

As early studies of sporadic E were done, for the most part, with ionosphere 
sounders, the terminology describing sporadic E describes its appearance on iono- 
grams. An idealized ionogram illustrating sporadic-E nomenclature is shown in 
Figure 1.3. The blanketing frequency (fbEs ) is defined as that frequency below 

l i 
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fmin 


fbEs foE* fxEt 


Figure 1.3 An idealized ionogram illustrating sporadic-^ nomenclature 
[URSI, 1972]. 
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which no traces are seen of reflections from above the sporadic-# layer. The 
top frequency [ftEs ) , sometimes called the critical or maximum frequency (fEe ) , 
is the frequency above which no echo is detected from the sporadic-# layer. If 
the maximum frequencies of reflection of ordinary and extraordinary waves can be 
resolved, they are designated by foEs and fxEs respectively. It must be empha- 
sized that these conventions were developed to document features on ionograms, and 
do not explain the physical features of the sporadic-# layers. 

Other nomenclature has been adopted to describe the appearance and time 
development of the echoes from sporadic # [Smith, 1957]. The various types of 
sporadic # are grouped into three major categories: auroral, equatorial, and 

midlatitude. Although ionograms of these three types of sporadic # share many 
characteristics, the physical phenomena involved are quite different. Auroral 
sporadic # is caused by ionization of the atmosphere by auroral particle precipi- 
tation [ Hunsucker , 1972] , while equatorial sporadic # is the result of the 
scattering of radio waves by the electron-density irregularities of the unstable 
plasma of the electrojet [Balsley et al 1976], In this study, only sporadic # 
occurring at midlatitudes will be considered. 

The detailed vertical structure of sporadic-# layers has been observed 
using sounding rockets. This is usually done by measuring the current to a 
Langmuir probe operated at a fixed voltage [Smith, 1969]. The layers are 
usually Gaussian-shaped, but more complex structures are sometimes observed 
[Smith and Meehtly, 1972] . 

In spite of the many observations of sporadic # by these and other tech- 
niques, many features of sporadic-# layers remain unexplained. It is the purpose 
of this work to examine the structure and motion of sporadic-# layers in detail, 
both in the vertical and horizontal planes, in the context of the present theory 
of the formation of these layers. 



1.2 Partially Transparent Echoes from Midlatitude Sporadic- E layers 

On an ionogram there is a frequency range between fbEs and ftEs where the 
sporadic-E layer is partially transparent to radio waves. Often this range is 
too large to be explained on the basis of the different reflection frequencies 
of the ordinary and extraordinary modes. Models proposed to account for this 
effect have been discussed by Reddy and Rao [1968] and by Whitehead [1972a]. 

In one model the echo is attributed to partial reflection at the large vertical 
gradients of electron density that characterize the upper and lower boundaries 
of a sporadic-E layer. Calculations by Reddy [1968] and by Chessell [1971a, b] 
indicated that partial reflections from thin layers can occur over a large range 
of frequency and provide an explanation for the effect. However, re-examination 
of this model in the present work using observed electron-density profiles of 
sporadic-E layers shows that this explanation is, in fact, inadequate. 

In another model, the partial transparency is attributed to high electron 
density in regions of limited horizontal extent in the sporadic-E layer. In 
this model the maximum frequency of the sporadic-E reflection, ftEs, corresponds 
to total reflection from the regions of enhanced ionization illuminated by the 
ionosonde. The blanketing frequency fbEs then corresponds to the peak plasma 
frequency of the layer in which these regions are embedded. 

The presence of irregularities of this type within sporadic-E layers has 
not been substantiated by rocket observations. Reddy and Rao [1968] showed that 
the peak plasma frequency of sporadic-E layers measured from rockets usually 
approximated fbEs, and was never as large as ftEs. Further, it was shown by 
Smith and Mechtly [1972] that the sporadic-E electron -density profile measured 
during the descent of the rocket usually matches the profile from the ascent in 
both shape and altitude, even though a distance of 100 km or more often separates 
the two points on the trajectories. However, this in itself is not sufficient to 


exclude the idea of horizontal variability in sporadic-E layers. It does 
support the idea that, if there are irregularities in sporadic-E layers with 
the required electron densities, they must be of small horizontal extent and 
sparsely scattered throughout the layer; the probability of a rocket passing 
through an irregularity must be very small. 

In a variation of this second model, the irregularities are small-scale 
fluctuations appearing in the vertical electron-density profile. Data of 
nighttime sporadic E taken by rocket-borne probes have indicated the presence 
of irregularities such as these which could scatter radio waves with 
frequencies much higher than would be reflected. Bouhill [1966] has shown 
that the partially transparent sporadic-E echo observed by the ionosonde on 
those occasions can be accounted for theoretically by scattering from the 
irregularities on the electron -density profile. However, this fine structure 
is not generally observed by rocket probes when ftEs » fbEs. 

1.3 Objectives of this Study 

The main objectives of this investigation are to explain the partial 
transparency of midlatitude sporadic-E layers to radio waves, and to evaluate 
the results with respect to the dominant physical processes involved in the 
formation and maintenance of sporadic-E layers. 

Chapter 2 starts with the development of the general theory of the 
formation of sporadic-E layers by wind shears. This includes a critical 
review of the important physical and chemical processes in the lower E region 
of the ionosphere as they apply to sporadic E. 

Because of their apparent importance to sporadic E, a section of this 
chapter is devoted to consideration of metallic atoms and ions. The basic 
chemical reactions involved are discussed and an estimate made of the product 
rate of metallic ions. A new simplified description is developed of the 



diurnal variation of the metal content at sporadic-# altitudes. 

Chapter 3 deals with the first of the models for the partial transparency 
of sporadic-# echoes. This is the model which attributes the partial reflec- 
tions to the large electron-density gradients that characterize the vertical 
profile of the layers as recorded in rocket flights. The theory of reflection 
at gradients is developed and evaluated in the case of five sporadic-# layers. 
It is found that this mechanism is not sufficient to explain the phenomenon. 

The alternative model, that sporadic-# layers are not uniform in the 
horizontal plane but contain b ilized regions of enhanced electron density, 
is then investigated. As described in Chapter 4, this is done using data 
taken by the incoherent- scatter radar at the Arecibo Observatory in Puerto 
Rico. The structure of many of the layers observed is complex, and includes 
irregularities of the magnitude required to account for the partial trans- 
parency of the sporadic-# layers. 

Various interpretations of the observed irregularities in sporadic # are 
discussed in Chapter 5. The complex structure of sporadic-# layers is 
interpreted in terms of the effects of motions of the neutral atmosphere. 
Evidence is shown that internal gravity waves play a major role in the develop- 
ment of this structure and are potentially able, in a region of high wind 
shears, to converge ionization horizontally and form highly enhanced irregu- 
larities in the electron density. The action of wind shears on meteor trails 
in the vicinity of sporadic-# layers is also shown to be one possible source 
of isolated regions of enhanced electron density. 

It is shown that sporadic-# layers do contain highly enhanced electron- 
density irregularities which are apparently responsible for the partial 
transparency of the layers to radio waves. The complexity observed in the 
electron- density distribution is interpreted in most cases as a complexity in 
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the motion of the neutral atmosphere. It is also concluded that in spite of 
the complexity of the sporadic-F layers, the features observed can be 
explained by applying wind*shear theory to metallic ions. 




11 


2. THEORY OF FORMATION OF SPORADIC-E LAYERS 

2.1 Introduction 

Nearly all sporadic-E layers that are observed occur in the altitude region 
between 95 and 120 km. This is an interesting region of transition wher<" the 
relative importance of many processes is changing. In the context of electri- 
cal conductivity it is sometimes called the dynamo region. It will be shown 
later that the mechanism producing vertical ion drifts, which lead to the forma- 
tion if sporadic-E layers, is based on the same concepts as are used in 
deriving the conductivity tensor. 

The atmosphere of the earth is divided into natural regions by various 
parameters such as ionization density and temperature. The relative locations 
of several of these regions are shown in Figure 2.1. The lower thermosphere, 
where sporadic E is observed, is a region of transition between dynamic regimes. 
The turbopause, shown in Figure 2.1 at 105 km, marks a division between turbu- 
lent mixing of the gases below, and non-turbulent gravitational (diffusive) 
separation above. 

The ionospheric D, E, and E regions are determined by layers of ionization 
which are not always well defined. The layering is the result of the photo- 
ionization of different atmospheric species. As the optical depth of the 
atmosphere increases and the composition changes, the primary species being 
photoionized and the ion production rate also change. The ionization of the 
daytime E region is predominantly 0^ and N0 + , with 0 reaching comparable 
concentrations above 110 km [ Whitten and Poppoff , 1971]. The basic theory 
of the formation of a layer by photoionization was first developed by 
Chapman [1931]. 

The lower E region also contains metallic ions and atoms. The average 
concentration of ionized metals is small compared with O 2 and NO 
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concentrations. However, the recombination rate of metallic ions is very 
small, allowing the convergence of the metallic ions into thin layers with 
peak densities much greater than that of 0^ + and N0 + . At night, the decrease 
in ionization density from the recombination of the molecular ions reveals 
layering in the electron-density profile characteristic of metallic ions on 
essentially every occasion [Smith and Mechtly 3 1972]. 

The formation of sporadic-# layers by wind shears was first suggested by 
Dungey [1956; 1959]. The theory has been developed by several authors, 
including Whitehead [1961], Axford [1963], MacLeod [1966], and Kato [1966], 

The layering of metallic ions is the subject of discussion of this chapter. 

The theory of formation of sporadic-# layers by wind shears will be developed, 
together with a discussion of the source and characteristics of the metallic 
ions. We begin with a discussion of the motions of the neutral atmosphere, 
which is the driving force for the formation of these layers. 

2 . 2 Wind- Shear Theory 

2.2.1 E-region winds. The lower E region has been found to be 
characterized by variable winds and large wind shears [Bedinger et at. 3 1968; 
Rosenberg and Zimmerman 3 1972]. From a series of vapor-trail releases from 
rockets, Rosenberg [1968] found an r.m.s. wind speed at sporadic-# altitudes 
of 70 m s *, and an r.m.s. wind shear of 120 ra s * per scale height, or 
approximately 0.02 s 1 . Results of the Rosenberg study, plotted in Figures 
2.2 and 2.3, give an indication of some of the characteristics of E- region 
winds. The mean vertical wavelength of the dominant w;:" was found to be 3.4 
scale heights, or about 21 km at an altitude of 100 km, with no vertical 
dependence of the amplitude. No significant seasonal dependence of either 
the speed or the vertical shear was detected. 

It is difficult to determine the prevailing winds because of the 
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variability of this region of the atmosphere. Meteor radar experiments 
indicate that prevailing winds are eastward in midlatitudes, increasing 
with altitude in the summer and decreasing with altitude in the winter. There 
is evidence that at certain latitudes the prevailing winds become westward at 
higher altitudes in the winter [Hines, 1963] . 

The motions of the atmosphere at sporadic-E altitudes are found to 
include tidal oscillations and gravity waves [Hines, 1974; Francis, 1975]. The 
two types of oscillations are actually very similar, differing mainly in period 
and horizontal wavelength. The main energy source for tides is solar heating, 
with some small effect contributed by the lunar gravitational field. Sources 
of gravity waves are not well known, but are more localized and are thought to 
include geographic features, shear instabilities in the atmosphere, turbulence, 
auroras, and solar eclipses [Hines, 1974--Preamble to paper 22]. 

Expressions describing both tides and gravity waves are derived [Rishbeth 
and Garriott, 1969] by combining the equation of motion of the neutral 
atmosphere, the continuity equation, the adiabatic gas equation, and the 
perfect gas law. In the case of tides the equations are solved in spherical 
coordinates. It is found, for example, that tidal oscillations of solar origin 
can be written in polar coordinates (r, 0, <}>) as 


= I E 
m n 
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mn 
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( 2 . 1 ) 


where x is the velocity divergence. The angular functions in equation (2.1) 

8 

are Hough functions and can be written in terms of spherical harmonics (Y^) as 

%<«**"* - llm hm t2 ' 2 > 

The radial term, E(r) , can be approximated locally by 

E (r) ' exp[ik r) (2.3) 
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where the wave number k is given by 

nn 

The scale height of the neutral atmosphere is represented by H, and h is 
the "equivalent depth" given by the eigenvalues of the various modes. If 
is real, the wave is internal and propagates vertically. If k is 
imaginary, the wave is evanescent, k ^ is imaginary for all waves in which 
the subscript n is less than zero. For n > 0, k can become imaginary at 
certain heights, for example where dH/dh becomes sufficiently negative, 
indicating that the wave is reflected. 

The energy density of the wave can be shown to be independent of height. 
Since the energy density can be expressed as 

E = h d\u\ 2 (2.5) 

where d is the mass density and u is the velocity amplitude, the amplitude of 
the tidal oscillations is inversely proportional to the square root of the 
density. This increase in amplitude with height is known as tidal amplification. 
A similar effect can also be derived for gravity waves. However, this increase 
of amplitude with height is not evident in the wind reported by Rosenberg [1968]. 

The behavior of internal gravity waves is described by the same equations. 
However, since the periods of gravity waves are much less than the earth's 
rotational period and their wavelengths are much less than the radius of the 
earth, both tidal and Coriolis effects can be neglected. The dispersion relation 
can then be expressed as 

2222 222 22 
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where k and k are the horizontal' and vertical components of the wave number 
x z 


and 


w a = yg/2c 


u g = (Y-l) g/c 


(2.7) 

( 2 . 8 ) 


c = (y gH) (speed of sound) (2.9) 

The dispersion relation, equation (2.6), shows that propagating solutions can 

exist only for waves with frequencies either greater than w , in which case 

d 

they are known as acoustic gravity waves, or less than w , which are known as 

y 

internal gravity waves. At sporadic-^ altitudes, corresponds to a period 

of about 4.5 min, and u to a period of about 5 min. 

& 

Although gravity waves have greater downward phase velocities than tides, 
the vertical wavelength of both tides and gravity waves are similar. This is 
also true of velocity amplitudes. Hines [1966] found a peak amplitude of about 
40 m s’ 1 . Kochanski [1964] estimates the amplitudes of gravity waves to be 
55 m s’ 1 . Whitehead [1972bJ considers this value of the amplitude of gravity 
waves to be an underestimate, since meteor-radar observations show the amplitude 
of the irregular winds to be greater than that of either tidal modes. In any 
case the amplitudes are roughly comparable, and very complex wind patterns and 
changes can arise from the superposition of the two types of oscillations. 

2.2.2 The equation of motion of a charged particle. The motion of a 
charged particle in the E region is, to a large degree, influenced by the effects 
of its collisions with neutral atoms and molecules, by electric fields and by the 
geomagnetic field. The equation of motion can be written 


dv -► ♦ -*-i 

m * mv(w- 1>) ♦ qv*B + qE + mg - — Vp 


( 2 . 10 ) 


where 


q 3 m t v - charge, mass, and velocity of the charged particle 
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■f 

w = velocity of the neutral wind 
v = effective collision frequency 
E = electric field 
B = geomagnetic field 
g = gravitational acceleration 

n, p = number density and partial pressure of the particular 
charged species. 

The electric field, , is the sum of any external electric fields and the 
internal (polarization) electric field resulting from forces acting with 
different effects on the electrons and ions. Thus, 


E = t 


ext 



( 2 . 11 ) 


In this analysis, the magnetic and gravitational fields are assumed to be 
constant for a limited range of altitude. The geomagnetic field is written 

$ = B q r (2.12) 

A 

where T is a unit vector in the direction of the field. 

The effective collision frequency is a simplified expression, and is used 
in the sense that mv(w-y) represents the force per unit volume exerted on the 
ions by collisions with the neutral gas. The large ratio of neutral-to-ionized 
particle concentrations in the dynamo region justifies neglecting collisions 
between ionized particles. The force due to the partial-pressure gradient is 
also a macroscopic property of the gas, and can be written as n *V(nW), where 
k is Boltzmann's constant and T is the temperature of the ionized species. 

This term is much smaller than the other terms in equation (2.10) and is 
usually neglected [MacLeod et al. 3 1975], 

It is convenient, in solving equation (2.10) for the drift velocity of 
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the ionization, to define 


p = V/O) 


'H 


(2.13) 


where is the gyTofrequency, given by 
qB 

- 0 
= ~ 


(2.14) 


and to combine all terms not containing v in the following manner 


u = w + 


p qB 


| - [q£ + mg - (nW) ] 


(2.15) 


Equation (2.10) can now be written, assuming steady state conditions. 


0 - p(u-y) ♦ yxf 


(2.16) 


which can be solved for y by taking, respectively, the dot product and the 
cross product with r and then combining the results, giving 

y = — [p 2 u * pu*r + (u-r)r] (2.17) 

1+p 

Equation (2.17) represents the most concise formulation of wind- shear 
theory, including the effects of the neutral wind system and of external 
electric fields. From equation (2.15) it can be seen that the electric field, 
if known, can be easily incorporated into the calculation of ion velocities. 

This expression shows the relative importance of the various components 
of the velocity of the neutral wind and of the electric field at different 
altitudes. In the case of ions p is equal to one near an altitude of 120 km 
[Wright and Fedor , 1970], At this altitude, the three terms in equation (2.17) 
are of equal magnitude. At lower altitudes, collisions with the neutral gas 
become more important, and the first term dominates. At higher altitudes, 
the third term dominates and the ions are constrained to move parallel to the 
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magnetic field. In the case of electrons, p is much less than one throughout 
the region of interest, and the electron velocity is dominated by the third 
term in equation (2.17). 

The analysis so far has been perfectly general. If winds are neglected 
(i.e. w = 0 in equation 2.15) and only the electric field is considered, 
coefficients of equation (2.17) give expressions proportional to elements of 
the conductivity tensor. Thus, equation (2.17) can be written 

y = ^- [o p ^ + o/*r + (° 0 -°p)(^*r) r] (2.18) 


where [Whitten and Poppoff, 1971] 

(Pedersen conductivity) 
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(Hall conductivity) 


(2.19) 
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(longitudinal conductivity) 


and Up is the plasma frequency. By choosing a coordinate system with the 
magnetic field oriented along the a-axis equation (2.18) leads to the 
familiar form 
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where j is the current density. 

2.2.3 The Polarization electric field. The equations derived above have 
neglected any interacticr-s between electrons and ions. However, when 
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equation (2.17) is applied separately to ions and electrons it is seen that 
they move in different directions creating a polarization electric field. 

This is usually neglected in calculating vertical drifts for midlatitude 
conditions because of the greater mobility of electrons parallel to the 
magnetic field; when ions move vertically, the electrons are able to follow. 

Explicit expressions for the components of the polarization field are 
given in Appendix I. The field is shown to be negligible at midlatitudes. 
However, the effect of the polarization field becomes appreciable at low 
latitudes, preventing the formation of sporadic-^ layers by vertical shears 
in the horizontal wind at the geomagnetic equator. 

The correlation of sporadic E with latitude was suggested by Smith [1957] 
and shown by Heislev and Whitehead r 1 960; 1964] using data from the study 
by Smith [1957]. The absence of blanketing sporadic E at the magnetic 
equator was noted during the I.G.Y. by Bowles and Cohen [1962] from vertical 
incidence ionograms. 

2.2.4 The formation of ionized layers. Vertical ion drift velocities 
can be calculated using equation (2.17) if the neutral wind velocity and 
the external electric field are known. The neutral wind is usually determined 

by the motion of vapor trails from rockets, as described in Section 2.2.1. 

The electric field in the midlatitude ionosphere has been measured by balloon- 
and rocket-borne probes [Moser, 1972; Sahutz et al. , 1973], and inferred from 
incoherent scatter observations [Harper et al., 1976]. The field is horizontal 
and equal to or less than 1 or 2 mV m ^ 

Equation (2.15) shows that an external electric field has the same effect 
on the ionization as does a neutral wind with a velocity of £/p B ^ For an 

electric field of 1 mV m \ this term would be equal to 25/p ms *, or 8 m s 1 

at an altitude of 110 km. Although the external electric field is small, it can 
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change the altitude of a sporadic-E layer by a few kilometers [Macleod et al. , 
1975]. 

Once the vertical ion velocities are found, the ion-density profile can 
be determined using the continuity equation 

3«/ 3t = Q - L + V*(EVn - *y) (2.21) 

where D is the coefficient of diffusion, and Q and L are production and loss 
terms for the species being considered. In the case of metallic ions, the 
production and loss by chemical reactions are small compared with the effects 
of diffusion and transport [Ferguson, 1972; Brown, 1973]. Thus, for 
steady-state conditions, equation (2.21) can be written 


lvn = v/D ( 2 . 22 ) 

If, for example, horizontal stratification is assumed and the vertical ion 

drift changes linearly with altitude as 

v = a(z-zj (2.23) 

z o 

and assuming the diffusion coefficient to be constant through the region, 
the solution of equation (2.23) is 


n - n° exp 




(2.24) 


Equation (2.5) describes an electron-density profile with a simple Gaussian 
shape, which is similar to many of the observed sporadic-E profiles. 

Because of the complex wind structure and chemistry occurring in the 
lower E region, ionization-density profiles are usually found numerically 
using a digital computer. In this case, equation (2.21) can be solved for 
all species, including those for which production and loss terms arc important 
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[MacLeod et at., 1975]. Axford and Cunnold [1966] discussed the effect of 
a wind shear on two ionized species with differing recombination coefficients. 
They showed that, while the slowly-recombining species formed a layer at the 
wind shear, the density of the ion with the fast recombination rate was 
actually depleted at the wind shear. This is attributed to the increased 
recombination loss rate resulting from the enhanced electron density at 
the layer. The effect has been observed in im mass spectrometer measurements 
[Narcisi, 1968; Zbinden et al., 1975]. 

The requirement that the ionized species from which sporadic-E layers 
are formed have a small recombination rate lead.; to the consideration of 
metallic ions. The presence of metallic ions in the E region, and their 
role in the formation and maintenance of sporadic-E layers will be 
discussed in the following section. 

2.3 Metallic Iona 

2.3.1 Introduction. In the previous section the theory of formation 
of sporadic-E layers by wind shear was described. One of the major problems 
encountered in the early development of wind-shear theory of sporadic E 
was the lifetime of the ions in the layer [Layzer, 19b4a]. It was assumed 
that sporadic-E layers would contain the same ionic constituents as the rest 
of the lower E region, predominantly 0^* and N0 + . However, the rapid 
recombination rates of 0 o and NO would prevent the peak electron density 
in the layer from reaching the observed values. In response to this problem, 
a model was developed of sporadic-E layers formed of metallic ions [Axford 
and Cunnold, 1966]. The presence of metallic ions in a sporadic-E layer 
had been shown by Istomin [1963] using a rocket-borne mass spectrometer. 

In this section, various sources and sinks of metallic ions are discussed. 
It is suggested that direct ionization of metal atoms at the time a meteoroid 
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enters the earth : s atmosphere is an important source of metallic ions. A model 
is developed describing the diurnal variation in metallic ion concentration 
to conform with a similar behavior in sporadic-# peak electron density. The 
model considers charge-exchange reactions between metal atoms and the major 
ionospheric constituents, as well as the direct contribution by ablating 
meteors. The ionization is removed by transport downward out of the lower 
E region. 

2.3.2 General characteristics of metallic ions. Mass-spectrometer 

measurements have shown the presence of layers of metallic ions in the 

region between about 85 and 130 km [ Istomin , 1963; Narcisi and Bailey , 1965; 

Narcisi, .1968, 1972; Young et al., 1967; Anderson and Barth , 1971; 

Johanessen and Krankovsky, 1972; Zbinden et al ., 1975; Goldberg , 1975]. 

+ + 

The most common ions in these layers are Fe and Mg . Also present arc 

+ + ♦ + + + 

Si , Na , Ca , K , and traces of others. Si has the same mass as N 2 

but is identified above 100 km by the low recombination rate required for 

the formation of layers and at lower altitudes by its chemical properties 

[Goldberg, 1975]. 

It is generally accepted that the metal ions are of meteoric origin, 
although the relative abundances of the metals in meteors are not. reflected 
in the metallic ions. Baggaley and Curmack [1974] calculate the relative 
abundances of the four main ions expected to form initially by collisions 
with the atmosphere as Fe + , 48%; Si + , 34%; Mg + , 10%; and N'a + , 7%. Other 
theories now discounted of the source of the metal ions include the transport 
to high altitudes of particles initially found in sea spray r Allen, 1970], 
and the evaporation from aerosols of both terrestrial and extraterrestrial 
origin in the ionosphere [Donahue and Meier . 1967] . 

Considering ablation from meteors as the source, metal atoms are ionized 
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by throe processes: collisions with the ambient atmosphere as the meteor 

impacts; charge exchange with and NO ; and photoionization. Detailed 
chemical reaction schemes have been outlined in several papers [Swider, 

1969; Ferguson, 1972; Brown , 1973; Baggaley and Cummack, 1974; Poole and 
Nicholson, 1975]. Using iron as an example, the important reactions 
involving metals at sporadic-^ altitudes, with rate constants as quoted in 
Brown [1973], are 


Fe + h\j Fe + e 

J = 5 * 

IQ' 7 s'* 

Fe ♦ -*■ Fe + + 0., 

k = l.l 

x 10" 9 cm -5 s" 1 

Fe + N'0 + -*■ Fe + + NO 

*2 = 9.2 

x 10" 10 cm -3 s" 1 (2.25) 

+ — 

Fe + e -*■ Fe + hv 

*3 = 1 X 

in -12 3 - 1 
10 cm s 


Rate constants for magnesium are approximately the same as' those for iron 
[Brown, 1973] . 

2.3.3 Production rates of metallic atoms and ions. The production rate 

of metal ions by the initial collisions with atmospheric molecules has been 

computed by Lebedinets and Shushkova [1970]. They find a peak production 
-1 -3 -1 

rate of 10 cm s occurring at about 105 km. 

The ionization rate due to meteor trails can be estimated from meteor 

-2 -1 

radar observations of line electron density. The number N(m s ) of 
meteor trails exceeding a line electron density A (m *) is represented by 

A/(>A) = U~ l (2.26) 

2 -3 -1 

with K * 1.3 x 10 m s [Sugar, 1964]. This is established over the range 
10 10 < A < 10*** m 1. (jiff eren ti a i flux is given by 

dN/d A = -K A -2 (2.27) 
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Assuming all trails are vertical, the ionization production rate 
* s *) is given by 





A dN = K 




(2.28) 


2 -3 -1 

Thus, for each decade in line electron density, ^ = 3 * 10 ji s . For 

10 18 3 -3 -1 

the range of A from 10 to 10 then, ^ s 2.5 x 10 m s 

In addition there is an obliquity factor resulting from the distribution 

of zenith angles of the trails. Assuming this to be of the order of a 

factor of 10, the ionization rate due to metal ions formed directly in 

meteor trails is ^ = 3 x 10 4 m 3 s * or 0.03 cm 3 s *. 

A check on this estimate of is found in data from the daytime 

equatorial ionosphere. Woodman [1970] quotes the upward vertical ion 

velocities in the upper E region as approximately 20 m s *. Aiken and 

Goldberg [1973] show daytime concentrations of individual metal ion 

_3 

species at heights above the meteor region of about 10 cm . If all metal 
ions produced in the meteor region are transported upward, the total 
production rate of metal ions in that region can be found. From the 
continuity equation (2.21), 


nv _ _ . 
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total 

which is consistent with the production rates given above. 

Another estimate may be obtained from the total amount of metal 

deposited in the ionosphere by meteors. The amount is not well known. A 

representative estimate of the total mass of the matter deposited by 

ablation from meteoroids is 9 * 10^ g yr ^ [ Gadsden , 1968] which corresponds 

6 - 2 

to a column influx of metallic atoms and ions of approximately 2 * 10 cm s 
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[Broun, 1973]. If the effective depth of deposition of meteoric material 
is 30 km, this estimate gives a combined production rate of atoms and 
ions as 0.7 cm ^ s *. Comparing this to the production rate of metallic 
ions discussed above gives the fraction of metallic atoms which are ionized 
directly in meteor trails as approximately 5%. 

Little is known of the neutral metal concentrations. Sodium has been 
measured more successfully than other metals, but only at altitudes near 90 km. 
There are many opinions of the relative production rates of atoms and ions 
by meteor ablation. Based on measured cross sections for ionizing collisions 
between neutral species. Brown [1973] concluded that thermal ization of the 
atoms by kinetic energy transfer is much more important than ionization by 
the collisions. However, Lebedinets and Shushkova [1970] found that, based 
on their calculations, ionization by meteors is an important source of 
metallic ions in the upper atmosphere. The ionization of metal atoms in 
meteor trails by enhanced N0 + concentrations and increased radiation levels 
may also contribute significantly. Anderson and Barth [1971] compared day 
airglow measurements of Mg and Mg + , and concluded that the ratio, Mg + /Mg, 
has a lower limit of 22.2. Gadsden [1972], on the other hand, held that 
the measurement of Anderson and Barth might be low by as much as a factor 
of 20. If meteor ablation and transfer processes are neglected, it can 
be shown using the reactions of equation (2.25) that over 99% of the metal 
atoms would be ionized under either daytime or nighttime conditions [Brown, 1973]. 

2.3.4 Lifetimes of metallic ions. Recombination rate constants of 

-12 3 -1 

the metallic ions are on the order of 10 cm s The time rate of 
change of the concentration of metal ions by recombination is given by 


^ to 


( 2 . 30 ) 
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where [M + ] is the metallic ion concentration and k ^ is the recombination 
coefficient. For the situation in which [M + ] << [ e ~ ] this can be solved 
for [M + ] giving 


[M + ] = [M*] q exp(fe 3 [e _ ] t) 


(2.31) 


The time constant for recombination is then expressed as 


1 

t r = k^= T 


(2.32) 


Assume, for example, that [e ] is constant at the daytime value of 
10 5 cm 3 and [M + ] is 10 3 cm 3 , then it is seen that x^ = 10 7 s, or 
approximately 16.5 weeks. 

The residence time of metallic ions due to transport processes, x^,, 

can be estimated by the ion production rates discussed previously. If an 

3 -3 

average metallic ion concentration of 10 cm is assumed, then if the ion 

-3 -1 

production rate is 0.025 cm s , x- is 11.1 hr; and for a production rate 
of 0.01 cm 3 s * , 27.8 hr. 

The only effective loss process at midlatitudes is thus the downward 
transport of the ions by the effects of the neutral wind. Long-lived ions 
trapped at a null point in a gravity wave follow the downward-moving phase 
of the wave until the ions are "dumped" in the region below 90 or 95 km 
where collisions with the neutral atmosphere inhibit their motion and other 
chemical reactions become important [Chimonas and Axford, 1968]. The 
predominance of either the diurnal or semidiurnal tide in the transport 
of metallic ions from the E region would be allowed by the range of 
residence times of metallic ions discussed above. 

2.3.5 Diurnal variation. The diurnal variations in the concentrations 
of metallic atoms and ions may be obtained as the solutions of 
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and 


= q z + K[M] - T[M*] 


(2.33) 


(2.34) 


where 


K = 3 ♦ fe i[° 2 + ] + ^2 [N0+] (2.35) 

In these equations, [M] and [M + ] are the concentrations of metallic atoms 
and ions, and q A and are the rates of production of atoms and ions by 
meteors. The term K[M ] represents the combined effects of photoionization 
and charge-exchange with the ambient ionization. For a metallic ion profile 
that decreases exponentially with height, the loss process (downward transport) 
is proportional to [/7 + ] , T being the constant of proportionality. 

The diurnal variation is mainly the result of the variation of the 
charge-exchange process contained in the coefficient K. This varies 
continuously during a 24-hour period in a way determined by the concentrations 
of 0 2 + and N0 + . A simple model is used in the following analysis in 
which K is approximated by a constant nighttime value between sunset (18 hr) 
and sunrise (06 hr) and is set equal to a constant value, representative 
of average daytime conditions, between sunrise and sunset. The time constants 
are such that this model gives a representation of the day-to-night change 
in concentrations of metallic ions and atoms but does not reproduce 
accurately the transition periods near sunrise and sunset. 

For the situation discussed above the concentrations of metallic atoms 
and ions can be found as the solution of equations (2.34) and (2.35) 
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[M] = {[M] o -q A /K)e~ Kt ♦ q^K 


(2.36) 


[W + ] = 


* [A V«4\ -Kt 


([^ + ] 0 - 








(2.37) 


where [Afl and [M ] are the concentrations of atoms and ions at the 

beginning of the time period being considered. The concentrations at the 

end of the half-day period are found by putting t = 12 hr in these equations. 

The four equations for the concentrations at sunrise, [M] ^ and [M ] ^ , and 

at sunset, [M ^ ma Y then be solved simultaneously in terms of 

the quantities q^, q^, K(day), /((night), and T. 

Possible values of the sunrise and sunset concentrations are illustrated 

by the following examples. From equations (2.25) values of K can be 

found by estimating concentrations of 0_, and NO*. These concentrations 

are quite uncertain at the present time. For this example [0 ? + ] and [NO*] 

4-3 4-3 

will be taken to be 4 * 10 cm and 2 x 10 cm in the daytime; 

2-3 2-3 

3 x 10 cm and 6 k 10 cm at night. The corresponding values of K 

are 6.3 * 10 ^ s 1 in the daytime and 8.8 * 10 ^ s ^ at night. The 
concentrations of metallic atoms at sunrise and sunset are then 


[tf], = 6.05 x 10 q 
1 A 


(2.38) 


[M] = 1.88 x 10* q 


(2.59) 


The concentration of metallic atoms is not affected by the residence 
time of the metallic ions. If this residence time (T ) is taken to be 


12 hr, then the concentrations of metallic ions at sunrise and sunset are 
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[ M *] 1 = (2.19 x 10 4 )^ + (4.32 x lO 4 )^ (2.40) 

[W + ] 2 = (5.67 x 10 4 )^ + (4.32 x 10 4 )q J (2.41) 

For T~ l = 24 hr, 

= (6.30 x 10 4 )^ + (.8 . 64 x 10 4 )q r (2.42) 

[M*] 2 = (1.02 x 10 S )q A + (8.64 x lO 4 )^ (2.43) 


Curves showing the dependence of the diurnal variation in the 

metallic ion concentration on q ^ and q^ are plotted in Figure 2.4 for 

T -1 = 12 hr and in Figure 2.5 for T -1 = 24 hr. The observed diurnal variation 

of [M + ] and the value of matching the estimates of the preceding discussion 

-1 -3-1 -2 -3-1 

are approximated where q^ =10 cm s and q^ - 10 “ cm s . These 
values for q. and are consistent with calculations of the fraction of 
atoms ionized during ablation [Sida 3 1969; Lebedinets, 1973], 

The estimated values for q^. and q^ are used in equations (2.36) and 
(2.37) with values of [Af] j , [W] 2 » and [M *] 2 based on equations (2.38) 

through (2.43) and using values of T * of 12 and 24 hr, representing the 
uiurnal and semidiurnal tides. The resulting diurnal behavior of [M ] and 
[M] is plotted in Figure 2.6. 

In the model described above, curves similar to Figure 2.6 can be 

found using any combination of q-, q^ 3 and T. The basic diurnal trends 

are not dependent on the current estimates of the relative importance of 

the three processes since changes in their relative values would affect 

only the magnitude of the diurnal variation. If q^ >> q the diurnal 

+ + — 

variation would be negligible; while if q^ << q^, [M ]^/[M ] 2 = e , 

where t = 12 hr. Also, larger values of T result in larger diurnal 






Figure 2.5. Diurnal variation of the metallic ion concentration for various 

values of < 7 ^ and q ^ with a residence time due to transport of 24 hr 
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Figure 2.6. Time variation of metallic atom and ion concentrations 
shown for transport lifetimes of 12 hr. and 24 hr. 

The concentration of metallic atoms is not affected by 
the transport lifetime. 
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variations of [M + ], but do not affect [A/]. A model similar to this was 
developed by Gadsden [1970] to explain the observed distribution of sodium. 

2.3.6 Effect of heavy ions. Meteor trails probably contain a continuum 
of particle masses. In studies of twilight luminescence. Link [1973] found 
idence of meteoric aerosols at altitudes near 100 km. These metallic 
particles can also become ionized and would be acted upon by the neutral 
wind as described by wind-shear theory. 

The efficiency of the wind-shear mechanism in converging ionization 

increases as the mass of the ions increases. An estimate of the effectiveness 

of a wind shear in maintaining a density gradient against diffusion is given 

by the ratio v /D (see equation 2.22). 
z 

The effect of the mass of an ionized particle on the vertical drift 
induced by the motion of the neutral wind is found by equation (2.17). The 
ion-neutral collision frequency is proportional to M [Ratcliffe, 1972], 
where M is the mass of the ion. Thus the ratio of the ion-neutral collision 
frequency to the ion gyrofrequency can be expressed as 

p = (ca /) 1/2 (2.44) 

where C is a constant. 

The molecular diffusion coefficient is given by 

D = kT/m (2.45) 

so that, for an east-west horizontal wind, equations (2.17) and (2.45) give 
the dependence o' v /D on the mass of the ionized particle as 


v 

z 


D 


a 


M 


1 *CM 




(2.46) 


The dependence of v n /D on the mass of the particles is shown in Figure 

2.7. The altitudes indicated along the curve mark the location where, for 

the atmospheric conditions at that altitude, M is equal to an ion mass of 

30 amu. This figure shows that the effect of the wind-shear mechanism is 

independent of mass at altitudes below about 105 km. At an altitude of 

120 km, if the mass of the ions is increased by a factor of 10, the ratio 

v /D increases by a factor of 1.8, from 0.5 to 0.9 in the units shown in 
z 

Figure 2.7. This effect would be enhanced if the particles were more than 
singly charged. 


7 l-ffect of ionic mass on the efficiency of the wind- 
shear mechanism. Altitudes indicated along the curve 
mark the location where, for atmospheric conditions 
at that altitude, M is equal to an ion muss of 30 amu 
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5. REFLECTION OF RADIO WAVES BY SPORADIC-*/ LAYERS 

3.1 Introduction 

I.i this chapter, the observation of sporadic-*/ layers by reflection of 
radio waves is discussed. A full-wave analysis of the reflection coefficient 
is developed and applied to electron-density profiles measured by rocket 
probes. It is shown that partial reflections from the steep gradients at 
the upper and lower boundaries of sporadic-* 1 layers cannot account f. r 
the partial transparency. 

3.2 Ionosonae Observations of Sporadic E 

Most observations of sporadic-* 1 layers are made by the reflection of 
radio waves from free electrons in the layer. However, information on the 
small-scale structure of sporadic-* 1 layers gained by the reflection of radio 
waves is limited by poor horizontal resolution, and by the fact that only 
the part of the layer below the electron-density peak can be investigated. 
These limitations have thus far prevented the resolution of the source of 
partial transparency at frequencies between fbEs and ftEs. 

Vertical -incidence ionosphere sounders, or ionosondes, hace been used 
extensively in the study of the ionosphere. The first measurements using 
an ionosonde-1 ike device were made by Breit and Tuve [1926]. The technique 
developed rapidly after the initial experiments, and by the 1950's, 
ionosondes were being operated routinely in most parts of the world [Smith, 
1957]. 

Today there are several variations of the basic ionosonde. Cne employs 
three or more spaced receivers to record phase information in order to 
determine the position and velocity of the ionization at the point of 
reflection [ Wright , 1968; Wright and Fedor , 1969]. Scale sizes of iono- 
spheric irregularities are estimated from the diffraction patterns of the 
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reflected radiation [Beynon and Wright , 1969; Armstrong and Lue , 1970; 

Goodwin and Summers , 1970] . Another variation involves the separate 
detection of the ordinary and extraordinary reflections [ Whitehead , 1972a; 
Chessell et at., 1973]. 

Ionograms arc recorded routinely at many stations, usually at intervals 
of 15 min or 1 hr. Examples were shown in Figure 1.2. Various parameters, 
including fbEs and ftEs, are scaled according to international ly accepted 
guidelines [V'fJJ, 1972] . The parameters relating to sporadic E 
which are scaled from ionograms are sketched in an idealized ionogram in 
Figure 1.3 

The recording of ionograms at regular intervals has made possible 
statistical studies that have revealed many of the characteristics of 
sporadic E. World-wide maps of the frequency of occurrence of intense 
sporadic E show a strong maximum in the area of Japan and the Philippines, 
and a weaker minimum over South Africa [Smith, 1957; Leighton et al., 1962]. 
Diurnal and seasonal variations are apparent, with maximum valuer of fbEs 
and ftEs occurring at midday and during local summer [Thomas and Smith, 1959; 
Reader, 1962]. In an observation quite relevant to this study, Heisler and 
Whitehead [1962] found that ftEs increased before and decreased after corre- 
sponding changes in fbEs. They also noted that large differences between 
ftEs and fl.Es were most likely to occur during periods where either or both 
of the frequencies were changing. 

Ir chapter 1 two models were suggested as possible sources of the 
partial transparency of sporadic-c layers. In one the source is considered 
to be reflection or scattering from irregularities embedded within the 
layers; the other holds that gradient reflections from thin, horizontally- 
stratified layers are responsible. There is evidence supporting both 
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positions [Smith, 1968], 

Rawer [1962] concluded from a statistical study of sporadic-^ that ftEs 
represents the plasma frequency of the maximum electron density of the layer, 
while fbE8 is the plasma frequency of the smallest value of the peak electron 
density. He determined on the basis of the 1962 study and ea.lier calculations 
[Rawer, 1939] that partial reflections can be neglected in interpreting 
ionograms. The same conclusion was reached by Whitehead [1972a] through an 
evaluation of the reflection characteristics of the ordinary and extraordinary 
components of the reflected wave. Spaced ionosondes have shown evidence of 
irregularities with scale sizes varying from one to a few tens of kilometers 
[DeGregorio et al>, 1962; Armstrong and Lue, 1970; Goodwin cc\d Summers, 

1970]. GoocLrin and Summers interpreted the irregularities to be frontal 
irregularities in the form of ridges and troughs usually travelling to the 
north-west or south-east. However, these measurements may be affected 
by gravity-wave oscillations in the height of the layer focusing or 
defocusing the reflected radiation and giving the illusion of movement 
[Hines and Reddy, 1967]. 

Reddy and Rao [1968] observed that fbEs more closely matches the peak 
plasma frequency of sporadic-^ layers measured by rockets than does ftEs. 

They also observed that while fbEs is generally smaller, by about 10%, 
than the peak plasma frequency indicated by the rocket measurement, foEs 
is much larger. This behavior is characteristic of reflections from thin 
layers. Chessell [1971a] showed that the frequency region in which a thin 
layer is partially transparent is increased by magnetoionic coupling of the 
two characteristic modes of propagation at frequencies below the peak 
plasma frequency of the layer. 

The coherent scattering of radiation from small electron-density 
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irregularities has been considered by Tao and Wait [1962] and Reddy [1968] 
and will not be treated in this work. Bowhill [1966] showed that scattering 
can account for the partially transparent echo on some occasions. 

The status of the two models discussed above when this work was begun 
was that the validity of neither could be established definitively. As 
will be shown later, previous studies of reflections from sporadic-# layers 
had deficiencies which left the question open. The purpose of the following 
analysis is to determine conclusively whether gradient reflections from 
sporadic-# layers are of sufficient magnitude to be the source of the observed 
partial transparency of the layers. 

3.3 Reflection Coefficients 

3.3.1 Introduction. This section presents a discussion of the propa- 
gation of electromagnetic radiation through a magnetoionic medium. The 
basic equations governing propagation are developed in a form readily suited 
to the derivation of reflection and transmission coefficients. Approximate 
methods of derivation of reflection coefficients from ionized layers are 
discussed and shown to be inadequate for this study. However, the approximate 
expressions are used to check the validity of the full-wave analysis. 

A scheme is described whereby the full-wave reflection and transmission 
coefficients are programmed on a digital computer. This requires the 
description of sporadic-# electron-density profiles by mathematical functions. 
The results of this full-wave analysis are discussed in Section 3.4. 

3.3.2 Propagation of radio waves in the ionosphere. As most observations 
of sporadic # are made by experiments involving radio waves, it is appropriate 
that the effect that gradients in the ionosphere have on radio waves be 
examined. This development is an extension of the work of Budden [1961] and 
Pitteway [1965]. 
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The matrix differential equation describing a vert ica 1 ly-incident electro- 
magnetic wave in the ionosphere is [Budden , 1961] 


de 

-f- = -ikle 
dz 


(3.1) 


where k is the magnitude of the free space wave number. 


and 


e = 


x 


(y /e ) 2 H 

O O X 


(y /e 

K o o y 


(3.2) 


T = 


M M 

-M ♦ 

yx (l+A/ aJg ) 


1 + M 
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M M 
xz zy 


xy (1 +M ZZ ) 


0 1 


1 0 


0 0 


0 0 


(3.3) 


M is the susceptibility tensor of the ionosphere. 

By an appropriate choice of the orientation of the coordinate system, 

T can be written simply in terms of the actual components of M. The 
coordinate system will be oriented so that the geomagnetic field lies in the 
y-z plane, with z vertically upward and x in an eastward direction. The 
susceptibility tensor is written in this coordinate system as 
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M = - X 
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(3.4) 
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(3.5) 


qB uj 

Y = — E - — r 
mu co 





(3.6) 

and 

U = 1 - i v/ w 





(3.7) 

The 

transverse (Yy) and longitudinal (Y^) 

components 

of are 

given 


by Y and Y . r is a unit vector in the direction of the magnetic field. 

y 2 

q and m are the charge and mass of the electron, N is the electron number 
density, and v is the electron-neutral collision frequency. The plasma 


frequency of the medium and the electron gyrofrequency are, respectively, 
and u^. 

The solution of equation (3.1) provides expressions for the components 
of the electromagnetic fields in the ionosphere as a function of altitude 
for a given electron-density profile. This is called a '‘full -wave" solution, 
and involves the simultaneous solution of four coupled complex differential 
equations for each mode. A method of solution using a digital computer is 
described later in this chapter. Once the total electromagnetic field is 
known, it can be separated into upward and downward travelling waves from 
which the reflection coefficient can be calculated. 
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In order to derive reflection and transmission coefficients using the 
equations of this section, the polarization of an electromagnetic wave and 
the refractive inuex of the medium must be known. The polarization, p, of 
a progressive wave travelling in the +2 direction is defined as 

p = E v /E x (3.8) 

In terms of the magnetoionic parameters defined in equations (3.5), (3.6), 
and (3.7), the polarization of the wave in the ionosphere is 

- U T 

p = 2 Y l {U-X) + 1 

The refractive index of the ionized medium is given by 


2 2 
{u-xy 


+ 1 


(3.9) 


n = 1 
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11 r t 
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2 {U-X) ~ 

2 * T 

L4 (£/-*) J 



(3.10) 


It will be noticed that (3.9) and (3.10) are double-valued. In each 

case the upper sign applied to the ordinary mode of propagation and the 

lower sign to the extraordinary mode. The refractive indices will be 

2 

designated n and n ^ respectively. The real part of n is plotted in 
Figure 3.1 for conditions at an altitude of 100 km. The values used are 


v/2-t = 0.02 MHz and uJ2v = 1.4 MHz. 

H 

Figure 3.1 shows that there are values of X and corresponding electron 
densities at which the refractive index vanishes. The effective wavenumber 
( kn ) also vanishes and the wave is reflected. In the absence of collisions, 
these points occur at X = 1 for the ordinary mode, and X = 1 * Y for the 
extraordinary mode. At greater electron densities the refractive index is 


0 2 4 6 8 10 1.2 1.4 1.6 1.8 20 

X 

2 

Figure 3.] Variation of the real part of n with .V, where n is the refractive 

index of a magnetoionic medium and X is the square of the ratio of £ 

the plasma frequency of the medium to the incident frequency of an 
electromagnetic wave. 
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imaginary, and the wave is evanescent. Figure 3.2 shows a representative 
electron-density profile for a sporadic-# layer. Figure 3.3 shows, for this 
layer, the altitudes at which the index of refraction vanishes for the modes 
of propagation. 

l.y 2 

There is a singularity in the refractive index at X = =-, which is 

just less than X = 1. In this region, the change in the elliptical polari- 
zation of the wave is very rapid, resulting in strong magnetoionic coupling 
between the modes [Budden, 1961]. The location of this singularity is 
indicated in Figure 3.3 by a dashed line. In the region between the infinity 
in and the zero in n both refractive indices are real. A wave with 
extraordinary polarization resulting from magnetoionic coupling in this region 
would propagate past the point where X = 1 and, depending on its frequency, 
reflect where X = 1 + Y or propagate through the layer. The so-called Z-trace, 
seen at times on ionograms, is the result of this reflection. 

3.3.3 Approximate solutions. Several approximate expressions for reflec- 
tion coefficient are used in optics and electrodynamics. Three of these 
expressions provide useful checks on the accuracy of the full -wave solutions. 
Derivations of the three approximate solutions are presented in Appendix II. 

The most direct approach to an estimation of the reflection coefficient 
is to approximate the ionized layer by a series of vertically thin layers 
(Section II. 1). The electron density is assumed constant within each layer. 


Denoting the refractive index in the regions below and above the series 
of layers respectively as n. and n , the reflection and transmission 
coefficients are given by, from equations (11.15) and (11.16), 



A Z (Km) 



Z(km) 



Figure 3.5 Altitudes at which the index of refraction vanishes 
for the nodes of propagation, for the case of the 
electron-density profile of Figure 3 . J . 
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an. * bn .n , + c + dn, 
z t t t 


where a, b, c, and d are constants that depend on the thickness and 
refractive indices of the layers in the reflecting medium. In the case of 
a single layer these become the Fresnel formulas. 

These equations are invalid where the frequency of the wave is at or 
below the plasma frequency of the medium. Also, reflections generated by 
the discontinuities at the boundaries between levels make it impossible to 
examine small reflection coefficients from a continuous medium using the 
approximation. 

Another solution, valid in a slowly varying medium where the frequency 
of radiation does not approach the plasma frequency, is the W.K.B. approximation 
(Section 1 1. 2). In this method, the wave equation 


2 2 

+ k n E = 0 


(3.13) 


is solved by substituting an arbitrary function of height for knz , where k 
is the wave number of the vertically-propagating wave, and n is the height - 
dependent refractive index. The refractive index is assumed to vary slowly 
with height, and several approximations are possible by virtue of its small 
secon derivative. The reflection coefficient for an ionized layer is 
found to be, equation (11.32) 


r _U An f z 

R = -'■$ n ^ ~dz ex P rt(a') 

-oo n 


dz ' dz 


(3.14) 


This approximation is valid in sporadic-^ layers unless the frequency of 
the incident radiation approaches the plasma frequency of the medium. However, 
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it is unacceptable for the present study because of the desire to examine 
the reflection coefficient at frequencies near and below the peak plasma 
frequency of the layer. 

The third approximate solution assumes the ionisation profile to be in 
the shape of a so-called Epstein profile [Budden, 1961], For this case, 
exact expressions for the reflection coefficient can be derived using solutions 
of the hypergeometric equation (Section 1 1. 3). 

For a symmetric layer of ionization, the magnitude of the reflection 
coefficient is found to be, equation (11.85) 


fl = 


-i *1 


cos (ttJ) 

2 2 2 2 
cos (tt<2) cosh (tt5) + sin (it d) sinh (ttS) 


(3.15) 


where 
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7~. 

J incident 


and m is a dimensionless parameter related to the full-width at half-maximum 
(As) of the layer by Aa = 3.526/m meters. 

The reflection coefficient of a transition region between low and high 
electron density is, equation (11.82) 


1*1 


sinh [~( !-*„,) ] 
sinh[^(l*nj] 


(3.16) 


Examples of symmetric and transition Epstein layers are shown in Figure 3.4. 


0 .1 .2 .3 4 .5 

ELECTRON DENSITY (I0 5 crr* 3 ) 


Figure 3.4 S>mjnetric and transitional Epstein profiles. 
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Reflection coefficients calculated by these approximate methods will 
later be compared with full-wave solutions in cases where the approximations 
are valid. 

3.3.4 Full-wave solution. The mathematical development of the equations 
describing the propagation of electromagnetic waves in an ionized medium 
were discussed in Section 3.2.1. It is the purpose of this section to 
derive reflection and transmission coefficients that accurately described 
the actual reflection and transmission by a sporadic-E layer. Because of 
coupling between modes of the incident radiation by the geomagnetic field, 
some transmission through the layer at frequencies below the peak plasma 
frequency is expected. Therefore, a full-wave analysis is necessary. The 
results of this analysis will be checked against approximate methods and 
also against independent derivations of full-wave solutions. 

The full-wave solution of the reflection and transmission coefficients 
builds directly on the derivation in Section 3.3.1, where it was shown that 
the electromagnetic fields in the wave and their derivatives are related by 

de 

= -ikTe (3.1) 


This equation is integrated through the region of ionization for each of two 
independent modes of e. The initial conditions for each mode specify an 
ordinary or extraordinary upward travelling wave above the layer of ionization. 
Thus, from Maxwell's equations and equation (3.8), 
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so that 

e - 


The polarization of the wave is specified in the expressions used for p and 
n, equations (3.9) and (3.10). 

By integrating equation (3.1) using the ordinary wave polarization in 
equation (3.18) as the initial condition, and again, using the extraordinary 
wave polarization, two independent solutions, and e^, are formed 

below the layer. Each solution can be separated into upward and downward 
travelling parts. Then, by taking the appropriate ratios, the reflection and 
transmission coefficients can be derived. 

The electric field can be expressed as the superposition of upward (£/) 
and downward (£>) travelling waves in the following manner: 

i - 1/* 0 * . t/ kn * = (3.i9) 


C V o U 0> 
(u c /c o^ " 


(3.18) 


from which 


j— = -iknd + iknfi 
dz 


(3.20) 


These two equations can be solved for the upward and downward travelling 


waves , 
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(3.21) 

(3.22) 


i 


I 


4 


55 


In general, the ordinary and extraordinary waves are elliptically 
polarized. The angle formed by their major axes in the x-y plane is bisected 
by a line rotated clockwise 45° from magnetic north [ Budden 3 1961], Thus, 
in order to compare electric fields of the two modes with each other the x 
component of one must be compared with the y component of the other [Chessell, 
1971b]. 

Using the definition of polarization, equation (3.8), 
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X ox *x 
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(5.23) 


where the subscripts o and x refer to the ordinary and extraordinary modes. 
The simultaneous solutions of equations (3.23) are 
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Then, by substituting equations (3.24) into (3.21) and (3.22), the 
desired components of the incident and reflected waves can be written 
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The transmitted waves have been specified by the initial conditions of 
the integration, equation (3.18), to be purely ordinary and purely extra- 
ordinary upward travelling waves, and will be represented by and V 
The reflection and transmission coefficients are determined by 


(3.26) 
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where, for example, fl describes the reflection of incident extraordinary 

o X 

radiation as ordinary radiation. The superscripts indicate the independent 
solutions of equation (3.1). The solution of equations (3.26) and (3.27) give 
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The expressions describing the reflection and transmission of radio 
waves by a magnetoionic medium are now in a form whereby they can be solved 
by computer techniques. This solution, described in the following section, 
involves the use of a Runge-Kutta integration routine to solve equation 
(3.1) for components of E and H. These components can be readily transformed 
using equations (3.28) and (3.29) into reflection and transmission coefficients 
of circularly-polarized radiation. 

3.3.5 Solution using a digital computer. The full-wave solution of 
the wave equation in an ionized medium requires the integration of four 
complex differential equations, equation (3.1), through the region of ioni- 
zation. A Runge-Kutta 4th-order integration routine developed by Gill [1951] 
provides a relatively simple but accurate way to perform the integration. 

The derivation of the Runge-Kutta algorithm can be found in several 
texts [ Stanton , 1961; Antosiewicz and Gaul .chi , 1962]. At each step in the 
integration, the value of the function is found from the previous value 
using a Taylor expansion 

,2 y 3 

e(;:+ft) = e(z) + he' ( 2 ) + yj-e ,, ( 3 ) + 37 e ' ' ' ( 3 ) + .... (3.30) 
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This is applied to each of the four differential equations (3.1), expressed 


as 


Tz * 


(3.31) 


with the result that, keeping terms up to fourth order in h, 
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in which the subscripts denote differentiation. This equation is to be 
expressed in the form 


A e = aK + bK_ + oK + dK 
u) 2 3 4 


(3.33) 


where 


(3.34) 


= hf(z s e) 

K 2 = hf(z+nih, e+nX^) 

= hf(z+ph, e+qK^+rK^) 

Kq = hf^z+sh, e+tK^+uK^+vK^) 

To do this, equations (3.34) are expanded about z and e again using a Taylor 
series. To simplify the resulting series, define 
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where a = 1, 2, 3; and 
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Using these definitions, (3.34) becomes 
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Now, equations (3.36) are substituted into (3.33), which is equated with 
equation (3.32). By matching coefficients of like powers of h, one can write 
thirteen independent equations and solve for the thirteen coefficients in 
equations (3.33) and (3.34). Thus, 
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and 


A'j = hf(z,y) 

A -> = hf(z + hh, y + Vfj) 

(3.38) 

Aj = hf(z * \h 3 y - ('* ~ 2~' 2 )X l + U - 2'V 2 
A 4 = hf[z + h, y - 2~ 2 * (1 + 2" 2 )A' 3 

The value of each component of e is given at each step by 

e(z) = e(s-A 2 ) + A e(z) (3.39) 

The above set of equations is in a form which is easily adapted for use 
on a digital computer. In this study, the version contained in the IBM 
Scientific Subroutines Package was used. This routine allows the simultaneous 
solution for both independent modes. Also, each step in the integration is 
done for two different step sizes, thus testing the accuracy. If the error 
excedes a given value, the step size is halved. This is repeated until the 
desired accuracy is attained. 

3.3.6 Description of a sporadic -E layer bu a mathematical function. 

The electron-density profiles of sporadic-^ layers have been measured by 
Langmuir probes carried by sounding rockets [Smithy 1969]. This technique 
provides a continuous measure of current to the probe. Over a limited height 
range the current is proportional to the electron density. The nature of 
digital computers requires that the electron density be either stored as a 
discrete array or generated at each step. Because of the uncertainty 
associated with the step size and the possibility of fictitious reflections 
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being generated by discontinuities in the slope, the latter option was taken. 
This meant that each layer to be studied first had to be described by a 
mathematical function. 

Functions describing the layers were found using a least squares technique 
which optimized the fit to a given series of functions by varying che parameters 
within these functions [Bevington, 1969]. A trial function was constructed 
for a layer, and the best fit of that particular function to the data was 
found. This function would then be subtracted from the data giving a new 
data set. A new function was then fitted to the new data set and the difference 
determined. This process was repeated until a good fit was obtained. The 
final electron-density profile would be the sum of the functions that best 
fit the original data. For example. Figure 3.5 shows the set of functions 
used to describe a sporadic-i’ layer observed at Wallops Island on 16 July 1970. 
This profile provides a good test for the method of curve fitting, as it is 
more complex than most sporadic-/? layers. It was possible by this method to 
model each of the sporadic-£ profiles used in this study to a high degree of 
accuracy. 

3.3.7 Performing the integration. The solution of equation (3.1) is 
found for various sporadic-£ profiles using the Runge-Kutta method of 
integration through the layer. The starting values for two independent modes 
are presented in equation (3.18), and describe upward-travelling ordinary 
and extraordinary waves. The integration is started and ended in regions of 
the profile where the gradient of electron density is small, thus avoiding 
the introduction of apparent reflections from discontinuities in the gradient. 
Double precision variables are employed throughout the computer program. 

Pittewau [1965] showed analytically that the independence of the two 
solutions should be preserved throughout the integration. It is found that, 
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during the course of the integration, the two modes, which were originally 

independent, lose their independence. The loss of their independence can be 

attributed to the finite accuracy of the digital computer. The truncation of 

the components of the solutions has the effect of introducing elements of one 

solution into the other, with the effect that as the integration proceeds, the 

two solutions will be dominated by the most rapidly growing eigenmode. 

Following the suggestion of Pitteway [1965], a constraint is introduced 

at regular intervals through the integration which insures that the two 

solutions are Hermitian orthogonal. This is done by replacing the second 
( 2 ) ( 2 ') 

solution, e , by a new solution, e , which is made up of a linear 
combination of the two solutions, but defined so as tc be Hermitian orthogonal 
to the first solution, . Thus, 
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(3.40) 


The rapid growth of evanescent modes was limited by the following 
procedure: if the value of any component of either solution became larger 
than ten, all components of both solutions and their derivatives were reduced 
by a factor of ten. The number of times this operation was done was recorded 
for subsequent use. 

Neither the orthogonality constraint nor the reduction of the size of 
the components and derivatives affects the reflection coefficient since for 
this calculation upward and downward travelling waves are compared with each 
other at the same level. All that is required is that there are two 
independent solutions. 

For transmission coefficients, upward travelling waves below the layer 
are compared with upward travelling waves above the layer. Because of the 
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addition of solutions to preserve orthogonality, information of transmission 
resulting from mode coupling is not preserved. Therefore, upgoing ordinary 
and extraordinary waves below the layer are compared respectively with 
ordinary and extraordinary waves above the layer to get and The 

waves below the layer are multiplied by the appropriate factors of ten to 
compensate for the procedure used to limit the evanescent modes. 

3.4 Results of Reflection Coefficient Calculations 

3.4.1 Comparison with independent solutions. The validity of the 
derived reflection and transmission coefficients discussed above and the 
accuracy of the computer program employed are now tested against some of the 
approximate solutions of Section 3.3.2. In each case the full -wave solution 
is found to be consistent with the result of the approximate method. The 
solutions are also compared with full-wave solutions published in the literature. 

As an initial test, the analysis was made using a constant electron 
density and an incident frequency of radiation well above the plasma frequency 
of the medium. To the accuracy of the computer, the result showed total 
transmission with no reflection. The electron-density profile was then 
modified so as to describe a semi-infinite region of constant electron density 
bounded below by free space. The reflection and transmission coefficients of 
a wave incident on this boundary, as given by the full -wave solution, were 
equal to the values given by equations (3.11) and (3,12). 

The next approximate solution was derived using the W.K.B. method. 

The reflection coefficient was calculated from equation (3.14) and compared 
with the full-wave solution. A symmetric Epstein profile used in the 
comparison is shown in Figure 3.4, and is described by 
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where X is the value of X at the peak of the layer and m = 0.01. The 
full-width at half -maxi mum of this Epstein profile is 352.6 meters. The 
results of this calculation, showing good agreement between the two methods, 
are plotted in Figure 3.6. 

The symmetric electron-density profile of Figure 3.4 is also used to 
test the full-wave solution against the solution derived from the hypergeo- 
metric equation. In this case, the reflection coefficient is given by 
equation (3.15). As is shewn in Figure 3.7, there is good agreement between 
the two methods. The result is also shown for a thicker profile with m = 0.00 
The reflection coefficients of transition layers with m = 0.01 and m = 0.02 
are shown in Figure 3.8. The difference between the value of the reflection 
coefficient from the W.K.B. method ?"d from the hypergeometric function is the 
result of neglecting the geomagnetic field in the latter case. 

Chessell [1971b] has calculated reflection and transmission coefficients 
of several types of profiles using a full wave analysis. Chessell's results 
were used to check the full-wave solution used in this study. This comparison 
is shown in Figure 3.9. 

The electron-density profile used in this comparison with Chessell's work 
is a parabola with a half-width of 477 m, measured to the point where the 
electron density is zero. This profile can be used effectively to study stron 
reflections from sporadic-F layers. However, the discontinuities in electron- 
density gradient at the upper and lower boundaries of the layer give rise to 
reflections that mask any gradient reflections. 

3.4.2 The effect of ionosonde sensitivity. The results of this study 
are to be compared with ionograms recorded at Wallops Island at the time the 
profiles were measured by rocket-borne probes. In this connection it is 
necessary to discuss the effect of ionosonde sensitivity on the ionograms. 
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Figure 3.9 Comparison of full-wave solutions of reflection and transmission 

coefficients developed in this study with full -wave solutions from 
Chessell [1971b] ("solid lines). The layer used was a parabola 
with a half-width of 477 m. 
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As the signal to an ionosonde becomes weaker, the gain is automatically 
increased. Also, the standard methods of recording ionograms by photographing 
oscilloscope screens or burning paper with a spark convey limited amplitude 
information. The effect is that a weak echo appears to be no different from 
a strong echo on the ionogram. 

The sensitivity of the ionosonde at Wallops Island was measured by 

attenuating the signals received by the ionosonde. The attenuation was 

increased until the trace disappeared from the ionogram. These tests showed 

that a layer at sporadic-# altitude with a reflection ceofficient greater 

than 10 ^ would be recorded on the ionogram. Other researchers estimate this 

-3 _ -> 

threshold reflection coefficient to be between 10 and 10 [Smith, 1957; 

Reddy , 1968]. 

Because of this automatic adjustment of sensitivity the partial 
transparency indicated in Figure 3.9, for example, at frequencies below 
the plasma frequency is not observed. This partial transparency in the 

frequency range f = f to f = f - results from magnetoionic mode 

P P Hf .. 

coupling (as discussed at the end of Section 3.3.2). Based on the calculations 
of this study and assuming similar transmission coefficients for waves 
travelling either direction through the layer, it appears likely that reflec- 
tions from layers above the sporadic-# layer would be recorded on the ionogram. 

For example, a signal reflecting totally at an altitude of 300 km which twice 

_2 

penetrates a layer with a transmission coefficient of 9.5 x 10 ‘ would have 
the same amplitude at the ground as if it had reflected at an altitude of 
100 km from a layer with a reflection ceofficient of 10 ^ . It is not likely, 


however, that this echo from the # region would be detected by an ionosonde. 
At frequencies just below the peak plasma frequency there would be a strong 
return from the sporadic-# layer, and the automatic gain control would set 
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the gain of the ionosonde at a low value so that the F - region echo would not 
be received. 

As noted in the previous section care has to be taken in the full -wave 
analysis to avoid discontinuities in the electron-density gradient. The 
apparent reflection coefficient from a discontinuous gradient will easily 

_3 

exceed 10 . If it were real, this reflection would be recorded on an 

ionogram. The problem of apparent reflections from discontinuous electron- 
density gradients has often been overlooked in the past in studies of 
reflection coefficients. For example, the conclusion of Reddy [1968] that 
gradient reflection is an important source of partial reflections from 
sporadic-# layers is based on apparent reflections from the two gradient 
discontinuities in the electron-density profile used. 

When actual electron-density profiles are used, it is usually impossible 
to begin and end the integration through the layer at altitudes where the 
gradient iz zero. In this study, the end points were carefully chosen to 
minimize this effect. Reflections from the two discontinuities interfere 
destructively and constructively. The result is that the reflection 
coefficient becomes periodic as a function of incident frequency, with the 
maxima giving a measure of the magnitude of the error introduced into the 
reflection coefficient. This periodicity makes for easy recognition, as can 
be seen in Figures 3.6, 3.7, and 3.8. 

3.4.3 Results of reflection coefficient only sis. A number of 
sporadic-# profiles have been measured by rocket probes. It is estimated 
that, excluding flights with the specific purpose of studying sporadic E, 

50*. of the daylight rocket launches and virtually all of the nighttime 
launches in the Aeronomy Laboratory program nave shown sporadic-# layers 
[Smith and Mechtly , 1972]. Four daytime profiles and one nighttime profile 
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have been selected for analysis of reflection coefficients. 

The first daytime profile is shown in figure 3.10 together with the 
fitted curve used in the full-wave analysis. The circumstances of the layer 
are described by Smith and Mechtly [1972]. Ionograms recorded during the 
launch at noon (EST) on 24 July 1968, showed fbEs to lie 6.70 MHz and ftEa 
to be 7.65 MHz. 

The reflection and transmission coefficients calculated from this profile 
are plotted in Figure 3.11. N'o echoes attributable to gradient rcfle«_L ons 
would be expected on the ionogram. The difference in frequency of fbEs and 
ftEs in this case (0.94 MHz) is accounted for by the different reflection 
frequencies of the ordinary and extraordinary waves. The calculated difference, 
based on the difference between the frequency at which becomes less than 

-3 

10 and the peak plasma frequency of the layer, is 0.90 MHz. 

The second daytime profile is shown in Figure 3.12, again with the 
fitted curve used in the analysis. For this layer, at 11:00 LSI on 16 July 
1970, fbEs was 5.35 MHz and ftEs was 10.70 MHz, a difference of 5.55 MHz, 
much greater than in the previous case. The calculated reflection and trans- 
mission coefficients are shown in Figure 3.13. The blanketing frequency agrees 
well with the reflection of the ordinary wave. The crid'-.nt of the lower 
boundary of this layer is rather steeper than the previous case and, indeed, 
scronger gradient reflections can be noted. At |P| = 10 * , the limit of 
sensitivity of the ionosonde, the expected difference in f: Es and ftEs is 
1.39 MHz, greater than for the previous example, but not large enough to 
account for the observations. 

Sporadic-^ layers are almost always present at night. 'V lij [1968] 
reports a tendency for the difference between f! and f: :: to be greater 
at night than in the day. In most instances nighttime sporadic-E layers arc 
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Figure 3.10 tlectron-density profile of sporadic-!-, layer observed 

at Wallops Island, 24 July 1068. The fitted curve used 
in the analysis is the dashed profile. 


% 


i 

I 


| 

i 








ALTITUDE (km) 



NIKE CAJUN 10.320 
1 1 00 EST 
16 JULY 1970 
WALLOPS ISLAND 



ELECTRON DENSITY (I0 5 cm' 3 ) 


i 


I 


I 


Figure 3.12 F lectron-dens i tv profile of sporadic-.- layer observed at 

Wallops Island, 16 July l!>"n, with the fitted curve used 
in the analysis. 
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thinner than daytime layers and have a lower peak electron density. iTie 
electron-density profile of a particularly intense nighttime layer is shown 
in Figure 3.14, again with the mathematical function used in the analysis. 

Its reflection and transmission coefficients are plotted in Figure 3.15. In 
this case fbEs was 1.7 Ml I z and ftEs was 4.6 MHz. Again, the top frequency 
of reflection cannot be explained on the basis of gradient reflections, the 
calculated value of ftEs - fbEs being 1.07 MHz. 

Two rocket launches have been made during this study of the partial 
transparency of sporadic-£ layers. The first was launched on 3 August 19"3 
into an intense sporadic-F layer when the range of partial transparency was 
small: fbEs was 5.5 MHz and ftEs was 6.3 MHz. At the time of the second 
launch, on 10 August 1973, partial transparency was detected over a large 
range of frequencies: ftjEs was 4.0 MHz and ftEs was 8.0 Ifllz. Ionograms 

recorded at the times of the two launches were shown in Figure 1.2. 

The electron-density profiles of these two layers are plotted together 
in Figure 3.16 together with the fitted curves. Initially the difference 
between the electron-density gradients in the two profiles encouraged the 
idea that the partial transparency results from gradient reflection. However, 
calculation of the reflection coefficients of the layers show that this is not 
the case. The reflection and transmission coefficients of the August 3 layer 
are shown in Figure 3.17, and of the August 10 layer in Figure 3.18. They 
appear to be very similar. In neither case can either ;7 rs or he 

explained on the basis of reflection coefficients. In the case of the layer 
observed on 3 August 1973, the value of ftEs - flEs based on the full-wave 
analysis is 0.95 MHz; while, for the 10 August 19-3, it is 0.9" MHz. 

The lack of exact coincidence between ftEs and the peak plasma frequency 
of the layers studies here is probably the result of the measurements being 
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Figure 3.1b Electron— density profiles of two sporaJie-h layers observed 
at Wallops Island, 3 August 11173 and 10 August 1973, with 
the fitted curves used in the analysis. Altitude scale is 
relative to the center of each layer. 
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Figure 3.17 Reflection and transmission coefficients of sporadic-E 
layer, 3 August 1973. 
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made at different locations. The F region echo from which fbEs is determined 
is reflected in the first Fresnel zone at the F region. At an altitude of 
300 km, the first Fresnel zone has a diameter of 8.5 km. Thus the radiation 
reflecting from the F region and detected at the ionosonde illuminates a 
circular section of the sporadic-!? layer which is only approximately 2.8 km 
in diameter. The conditions in this small region may be quite different from 
the conditions along the path of the rocket, which is typically separated by 
a distance on the order of 100 km. Smith and Meohtly [1972] also note cases 
in which simultaneous comparisons give disagreement, yet which can be 
reconciled by geographic considerations. 

The reflection coefficients of sporadic-# layers described in this section 
present convincing evidence that the partial transparency over a range of 
frequency greater than 1 MHz cannot be explained on the basis of gradient 
reflection. Although there is a small dependence on electron-density gradient, 
it is insufficient to account for the characteristics of the reflection from 


the layers studied. 
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4. OBSERVATION OF SPORADIC-E LAYERS BY INCOHERENT-SCATTER RADAR 
4 . 1 Introduction 

The previous chapter dealt with the possibility of gradient reflections 
being observed from sporadic-E layers. It was shown that the range of fre- 
quencies at which observed gradients of sporadic-E layers partially reflect 
radio waves is too small to account for the large differences between fbEs 
and ftEs that have been observed. 

According to a second model, partially transparent sporadic-E echoes 
are the result of reflections from regions of high electron density 
scattered throughout the horizontal plane of the layer. Measurements by 
radio reflection have been interpreted to show the existence of irregular 
structure in sporadic-E layers, but the implied structure has never been 
directly observed. 

Recent developments in incoherent-scatter radar have led to an in- 
creased resolution in range which now makes it possible to make effective 
measurements of sporadic-E layers. With the use of pulse -compress ion 
techniques, a resolution of 600 m in range is now being achieved; the 
thickness of a sporadic-E layer is generally between one and three kilo- 
meters . 

In this chapter, experiments will be described which exploit the technique 
of incoherent-scatter radar in the investigation of sporadic-E layers. The 
observations were taken during January, 1974, and July, 1975, at the Arecibo 
Observatory located in Puerto Rico. In one type of experiment, pulse 
compression techniques were used to gain the best possible space and time 
resolution of sporadic-E ionization. It will be shown that small regions 
of enhanced electron density do indeed exist at times in sporadic-E layers, 
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providing an explanation of the partial transparency of the sporadic-# 
echoes observed simultaneously by the local ionosonde . 

In the other experiment ion drift velocities were measured by a multiple- 
pulse technique. These data, in addition to indicating the motion of the 
layer, also allow investigation of the mechanism of formation of the layer. 

The main features of the observations are presented in this chapter. 

Their interpretation in terms of current theories of sporadic-# layers is 
considered in the subsequent chapter. 

A preliminary account of the experiments of January 1974 has been 
given by Miller and Smith [1975] . The material of that paper is included, 
properly identified, in this chapter. 

4 . 2 Incoherent Scatter 

4.2.1 Introduction. The development of high -power radar in the past 
three decades has enabled the scattering of radio waves from the ionosphere 
to become a powerful tool in the study of the atmosphere. The technique 
known as incoherent scatter provides the ability to measure continuously 
such atmospheric quantities as electron density, electron and ion 
temperatures, ion composition, and ion velocities. 

The rapid advance of radar technology during and immediately following 
World War II lead to the active study of scattering of radio waves by the 
lower ionosphere. The first efforts had the purpose of either reducing 
interference on military radars or developing long-range communication 
systems [ Booker and Gordon, 1950; Bowles, 1964] . Gordon [1958] proposed 
that radar technology was sufficiently advanced that incoherent backscatter 
could be detected from free electrons in the ionosphere. The first 
incoherent-scatter measurement was made in the same year by Bowles [1958] 
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at Long Branch, Illinois, using a 41-MHz radar. A review of the present 
status of incoherent scatter is presented by Evans [1975] . 

4.2.2 Scattering of radio waves. Incoherent scatter is used at fre- 
quencies much greater tha^ the plasma frequency of the ionosphere. The 
radio waves are scattered by variations in the permittivity caused by random 
fluctuations in the electron density. The scattered power received at the 
ground from incoherent scatter is shown by Evans [1969] to be 



where R is the range from the receiver to the region where the radiation 
is scattered, n is the electron density of the scattering medium, and C is 
a constant which is dependent on various parameters associated with the radar 
itself. The scattering cross section, o, is given by 

a 

a = 5 5 “ (4.2) 

(i+«)[HVr { )V] 

where a = 4 it D/\, D is the Debye length, and the ratio of electron tempera- 
ture to ion temperature, T /T ., is less than about 3. The radar cross 

2 

section of an electron, a , is equal to 4itr , where r is the classical 
electron radius. In practice, the power received is usually separated into 
range gates of length corresponding to the pulse length, and converted to 
electron density by calibrating it with the peak electron density of the 
F region, as measured by an ionosonde. 

The power received also contains an undesired noise component, P^. The 
signal-to-noise ratio, PJP„, is only of the order of unity in most inco- 

U If 

herent-scatter experiments. To improve the accuracy of the power measure- 
ment, a series of measurements are integrated [Evans and Loewen dial, 1964], 
The noise power can be estimated by measuring the power at a number (m) of 
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time delays from which no incoherent-scatter echo is expected, 

i m 

P„ = - l P. 

N m . L , v 


(4.3) 


This is then subtracted from the power measured at time delays where signal 
is expected, giving P . The accuracy of the measurement integrated over n 
sweeps of the radar time base is given by 


A P s = iUAP ) 2 *^) 2 ]* 5 


where 


A P = ± 


P +P 
S N 


n 


a p = ± 
a N 


’ N 




(4.4) 


(4.5) 


(4.6) 


Information can be deduced on electron and ion temperatures and ion 
composition from the shape of the frequency spectrum of scattered radiation 
[Ratcliffe, 1972; Evans, 1969]. The scattering of radiation is strongly 
controlled by fluctuations with wave number equal to twice the wave number 
of the incident radiation. Two types of longitudinal waves are mainly 
responsible for the scattering: ion-acoustic waves and electron plasma 

waves. Features in the spectrum of the returned signal can be identified 
with each type. 

For the incident wavelengths much longer than the Debye length used 
in incoherent-scatter measurements the contribution of the ion-acoustic 
waves is centered at the incident frequency, with maxima Doppler shifted 
to greater and lower frequencies by the mean speeds of the ions. An 
ion-acoustic wave of length A has a frequency of [Ratcliffe , 1972] 
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(4.7) 
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For thermal ion velocities of 
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(4.8) 


the Doppler shift of the backscattered signal is 

% 




y.iwyr.)] 


m . 


(4.9) 


where k n is Boltzmann's constant and m. is the ionic mass. 

B t 

The component in the incoherent-scatter spectrum due to electron plasma 
waves, has its frequency shifted by the plasma frequency, regardless of the 
incident wavelength. This component is called the plasma line, and has 
been used in measurements of electron density. 

An example of the scattered-power spectrum, adapted from Hagfors [1961], 
is shown in Figure 4.1. The shape of the spectrum is determined by the 
quantities T., T /T., and m. [Evans, 1969]. By assuming a value for one 
of the quantities, the other two can be determined by comparing the shape 
of the measured spectrum with theoretical spectra. 

Hie velocity of the ionization, if there is no net electric current, 
shifts the entire spectrum by 


2v, 

Lf d = ~r 


(4.10) 


where v ^ is the component of the ionization velocity toward the radar. If 
there is motion of the electrons with respect to the ions, the spectrum 
becomes asymmetric [Evans, 1969]. This effect is important in equatorial 
and aur'fal regions, but is not significant at midlatitudes. 

4.2.3 Autocorrelation analysis. In incoherent-scatter observations, 
it is more often the practice to measure the autocorrelation function of 
the signal than the frequency spectrum [Farley, 1969]. Autocorrelation 


i 

l 


I 


POWER W(Af e )/Nr e 2 




90 


* 


functions are measured Ly transmitting a series of pulses, either as a pair 
separated by a varying multiple of a time delay, x, or as a sequence of 
unevenly-spaced pulses from which different combinations give different 
time delays. The autocorrelation function is then constructed as 




t)> 


<|IW)| > 


(4.11) 


where V{z,t) is the complex signal (i.e. amplitude and phase information) 
received from a range z at time t. 

As in the case of the frequency spectrum, the shape of the autocorrela- 
tion function contains information on the several ionospheric properties 
(Figure 4.2). The time delay for the first zero-crossing in the autocorre- 
lation function is less for higher ion temperatures or lighter ions. The 
shape of the curve is affected by the electron-ion temperature ratio, ibn- 
neutral collisions, and the magnetic field. If there is either a net 
plasma drift or an electric current, the frequency spectrum is asymmetric 
and the autocorrelation function is complex. 

After the autocorrelation function is measured, it is compared with 
theoretical autocorrelation functions generated from various combinations 
of physical quantities [ Farley , 1969; Zamlutti and Farley, 1975]. 

The main advantage in measuring the autocorrelation function rather 
than the frequency spectrum is that the required accuracy can be attained 
using a shorter pulse. The increased height resolution resulting from 
the shorter pulse is gained at the expense of time resolution. Since the 
frequency spectrum can be obtained by a Fourier transform of the autocorre- 
lation function, information of atmospheric properties is not lost by 
the use of this technique. 



AUTOCORRELATION 


S 0 
< 


TIME DELAY (^s) 


Figure 4.2 Real and imaginary autocorrelation functions measured using 
the five-pulse technique with a 30-min integration time and 
a 16 us pulse. The solid curves show a fitted theoretical 
autocorrelation function corresponding to = 250 K, 
v.£ n = 300 s"l, and a radial drift velocity of 14 m s _ l 
away from the radar [Zarnlutti and Farley, 1975]. 
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4.2.4 Pulae oompreaaion. Another method of increasing height 
resolution without increasing transmitted power is known as pulse 
compression. Although this technique has been used in radio astronomy 
for several years, it has only recently been used successfully in incoherent 
scatter measurement [Ioannidis and Farley, 1972]. It is through the 
development of pulse- compression capabilities at the Arecibo Observatory 
that the height resolution was increased to the point where sporadic-# 
layers could be effectively studied. 

Using pulse compression, segments of a transmitted pulse are reversed 
in phase in such a way that the resulting frequency spectrum resembles that 
of a much shorter pulse. When the scattered pulse is received and decoded 
it is approximately the same as would be received from a short pulse with 
the energy of the complete long pulse that was transmitted [ Gray and Farley , 
1973]. 

At Arecibo a 13-bit Barker-coded pulse is used, shown with its 
autocorrelation function in Figure 4.3. The autocorrelation functions of 
Barker codes are characterized by uniform side lobes and a central maximum 
with an amplitude n times greater than the amplitude of the side lobes, 
where n is the number of bits in the code. 

A phase-coded pulse, b[t), after being scattered by a thin slab in 
the ionosphere and received at the ground, can be represented by b[t)f[t) 
[Gray and Farley, 1973]. When the signal is received it is cross -correlated 
with a replica of the transmitted pulse, so that the final signal can be 
written as 

G(t) = (4.12) 
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Then the ensemble average of the power is 


<|C(t)| 2 > 


r 


r 

.00 


(4.13) 


If t ^~ t is substituted for t 7 , becomes the autocorrelation 

function of the scattering medium, equation (4.11). Since the autocorrelation 
term at x = 0 is proportional to the electron density, equation (4.13) 
can be written 


< | G(t) | 2 > « n 


r oo 

/?(x)F(t, x)cix 

.00 


(4.14) 


where 


F(t,x) 


r 

b(t^)b [t^+T)b {t^-t)b{t^+x-t)dt^ 

.00 


(4.15) 


The value of F(t,x) is not affected by the characteristics of either the 
receiver or the scattering medium, but only by the particular pulse 
compression scheme used. 

For perfect correlation, /?(x) is equal to one and <|G(t)| > is propor- 
tional to the square of the autocorrelation function of the phase-coded pulse. 
This is approximately true also for cases where the length of the pulse is 
small compared with the correlation time of the scattering medium. The power 
measured at i = 0, i.e. the center peak of the autocorrelation function of 
the pulse, depends only on the total length of the pulse relative to the 
correlation time of the medium. 

One of the limitations in the use of this technique in incoherent 
scatter is that the pha;e coherence of the signal is partially lost during 
scattering. Gray and Farley [1973] showed that the loss of phase coherence 
is not significant as long as the total length of the pulse is less than the 
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time to the first zero crossing of the autocorrelation function of the 
scattering medium shown as approximately 300 ys in Figure 4.2. 

A detailed consideration of the measured signal in pulse compression 
experiments can be found in Gray and Farley [1973] . The range variation 
of > for various values of the ratio of total pulse length to 

correlation time of the medium is shown in Figure 4.4. The side lobes 
show the relative influence of echoes from other altitudes on the measurement 
of electron density. The contribution from other altitudes is not important 
if the electron density varies slowly on the scale of the altitude resolution. 

However, in the case of sporadic-E layers the difference in electron 
density from one altitude measured to the next can be quite significant. 

The effect is to decrease the apparent peak density of the layer and to 
increase the apparent density immediately above and below the layer. 

This smearing can be seen in an electron-density profile measured in 
this study that apparently includes a meteor trail, shown in Figure 4.5. 

As in the above derivation, the only scattering is from a very small 
altitude region. The positions of the side lobes with respect to the height 
gates compare with the range variation shown in Figure 4.4. In most sporadic-E 
layers, the peak electron density does not approach the value relative 
to the background ionization shown in Figure 4.5, and side lobes are not 
apparent in the data. 

The main effect of the contribution of side lobes to the scattered 
power is at altitudes immediately above and below the sporadic-E layer. 

At the peak of the layer, the relative contribution from other altitudes is 
decreased, with the result that the observed electron density relative to 
that of the E region is too low by approximately 10 % [ Gray and Farley 3 1973] . 

The reduction in signal from the sporadic-E layer depends on the thickness 
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4 The time (range) variation of <\G{t)\ > in dB below 
maximum value for the 13-baud Barker code for various 
T JT with T IT. = 1.0 [Gray and Farley, 1973], 
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Figure 4.5 The scattered-power profile of a meteor trail. Secondary 
maxima are artifacts of the pulse-compression technique. 
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of the layer and its position with respect to the range gates, and can not 
be determined accurately. For this reason, no correction was made for this 
effect in the analysis of the data. 

The electron-density profile is determined from the received power 
in the same way as in the single- and multiple-pulse experiments described 
earlier, with the height resolution being determined by the baud length 
of the phase-coded pulse instead of by the total pulse length. Theoretically, 
it would also be possible to measure other atmospheric properties by 
forming autocorrelation functions from sets of phase-coded pulses. The 
pulse -compression technique has been used to estimate the concentration of 
Fe + in a sporadic-# layer [ Behnke and Vickrey, 1975]. However, more 
development is needed to extend the technique to other properties. 

4.3 The Areciho Observatory 

The incoherent-scatter facility used in this study is located at the 
National Astronomy and Ionosphere Center, 15 km south of Arecibo, Puerto Rico. 
The geographic coordinates of the observatory are 18° 21' North, 66° 45! West. 
It is at a geomagnetic latitude of 30° North. The dip angle of the magnetic 
field is 50°; its magnetic declination is 8° 30' West. 

The antenna system is a fixed spherical reflector below a movable 
feed. The aperture has a diameter of 305 ra and a radius of curvature of 
265 m. A 430-NHz transmitter is used in the experiment, giving a beam 
width of approximately 300 m at an altitude of 100 km. The power output 
of the transmitter, when pulsed, is 2.5 MW. The maximum average transmitted 
power in a pulsed mode is 150 kW. The pulse width has a range of from 
2 us to 10 ms, with a pulse repetition period of from 1 ms to 1 s. For all 
experiments described here, the transmitter is controlled by an on-line digi- 
tal computer. 
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The spherical shape of the reflector causes the radiation to be 
focused into a line. For the 430-NHz transmitter a 30 m-long line feed 
is used. The feed is a waveguide suspended from a carriage house mounted 
on a curved track with a radius that coincides with the radius of the 
reflector. The position of the feed on the curved track determines the 
zenith angle. This track is suspended from a circular horizontal track, 
so the line feed can also be moved in azimuth. 

The steerability of the radar beam is of major importance to the 
experiments described here. The maximum slew rate of the beam is about 
23 deg min'* in azimuth, and 2.5 deg min * in elevation. Limits of the 
feed are 0 to 540 deg in azimuth and 0 to 20 deg in zenith angle. Because 
of the change in gain when the elevation is changed, the zenith angle 
of the beam was held constant in each experiment. The feed illuminates 
the complete surface of the reflector when it is vertical, but at a 
zenith angle of 20 deg the signal is reduced by about 4 dB. This reduction 
in signal becomes important only when the zenith angle approaches 20 deg. 

The gain does not change significantly with azimuth. 

When the scattered signal is received, it is first heterodyned to a 
center frequency of 30 MHz at the receiver located in the carriage house 
and then sent to the operations building, where further analog processing 
is accomplished. Depending on the requirements of the experiment this 
may include decoding and autocorrelation. 

The signal is then sampled by analog-to-digital converters, and 
immediately processed by the computer. Finally, the results, usually in the 
form of autocorrelation functions or scattered-power profiles, are stored 
on magnetic tape. 
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4.4 Experimental Techniques 

4.4.1 Electron density. Data were taken in the form of scattered- 
power profiles using a 13-bit Barker-coded pulse with a 4 -nsec baud length 
[Gray and Farley, 1973]. The power profiles are converted to electron- 
density profiles by calibration of the radar to the electron density at the 
peak of the F region obtained by the ionosonde . In this calibration the 
echo power is corrected for the effect of the value of T /T. being different 
in the E and F regions. For the set of experiments done in 1974, the 
value of T /T. in the F region is taken from Evans [1967]. In 1975, 

T and T. were measured in the multi -pulse experiment described above and 

& 7s 

assumed to have the same variation from one day to the next. In the E 
region adjacent to the sporadic-# layer T /T. is assumed to be equal to one. 

The data must also be corrected for range effects which include both 
the decrease in signal strength with the square of the range, and, for 
D- and E-region observations, a near-field correction to the gain [Rowe, 
1974]. 

Data were taken on one occasion with the beam held in a vertical 
position and on other occasions with the beam scanning in azimuth at a 
fixed zenith angle. The two modes were used to evaluate the potential of 
the radar for studies of horizontal structure. 

The 600-m resolution in range made possible by the Barker coding is 
important because of the anticipated thinness of the layers. Intense 
daytime sporadic-E layers observed in rocket measurements are typically 

1.5 to 3 km in total thickness while layers at night may be 1 km in total 
thickness [R-ri th, 1966; Srtr.th ana '4echci,y, 1972]. 

The height resolution of 600 m attainable using Barker-coded pulses 
is not sufficient to determine accurately the peak electron density of a 


sporadic-# layer. The result is that the peak electron density measured 
by incoherent scatter is consistently below the density expected, based on 
ionosonde measurement. Estimates of the effect of this volume averaging 
were obtained using electron -density profiles of sporadic-# layers from 
rocket-based measurements. This is illustrated in Figure 4.6 by the electron- 
density profile of a sporadic-# layer measured from a sounding rocket on 
10 August 1973. To simulate the effect of volume averaging the electron 
density has been averaged over altitude regions of 600 m, giving approxi- 
mately the same effect as the incoherent-scatter measurement. The 
averaging is done for two extreme cases, and a parabolic fit is made to 
the highest three points. It is seen that, even in the best case where 
the 600-m gate is exactly centered on the layer, the value of electron density 
given by incoherent scatter is less than the actual peak electron density. 

In the four layers examined, the ratio of the peak electron density obtained 
by the parabolic fit to the true peak electron density of the layers 
varied from 0.6 to 0.9. 

Since the horizontal variation of the electron density is of more 
importance to this study than the peak electron density of a layer, the 
problem of smearing the profile by the 600-m resolution in range is not 
critical. However, two effects of this smearing should be noted. First, 
if there is a change in the altitude of a layer, the apparent peak density 
will be affected. This can be seen by comparing the peak densitites given 
by the two averaging configurations in Figure 4.6. The farther the peak 
of a layer is from the center of a range gate, the smaller the apparent 
peak electron density will be. Also as J. D. Whitehead [personal communica- 
tion, 1974] has pointed out, regions of small horizontal extent would be 
unresolved if they fail because of size or position to fill the radar beam 
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ELECTRON DENSITY (l0 5 crrT 3 ) 


Figure 4.6 Effect of volume averaging of the electron density in 
600-m range- gates. Parabolas have been fit to average 
densities for two cases. 
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throughout the entire integration time. 

Another major difficulty in using the Arecibo Observatory to study 
sporadic-# by this method is the tum-on time delay of at least 15 min. 

When the detection of an intense sporadic-# layer by the ionosonde was 
used as a criterion for making an observation, the transmitter could not 
in many instances be brought up in time to observe the peak of the event. 

4.4.2 Ion drift velocity. Ion drift velocities in the altitude 
region 100 to 167 km were measured by incoherent scatter using a five -pulse 
technique described by Z amlutti and Farley [1975]. Autocorrelation functions 
are found from a series of five 8-us pulses at relative times of 0, 120, 

160, 360, and 440 sec. Various combinations of the pulses gave points of 
the autocorrelation f action at intervals of 40 us to a maximum time delay 
of 440 us, excluding the point at 400 us. Figure 4.2. Above 167 km, a 
single-pulse experiment with a pulse length of 296 us was used. 

In order to construct an ion drift velocity vector, measurements were 
made at three positions of the radar beam. The positions used were at 
180, 270, and 0 deg azimuth. The components of the ion velocity can be 
found from the line-of-sight velocities by 


V V 180' 2t, 270 
2 sin 6 


(4.16) 


V V 180 

V N = 2 sin 0 


(4.17) 


V 


z 


y^iso 

2 cos 0 


(4.18) 


where 0 is the zenith angle of the radar beam. For these observations, 
0=15 deg. 


The measured autocorrelation functions from each position of the beam 
arc fit by a least-squares method to a library of functions stored in the 
observatory computer system giving electron and ion temperature, collision 
frequency, and line-of-sight drift velocities. Error estimates are based 
on the accuracy of the least-squares fit. Estimated errors in T & and T \ 
are of the order of 5%. Drift-velocity errors are typically between 5 and 
10 m s" 1 . Errors were found to become very large for measurements taken 
in early morning or late evening, when the electron densities are small. 

No measurements are possible at night. 

Scattered-power profiles were measured during the five minutes 
required to swing the beam through the 90 deg between positions. The 
measurement used single 8-us pulses, giving a height resolution of 1.2 km. 
The profiles were stored on magnetic tape after integrating 300 samples, 
or every 3.5 sec, for a resolution in azimuth of about 1 deg. 

The measurement of winds at sporadic-# altitudes by the method of 
Zamlutti and Farley [1975] is limited by the minimum altitude at which 
measurements can be made. The correlation time in the E region is of the 
order of 500 us. Because of ground clutter, no signal can be used until 
the echo height of the last pulse of the set is above about 30 km, a time 
delay of 200 us. This means that to measure autocorrelation functions 
with time lags of up to 500 us the minimum range that can be used is 
105 km, representing a time delay of 700 us. Even with the radar beam at 
a zenith angle of 15 deg, the minimum altitude is higher than the altitude 
of occurrence of many sporadic-# laye7S. 

The correlation of wind data with sporadic-# layers was also hampered 
by the inability to make simultaneous measurements of the electron density 
and the ion velocities. The electron densities were measured with good 
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height and time resolution using Barker-coded pulses and short integration 
times. The ion velocity measurement, on the other hand, used a five-pulse 
autocorrelation scheme which required long (-15 min) integration times 
during which no electron-density profiles were recorded. Various combinations 
of these two experiments were used in an effort to correlate the data. A 
description of these experiments is given in the following sections. 

4.4.3 Data collection and presentation. Observations were made at 
the Arecibo Observatory during two periods: 11-28 January 1974, and 

10-27 July 1975. Measurements of electron density using the pulse-compression 
technique were made in 1974. The same experiment was used again in 1975 
together with a multiple-pulse experiment to measure ion drift velocities. 
Schedules of the two periods of observation, including the general features 
of the sporadic-# layers, arc presented in Table 4.1 for 1974 and Table 4.2 
for 1975. 

The data-taking program and observational procedures were tested on 
11 and 13 January 1974. Problems in recording the noise level in the 
scattered power were encountered and corrected. It was found that the 
system used for tracking astronomical objects with the radio telescope 
could be used to scan the radar in azimuth at a predetermined rate. 

Controlling the horizontal velocity of the point of intersection of the 
radar beam and a sporadic-# layer by adjusting the zenith angle of the 
radar was also considered. However, it was decided that required integration 
times were sufficiently short that the maximum angular velocity of the 
radar beam (-23 deg min”*) at the maximum usable zenith angle was feasible. 
Considering the reduction in signal of Ue radar at large zenith angles, 
and on the advice of the observatory staff, r ixed zenith angle of 15 deg 
was used. Sporadic-# layers are expected to move with a velocity oi the 


TABLE 4.1 


Observations during January, 1974, at Arecibo, P. R. 
Zenith 



Starting 

Ending 

Angle 

Speed 

Sweep 

Max 

Max 



Time 

Time 

of beam 

of beam 

Pattern 

fbEs 

ftEs 


Date 

(AST) 

(AST) 

(deg) 

(deg/min) 

(Azimuth) 

(MHz) 

(MHz) 

Remarks 

11 

1620 

1708 

15 

18 

60 o -120°-180 o -240° 

2.5 

3.2 

test of procedure 

11 

1823 

1858 

8.5 

23 

180° (Hold) -270° -90° 

2.5 

3.5 

test of procedure 

13 

1600 

1717 

15 

20 

60°-120° 

2.5 

3.2 

test of procedure 

13 

1719 

1743 

12 

23 

360° -540°- 180° 

0 

0 

test of procedure 

14 

1411 

1503 

12 

23 

90° -630° -90° -150° 

7 

9 

Scan full circle 

17 

2006 

2215 

12 

23 

60°-120° 

1.8 

2.8 

Night - 60° scan 

18 

1238 

1402 

12 

23 

60°-120° 

5.2 

8.0 

Scan 60° arc 

19 

1253 

1305 

12 

23 

150°-210° 

4.9 

6.0 

Scan 60° arc 

19 

1421 

1456 

12 

23 

1S0°-210° 

5.0 

6.0 

Layer died away and 









built back up 

24 

1803 

1901 

0 

— 

Stationary 

2.7 

4.7 

Two layers on ionogram 

28 

1237 

1510 

12 

23 

210°-270° 

6.0 

8.9 

Recorded layer build-up 


Table 4.2 


Observations during July, 1975, at Arecibo, P. R. 


DATE 

TIME (AST) 

MODE* 

MAX 

fbEa 

X uS 

REMARKS 

10 

11:34-13:33 

N 

5.2 

5.7 

Test of programs. 

16 

12:14-18:58 

W 

8.4 

11.0 

Kinds were measured that morning 
by R. M. Harper. 

17 

03:33-06:16 

P 

2.9 

6.9 

Scan between 270 and 310 deg. 
Intense layer formed at 04:30. 

17 

06:40-08:20 

W 

5.3 

7.1 

Large uncertainties, small 
electron density. 

19 

09:05-10:55 

w 

5.0 

6.3 

Problem with set-up. Measurement 
at 270 and 360 deg only. 

19 

11:45-14:22 

p 

> fotl 

12.8 

Scan between 265 and 275 deg. 
Observed irregularity with [e'J 
> 5.4 x 10 6 cm' 5 . 

20 

10:55-12:07 

w 

5.0 

7.0 

Large zonal component in ion 
drift . 

20 

12:17-13:37 

p 

3.9 

4.6 

Scan between 90 and 130 deg. 

20 

13:51-15:23 

N 

4.3 

5.3 

Ea layer not at ■ 0. 

21 

19:52-21:00 

and 

21:38-24:00 

P 

2.4 

6.0 

Scan between 270 end 310 deg. 
Transmitter failure at 21:00 
Ea layer formed at 21:30. 

23 

Ub: 10-08:05 

N 

2.8 

3.3 

No power profiles measured 
between first two positions. 

24 

N/A 

N 



Noise measurement. Antenna 
on load. 

26 

09:07-09:42 

P 

3.4 

4.2 

Scan between 270 and 310 deg. 

Ea on ionogram but not on profile 

26 

10:10-11:25 

N 

5.1 

8.8 

Smoothly-spiraling hodograph, 
no Ea layer. 

26 

11:33-14:50 

P 

5.4 

6.2 

Scan between 90 and 130 deg. 
Periodic enhancement of layer. 

26 

15:03-16:47 

M 

3.4 

4.0 

Large northward velocities below 
111 km. 

26 

19:44-21:37 

P 

2.5 

3.1 

Scan between 180 and 220 deg. 
Two weak layers merge at 94 km, 

27 

11:25-15:02 

P 

3.9 

4.3 

Scan between 180 and 220 deg. 
Ea layers at 95 and 110 km. 

27 

15:12-16:56 

W 

3.6 

4.3 

large errors in data. Possible 


problem with patch panel. 


* H ■ Five-pulse measurement of Inn drift with power profiles measured by a 
single-pulse technique during scans. 

P • Power profiles measured using Barker-coded pulses. 


ins 


order of SO m see'*. At a zenith angle of 15 deg, the horizontal velocity 
of the intersection of the radar line of sight with a layer at an altitude 

of 100 km is 190 ra sec -1 , when the rad.*r is scanning in azimuth at 23 deg 

. -1 
nun 

It is apparent that due to the small beamwidth (-300 m at sporadic-# 
altitudes) only a small fraction of the sporadic-# layer is actually 
sampled in the experiment. Figure 4.7 illustrates how a circular patch 
10 km in diameter and an edge of a lai'ge patch might appear in the data. 

In this example the beam is scanning in azimuth at the maximum rate, 

23 deg min -1 , and is held at a constant zenith angle of 15 deg. The 
sporadic-# layer is at an altitude of 100 km and is moving at a speed of 
50 m sec -1 . The motion of the layer is tangent to the center of the scan 
in the upper diagram and radial to the center of the scan in the lower 
diagram. Tne drawings on the left show the patteins traced by the 
intersection of the beam and the sporadic-# layer. The diagrams on the 
right are a convenient representation of the data. Each scan is plotted 
as a function of azimuth with the points where the beam encounters a 
target represented by shaded areas. The direction of the beam is opposite 
in adjacent scans. It can be seen that the size and position of the patch, 
as well as its velocity, have a large effect on the appearance of the 
patch on the data map. 

Another method of representing the deta is shown in Figure -!.8. Arrays 
of successive electron-density profiles, similar to the one shown on the 
left of the figure, are printed by a computer. Contours of equal electron 
density are then drawn, giving the impression of a vertical cross section of 
the ionization. This plot can be converted to distance coordinates in both 
vertical and horizontal axes by noting that, at a zenith angle of 15 deg. 
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Projection of 
Radar Beam 
on Ionosphere 


AZIMUTH (deg) 


wind velocity = 50 m/sec 


AZIMUTH (deg) 


wind velocity = 50 m/sec 


An illustration of the method used in the experiment and the effects of 
winds on the data [Miller and Smith, 1975]. 
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Figure 4.8 Illustration of a scan of the radar by drawing 

contours on a computer array of electron densities 
printed as functions of altitude and azimuth. 
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one degree in azimuth represents 468 m of horizontal distance along an arc 
at an altitude of 100 km. The horizontal distance per degree increases 
linearly with height. The azimuth scale is found to be slightly nonlinear 
because of changes in the angular velocity of the radar. 

4.5 Observations 

4.5.1 Introduction. A feature that becomes immediately obvious in 
examining the data taken by the incoherent-scatter radar is that sporadic-^ 
layers are as variable in structure as they are unpredictable in occurrence. 
Some have a simple vertical structure and are uniform in the horizontal 
plane, consistent with rocket observations. Others have complex vertical 
and horizontal structure as has been implied by the radio reflection 
observations. Because of the variety of features which were found in a 
total of 29 observing periods in 18 days it has been decided to first 
present the individual periods in chronological order, pointing out the 
features of each, and later to discuss the implications of these 
observations in the formation of sporadic-/? layers. 

Where appropriate, data are presented using the representations 
described in Section 4.43, as well as more conventional methods. Time is 
given as Atlantic Standard Time (AST) which is UT minus 4 hr. When data 
from ionograms are presented, the reflection characteristics of the part 
of the signal with ordinary polarization are examined by plotting values 
of foEs instead of ftEs. In all cases, foEs is determined by subtracting 
half the gyrofrequency from ftEs. 

4.5.2 14 January 1974 t 14:11-15:02 AST. Several methods of scanning 
the radar beam were tested. It was decided that the optimum mode was to 
scan repeatedly through arcs of 60 deg or less. In the course of the 
observation, the radar was scanned through one and a half revolutions in 
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azimuth, starting at 90 deg and continuing clockwise. The scan 
was then reversed. 

A sporadic-# layer appears to have drifted into the field of view 
at the end of the first scan (270 deg azimuth) and was seen again in 
the second scan. Contours of electron density observed during parts 
of the two scans are shown in Figure 4.9 as a function of azimuth and 
altitude. Figure 4.10 shows that while foEs remained nearly constant, 
there was an increase in fbEs between 14:30 and 14:40 AST as the 
sporadic-# layer passed overhead. 

4.5.3 17 January 1974 , 20:06-22:15 AST. This observation was made 

to gain experience with the system at night. A map of the data, similar 
to the example of Figure 4.7, is shown in Figure 4.11. The plot is of 
the peak electron density of the sporadic-# layer during successive scans 
as a function of azimuth. In this case the electron density is determined 
by a parabolic fit to the three points defining the layer. To aid in 
visualizing the variations, contours of electron density have been 
drawn, giving the impression of a map of the layer. This map may be 
considerably distorted as a result of the movement of the layer. 

The structure of this sporadic-# layer is relatively simple 
compared with subsequent nighttime observations. Correspondingly the 
variation in fbEs and ftEs, plotted in Figure 4.12, is relatively smooth, 
with ftEs reaching a maximum of 3.3 Mil z at 20:50 AST. By allowing 
roughly 25 min for the sporadic-# ionization to move the 30 km from the 
field of view of the ionosonde (0 deg zenith angle) to that of the radar 
(15 deg zenith angle), the peaks in foEs (Figure 4.12) are found to 
coincide with the observation of enhanced electron-density irregularities 
(Figure 4.11) . 
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Figure 4.9 Electron density (10 cm ) of a sporadic-^ layer observed 
during successive 90 deg scans on 14 January 1974. The 
motion of the radar beam is clockwise in the first scan and 
counterclockwise in the second. 
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Figure 4.10 fl'Es and foEs on 14 January 1974. 
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Figure 4.11 Representation of maximum density of the sporadic-/ 1 .' layer 
observed on 17 January 1974. The altitude of the layer is 
97 km. The direction of odd numbered scans is toward lower 
values of azimuth; *.ven numbered scans toward greater azimuth 









Figure 4.12 fbEs and foEs on 17 January 1974. 
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4.5.4 18 January 1974, 12:38-14:02 AST. In this and the remaining 

cases of January 1974 the observation was initiated on the basis of 
indication of the presence of sporadic E by the ionosonde. The layer 
appears to be made up of discrete patches of ionization separated by 
regions where no sporadic-E ionization was observed. The passage of one 
of these patches is shown in Figure 4.13. The discrete nature of the 
patches may not indicate a discontinuous metallic ion layer, but rather 
a variable layer with ionization densities in small regions becoming 
large enough to be observed above the ambient ionization. 

What appears to be three patches on the right side of the fifth scan 
in Figure 4.13 is the result of the electron density being averaged 
into two height gates at certain altitudes (compare with Figure 4.6). 

This effect is expected for both a tilted layer and a descending horizontal 
layer. In this case it is identified as a tilted layer, as the implied 
descent velocity would otherwise be excessive. 

The range of peak plasma frequencies derived from incoherent - 
scatter data is compared with fDEs and foEe in Figure 4.14. Two sporadic-E 
traces were observed on the ionograms. The difference in the frequencies 
from the two measurements is interpreted as being the result of volume - 
averaging by the radar. In this figure, the similarities between the 
variation of the lower foEa trace and the upper edge of the shaded area 
and between the variation of fuEa and the lower edge of the shaded area 
are apparent. This is consistent with the interpretation of foEa as 
being the largest value and fbEa the smallest value of the peak plasma 
frequency in the sporadic-E layer above the ionosonde. 


4.5.5 19 January 1974 , 12:53-14:56 AST. Between 14:21 and 14:25 AST 





a double sporadic-E layer with peaks at about 100 and 101.5 km was seen 






Figure 4.14 fbEs and foEs on 18 January 1974. Shaded area represents limits 
of plasma frequencies derived from incoherent-scatter data. 
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(Figure 4.15). A transmitter failure occurred during this observation. 

The upper layer was not as intense as the lower. No relative vertical 
motion was apparent between the two layers. Neither layer was observed 
after the scan shown at the bottom of Figure 4.15, although a much weaker 
layer formed at an altitude of 99 km after 14:40 AST. Two sporadic-!? 
traces, closely spaced in virtual height, were recorded on iono grams 
during this observation. 

4.5.6 24 January 1974 s 13:03-19: 01 AST. On this occasion the radar 

beam was held in a vertical position and a sporadic-# layer observed. 

Figure 4.16 shows the excess electron column density of the layer. The 
length of the column (1.8 km) is defined by the three height gates centered 
on the layer. The electron densities above and below the layer show a 
smooth decrease in ionization throughout the experiment. They were 
averaged to determine the background ionization. The background is a 

small part of the total ionization in this region and varied from 

9-2 9-2 

0.94 x io cm at the beginning of the observation to 0.38 x 10 cm 

at the end. During the one-hour period of observation the solar zenith 

angle changed from 95 to 105 deg. 

The peak of the layer is observed initially to be :»t about 100 km 
and descends during the period by about 200 m, based on the parabolic fit. 

The signals have been averaged over 20,000 samples, a period of approximately 
64 sec. The error bars represent statistical errors based on the number 
of samples averaged and the signal -to-noise ratio, as described by 
equation (4.13). 

The variation in excess electron column density is interpreted tc be 
a spatial variation and not a variation in time, since no major 
redistribution in altitude of the ionization is apparent in the vertical 
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Figure 4.16 Electron column densities through a sporadic-!!, layer 
at an altitude of 100 km observed on 24 January 1974 
The radar beam was held in a fixed vertical position 

[Miller and Smith, 1975]. 
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electron-density profiles. If a wind speed of 50 m s 1 is assumed, 
this layer has a horizontal extent of 165 km. In this observation, 
very little small-scale structure is seen. 

p igure 4.17 shows the values of fbEs and foEs from ionograms taken 
daring this experiment. It may be noted that fbEs and foEs are 
approximately equal on this occasion. The temporal variation of the 
column electron density shown in Figure 4.16 is quite similar to that of 
the blanketing and critical frequencies shown in Figure 4.17. 

In examining the data from 24 January 1974, an interesting feature 
was noted which illustrates the variability that can be present in the 
E region at sunset without appreciably affecting a sporadic-F layer. 

A sequence of electron-density profiles of the E region is shown in 
Figure 4.18. These profiles show what appears tc be a layer of ionization 
ascending rapidly from 125 to 155 km. The time between profiles is 
approximately 3 min. During the observation, the solar zenith angle 
increased from 95 deg to 105 deg. There is also a rapid decrease in 
the ionization density of the E region which is better illustrated in 
Figure 4.19. The apparent upward movement of the layer below 150 km appears 
to be the result of both recombination and vertical transport, and 
does not necessarily represent only the bulk motion of the ionization. 

The display of incoherent*scatter electron-density data by staggering 
vertical profiles, as in Figure 4.18, is very popular in the literature. 

The comparison of these two figures emphasizes that care must be taken 

in the interpretation of features related tc the electron density in such plots. 

4.5.7 28 January 1974 , 12:52-14:88 AST. On this occasion a complete 
sporadic-F event was observed, including both build-up and decay. The 
layer was observed at an altitude of 105 km. A map of peak electron densities 
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Hgure 4.17 fbEs and foFs on 24 .January 1974 [Miller and Smith, 1975]. 
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based on a parabolic fit is shown in Figure 4.20. The contours are 

slightly oblique with respect to the scans, indicating that the movement 

of the sporadic-E layer is in a direction essentially tangent to the 

scan. Figure 4.20 shows two regions containing patchy structure with 

electron densities higher than the surrounding layer. The first, between 

13:45 and 13:50 AST includes the highest electron density recorded: 

5 -3 

3.1 x 10 cm ; the second, between 14:10 and 14:25 AST has a maximum 
electron density of 2.3 x 10 5 cm 3 . 

Vertical-electron-density profiles show that the column electron 
density of the sporadic-E layer increases when an increased peak electron 
density is observed, indicating that the irregularities are not caused 
by redistribution of the ionization in altitude. This is demonstrated 
by Figures 4.21 and 4.22, in which contours of electron density are 
plotted as a function of azimuth and altitude. The two regions where 
irregularities are observed are illustrated in these figures, which corre- 
spond, respectively, to scans in Figure 4.20 numbered 17-21 and 29-55. 

The thickening of the sporadic-E layers at the regions of enhanced 
electron density is apparent. 

The values of fbEs and foEs during the time of observation are 
shown in Figure 4.23. For part of the time two sporadic-E echoes are 
present; both are shown in Figure 4.23. The higher trace was seen to 
descend quite rapidly, possibly the result of off-vertical reflections. 
However, at no time was the virtual height of the higher trace observed 
to increase. The electron-density profile shows only 3ne layer to be 
present. Rocket observations by Smith [1970] similarly showed that not 
all sporadic-E echoes are associated with identifiable features of the 
electron-density profile. 
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Figure 4.20 Representation of maximum density of the sporadic-/? layer observed on 
28 January 1974 between 12:52 and 14:58 AST. The altitude of the 
layer is 105 km. The direction of odd-numbered scans is from left tc 
right, even-numbered scans from right to left [ Miller and Smith, 1975 
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Figure 4.21 Electron density (10 cm 15 ) above Arecibo, 28 January 1974, 
between 13:40 and 13:54 AST. Each scan represents approxi- 
mately 3 min. 
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Figure 4.23 fbEs and foEs on 28 January 1974. Shaded area represents limits 
of plasma frequencies derived from the incoherent -scatter data 
[Miller and Smith, 1975]. 
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The shaded area in Figure 4.23 is bounded by the plasma frequencies 
derived from the maximum and minimum values of the peak electron density 
of the layer during each scan in the incoherent-scatter data. The 
blanketing and maximum frequencies from the ionosonde data are both 
seen to be generally higher than the plasma frequencies derived from the 
incoherent scatter data, an effect of volume averaging. 

4.5.8 10 July 1975 , 11:31-13:32 AST. The first observation of the 

second series was devoted to a test of the data-taking program for the 
multiple-pulse ion drift experiment. Ion drifts were measured with the 
radar at azimuths of 180, 270, 360, and again at 180 deg. These positions 
will be identified as 1 , 2, 3, and 4, respectively. 

As described in Section 4.4.2 the three components of ion velocity 
can be obtained by combining measurements from positions 1, 2 and 3 or 
from positions 2, 3 and 4. Hodographs of the horizontal ion drift 
velocities and the corresponding profile of vertical ion drift deri.ed 
separately from measurements at the two sets of positions are shown in 
Figure 4.24. The results from both sets of data would be identical 
if the velocity components were constant and errors in measurement of 
velocity were negligible. It is seen that there was in fact a significant 
change in the ion drift velocity in the 1.5 hr required to make the 
measurements. As a result of this test observation an interpolation 
procedure was developed which allows for slow changes in the wind 
structure. This will be illostrated in subsequent examples. 

While moving between the four positions electron-density profiles 
were obtained with an altitude resolution of 1.2 km. These show a 
sporadic-# layer at a constant altitude of 109.5 km. As indicated in 
Figure 4.24 this layer occurs near a convergence point in the vertical 
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Figure 4.24 Ion drift velocities, 10 July 1975. The solid line shows 
velocities from incoherent-scatter measurements between 
11:31 and 12:39 AST, the dashed line from measurements 
between 11:56 and 13:07 AST. 
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ion drift velocity. During this observation, ionosonde measurements 
showed fi-Es to be approximately equal to foEa , and vary between 3.6 
and 5.2 KHz. 

4.5.9 16 July 1976, 12:14-18:58 AST. On this day, an extended 

observation of ion drift velocities was made. Values of fbEs and 
foE8 from the ionograms are plotted in Figure 4.25. An intense sporadic-# 
event commenced at about 11:45 AST. No incoherent-scatter data were 
taken until 12:14 AST. At 12:59 AST, the value of ftEs reached a peak 
of 11.0 NHz, but by the time a third wind vector was measured and ion 
drift velocities determined, it was 13:22 AST, and ftEs ha' 1 dropped 
to 7.0 MHz. The sporadic-!? layer persisted throughout the afternoon 
with values of ftEs ranging between 4.5 and 7.6 Miz. The experiment 
was continued until 18:58 AST. 

The ion velocities measured at each position of the radar were 
interpol '■ed to a common time in order to minimize the effect of the 
changing velocities. This is justified by the continuity apparent 
between measurements. Figure 4.26 shows the line-of-sight ion velocities 
at altitudes of 100.4, 104.0 and 107.5 km, measured at each position of 
the radar and plotted at .he corresponding time. The error bars are 
based on the fit of the measured autocorrelation function to a theoretical 
autocorrelation function. The relative magnitude of the error bars to the 
Cartesian velocity components is illustrated in the wind profile shown 
in Figure 4.27, which is based on values interpolated to 13:12 AST. 

Figure 4.28 shows the components of the ion drift velocity as a 
function of time. The changing position of the sporadic-# layer is 
indicated. In spite of the presence of a sporadic-# trace on the 
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Figure 4.25 fbFs and foEs on lb July 1975. 


it* 

Ln 










ALTITUDE (Km) 


) \ =30^ J 

J \ V -2S / 


/ \ 


\ •'■' 




• . ^--?o- 


ir»-'o 

o O 

°? 

s 

d*-> 


/ l \ -as 

/ - <C^ -'3\ -20 s * — ,g — 




3 ^ 

(G,) £ (o) (o) (Q)> 


^-$5- --t 

^ 

/ 

' M 




IlOh ~3Q 



2>3U SO\ 




<C _C!!> 

■-23 


O o 


~2s~ysb~ 

( 


) ) 

< C 


(O) 10! (O) 


I00 1 — 

1300 


29 


'-SO 


(O) . 23— O 

3(H 


TIME (AST) 


Figure 4.28 Ion drift velocities, 16 July 1975. The altitude of 
a sporadic-E layer is indicated by the circled points 
Parentheses indicate a weak layer. 
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ionograms, no layer was detected by the incoherent-scatter experiment 
between about 14:40 and 16:00 AST. 

It is also apparent in this figure that the sporadic-# layer is 
not at an altitude at which the vertical ion drift converges. The general 
downward velo-ity is probably the result of an inaccurate estimate of an 
internal frequency bias in the system [Harper et alo, 1975; 1976], Correction 
of this bias would not produce vertical convergence at the altitude of 
tiic layer. However it will be shown later that the altitude of the layer 
is consistent with wind shear theory. 

4.5.10 1? July 1975 t 02:33-08:20 AST. Electron density was 

measured using phase-coded pulses until 06:16 AST, followed by an 
ion drift velocity measurement. 

Figure 4.29 shows a sequence of scans where contours of electron 

density are drawn on a plot of altitude and azimuth. The contours 

4 5-3 

represent values of electron density of 1 x 10 and 1 x 10 cm . Each 
scan covers a range of azimuth of approximately 40 deg and is centered 
on 290 deg. The vertical and horizontal scales in each frame do not 
give a true spatial representation of the region. The 40 deg range in 
azimuth represents a distance along the arc of 18.7 km at an altitude 
of 100 km, and 24.3 km at an altitude of 130 km. Beginning with the first 
frame in the upper left-hand corner of Figure 4.29, the actual motion of 
the radar beam is from right to left in the odd-numbered scans, and from 
left to right in the even-numbered scans. 

The most striking feature of Figure 4,29 is the complex structure 
exhibited by the layers of ionization. At the beginning of the sequence, 
layers are observed at altitudes of 123 and 108 km. The upper layer 

remains at a constant altitude for approximately 20 min, and then begins 

l 


Figure 4.29 Electron density observed on 17 July 1975. Contours are at 104 and 10° cm 

Each scan of the radar beam was between approximately 270 and 310 deg azimuth 
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to descend dissipating at 04:35 AST at an altitude of about 116 km. 

If the descent is extrapolated, it intercepts the lower layer at 
04:40 AST corresponding with an intensification of the lower layer, 
which, however, remains at constant altitude. The peak electron 
density in this lower layer is dramatically enhanced, and the profile 
shows evidence that the layer becomes unstable. Another layer is seen to 
descend into the layer at 108 km between 04:12 and 04:26 AST, with less 
effect. The diffuse patch of ionization located just above the lower 
layer at 03:52 AST may also be noted. Each of these intensifications 
can be identified with increases in fbES and foEa as shown in Figure 4.30. 

An attempt was made to measure the ion velocities between 06:40 and 
08:12 AST to aid in the interpretation of the electron-density measurements. 
However, it was found that the low electron density at that time resulted 
in large uncertainties and improbably large velocities and no useful data 
were obtained. 

4.5.11 19 July 1975 , 11:05-14:22 AST . Perhaps the most dramatic 
evidence of patchiness in a layer associated with a high value of foEa 
is found in data taken this day using Barker-coded pulses . A particularly 
intense sporadic-^ layer was present for the entire period. The range of 
plasma frequencies derived from the maximum and minimum peak electron 
densities during each scan is compared with fbEa and foEa in Figure 4.31. 

In this experiment, the beam was scanned through approximately 10 deg in 
azimuth centered at 270 deg. 

The irregularity which is represented by the high plasma frequencies 
near 12:55 AST is shown in more detail in Figure 4.32. The electron- 
density configuration is shown for twelve successive scans. The azimuth 
scale has been adjusted so that horizontal and vertical distances 
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Figure 4.30 fbEs and foEs on 17 July 1975. 
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Figure 4.31 


fbEa and foEa on 19 July 1975. Shaded area represents limits of plasma 
frequencies derived from incoherent-scatter data. 
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represented in the plots are approximately equal. At the time of this 

observation, the electron density of the E region adjacent to the sporadic-# 

4 5 -3 . 

layer varied between 9 x 10 and 1.3 x 10 cm . In Figure 4.32 the 
dashed line represents an electron density contour of 5 x 10^ cm 
the solid lines are contours at intervals of 1 * 10^ cm 

Figure 4.32 shows a very localized irregularity either forming and 
dissipating rapidly, or drifting through the radar beam. The map of 
the irregularity as a function of azimuth and scan number is presented in 
Figure 4.33. Here the volume- averaged electron density recorded by the 
incoherent scatter is used, without attempting to estimate the maximum 
electron density. The short scan used in this experiment (-10 deg) shows 
the motion of the irregularity along the arc intersected by the radar 
beam. However, lack of detailed knowledge of ion drift velocities makes 
it uncertain whether the apparent motion is the result of a changing 
velocity of the irregularity in a direction tangent to the scan, or 
whether the velocity is radial and the irregularity has an elongated 
shape. 

The irregularity is not resolved by the radar and is no larger than 
either the width of the radar beam (300 m) or the range resolution (600 m) . 
The effect is that the peak electron density of the irregularity cannot 
be determined. The irregularity is recorded in two gates, part in one 
at 113.05 km and part in one at 113.63 km. The volume -averaged values at 
these gates at 12:57:42 AST are 5.41 x 10 6 and 5.18 x io 6 cm~\ respectively. 
The peak electron density must be greater than, and may be much greater 
than, the larger of these two values. 

4.5.12 20 July 1975 , 10:55-15:23 AST. Ion drift velocity measurements 
using the multiple-pulse experiment were made from 10:55 to 12:07 and from 



Figure 4.33 





13:51 to 15:23 AST. Hodographs of the horizontal velocities are shown 
in Figure 4.34, together with a profile of vertical velocities. The 
zonal (east-west) component is seen to be predominant, especially during 
the morning observation. 

Sporadic-# layers were observed at altitudes of about 102 km in 
the morning and 103 km in the afternoon. In both cases, the layeis were 
located on a strong east-west shear in the ion drift velocity. However 
as in the earlier case described in Section 4.5.9, their positions do not 
correspond to 1 ocations of vertical convergence . 

Electron-density profiles were observed by Barker-coded pulses 
between 12:17 and 13:37 AST. The sporadic-# layer was found to be centered 
at an altitude of 102.5 km. The structure observed in this layer is 

4 

illustrated in Figure 4.35. The interval between contours is 2.5 * 10 cm 
This observation is an example of rather minor irregularities in a 
spoiadic-# layer. 

4.5.13 21 July 1975 t 19:52-24:00 AST. The observation was begun 
at 19:52 with no sporadic E indicated in the data from the ionosonde, 
shown in Figure 4.36. The transmitter was inoperative between 21:00 
and 21:38 AST. During this interval ftEs built up to 6 MHz. By the time 
the observations were begun again, ftEs had dropped to 4.2 Mlz, and 
continued to decrease until the sporadic-# trace disappeared at 22:05 AST. 
The incoherent -scatter radar detected a sporadic-# layer with a thickness 
of about 4 km centered at an altitude of 112 km. The peak electron 
density of this layer decreased steadily until 22:06 AST after which 
it was not detected. 

Later, starting at 22:27 AST, a layer was observed which developed 
unusual structure. Consecutive scans are shown in Figure 4.37. The 
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Figure 4.37 Electron density (10^ cm'^) from a sequence of 
scans of the radar beam on 21 July 1975. 
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range of scan is 40 deg, centered on 290 deg. The layer of ionization 
had formed at an altitude of 127 km. The ionization at one end of the 
region being observed suddenly began to descend. By 23:12 AST, this 
ionization had re-formed into a layer at 106 km. At the same time, the 
part of the layer at the other end ascended to 135 km. 

4.5.14 23 July 1975 a 03:11-08:05 AST. Electron-density profiles 
were measured from 03:11 to 06:00 AST with no sporadic-E layer being 
observed. The ion drift velocities, measured between 06:10 and 08:05 
AST, have large uncertainties. However, the magnitudes of the velocities, 
shown in Figure 4.38, are not unreasonable. A sporadic-E layer was 
observed during the multiple-pulse experiment at an altitude of 101.5 km, 
which was also near a convergence point of the vertical motion. 

4.5.15 26 July 1975 , 09:07-16:47 AST. It had been noticed that iono- 
grams often showed sporadic-E traces in the midmorning as was the case on 
this day. The data on foEs and fbEs are shown in Figure 4.39. An incoherent - 
scatter-radar observation was made from 09:07 to 09:42 AST. Although ft-Es 
reached 4.2 MHz, no sporadic-E layer was detected by the incoherent -scatter 
radar. However, the vertical-electron-density profiles appeared to have a 
very irregular structure. A profile measured at 09:12 AST, a maximum of 
ftEs, is shown in Figure 4.40. Any irregularities observed in the vertical 
profile on this occasion were of limited horizontal extent, and do not 
represent layering. This is the only occasion when electron-density 
observations were attempted in midmorning and structure of this type was 

not seen at other times of day. 

Ion drift velocities were measured during two periods, 10:10 to 11:25 AST 
and 15:03 to 16:47 AST. Horiztonal and vertical velocities for both 
periods are shown in Figure 4.41. Data from the morning observation 
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Figure 4.38 Ion drift velocity, 2? July 1975 
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Figure 4.40 Electron-density profile measured by incoherent 
scatter radar at 09:12 AST on 26 July 1975. A 
three point weighted average was used to smooth 
the profile. 









shows a smoothly spiraling hodograph and no nulls in the vertical-ion- 
velocity profile. No sporadic-# layer was present in scattered-power 
profiles, although the ionosonde showed fbEa = 3.5 KHz and ftEs = 4.0 NHz. 

The afternoon ion drift velocity measurement shows strong velocities 
to the north at the lower altitudes, ceasing abruptly at about 111 km. 

A sporadic-# layer at 101.5 km is on shear in the east-west ion velocity. 
However, both the sense of the shear and the vertical velocity information 
indicate that the ionization should diverge at that altitude. 

Between the two velocity measurements electron density was observed, 
from 11:33 to 14:50 AST, using Barker-coded pulses. A weak sporadic-# 
layer was found at an altitude of 97 km. Plasma frequencies of the maximum 
and minimum peak electron densities are compared with fbEs and foEe in 
Figure 4.42. Enhancements of the layer with a periodicity of approximately 
15 or 20 min are apparent in the sporadic-# layer. The enhancement of 
the electron density of the layer is accompanied in most cases by development 
of wave-like structure. Six scans showing the enhancement at 13:39 AST 
are presented in Figure 4.43. 

4.5.16 26 July 1975, 19 :44-21: 37 AST. Electron-density profiles were 
again measured between 19:44 and 21:37 AST. A very weak sporadic-# layer 
was observed at an altitude of 93 km. A second layer begins to form at 97 km 
at 20:40 AST. The two layers, shown in Figure 4.44, merge at 21:09 AST and 
form a 4-km thick layer centered at about 96 km. The layer then descended 
slowly, until at the end of the observation it was centered at 94 km. 

4.5.17 27 July 1975 , 11:25-15:02 AST. Electron-density profiles 
recorded between 11:25 and 15:02 AST showed weak layers near 95 and 120 
km. The profile at 12:55 AST is shown in Figure 4.45. No significant 
small-scale horizontal structure was observed. Although a weak F-region 


13 = 


TIME (AST) 


j-igure 4.42 fbE 's and foFr; on 2h July 1975. Shaded area represents 
limits o f plasma frequencies derived from incoherent- 
scatter data. 
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Figure 4.44 Electron density from a sequence of scans 

of the radar beam on 26 July 1975. Contours 
are at intervals of 5 x 10 -5 cm - -*. 
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figure 4.45 Electron-density profile measured by incoherent- 
scatter radar at 12:55 AST on 27 July 1975. A 
three-point weighted average was used to smooth 
the profile. 
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trace on ionograms made determination of fbEs impossible, ftEs remained 
nearly constant at 3.9 Mlz. 

4.6 Comparison of lonoscnde and Incoherent-Ecatter-Ftadar Observations 

4.6.1 Effect of volume averaging. 'Hie horizontal variation seen 
in the sporadic-# layers discussed in Section 4.S shows a correlation 
between the presence of small regions of enhanced electron density and 
the partial transparency of sporadic-# layers at frequencies greater than 
fbEs when allowance is made for a time delay resulting from a separation 
of the regions of the layer being sampled by the incoherent-scatter radar 
and ionosonde (see Figures 4.12, 4.14, 4.23, 4.30, 4.31, and 4.36). In 
Figure 4.31, the F region was blanketed by the sporadic-# layer so that 
fbEs could not be determined during the time when the intense irregularity 
was observed. However, a large increase in foEs occurred at that time 
when the irregular structure was observed. 

In the observations presented here, the peak electron density given 
by incoherent scatter is less than the peak electron density implied by 
the ionosonde observations. In the case of the layer observed on the 
evening of 24 January 1974, Section 4.5.6, the magnitude of the electron 
density at the peak of the layer given by the incoherent-scatter experi- 
ment is 28 percent of the density derived from the blanketing frequency. 

This can be explained by postulating that the thickness of the layer 
is less than the altitude resolution of the radar (b00 m) . A triangular 
sporadic-# profile that is approximately 300 m in total thickness would 
reconcile the data from the ionosonde with the data from incoherent scatter. 
It is apparent that with such thin layers the present meth of investi- 
gation is unable to determine accurately the exact peak electron density 
of the layer. 
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The layers with large values of peak electron density were found 
to be substantially thicker than the evening layer just described, 
usually contributing to the elect ... sity in three successive 
height gates. Peak election densities appear to be about 70 percent of 
densities indicated by the ionoscnde aa.a. In Section 4.4.1, using 
electron-density profiles for rockets, it was shown that the effect of 
volume averaging in the incoherent-scatter data is to give peak electron 
densities which are between 50 and 90 percent of the actual peak electron 
density, the variation depending on the actual location of the range gate rela 
tive to the layer. This effect accounts for the range of plasma frequencies 
(i.e., the shaded area in the figures) being below the region bounded by 
fbEs and foEs. 

4.6.2 Effect of limited sampling. Although small regions with rela- 
tively high electron densities were observed, for example, on 28 

5 -3 

January 1974 (Figure 4.23), the highest density was 3.1 * 10 cm , 

5 -3 

whereas an electron density of 7.1 ' 10 cm is required to explain foEs 
on the basis of total reflection. This apparent low electron density is, 
in part, the result of volume averaging. Interpretation of the data must 
also take account of the percentage of the layer actually being sampled 
by the incoherent-scatter radar. An estimate of the motion of the layer 
observed on 28 January 1974 can be made from the rate of drift of the 
patches across the data map shewn in Figure 4.20 and from the time lag 
apparent in Figure 4.23 between features in the fields of view of the iono- 
sonde and the radar. The sporadic-E layer appears to be moving at approxi- 
mately 85 m s * on an azimuth of 190 deg. The scan pattern of the beam. 

Figure 4.46, is confined to a strip 17.5 km in width. The beam, 300 m in diam 
eter at 100 km altitude, samples 4.5% of the total area of the strip. 


Figure 4.46 Beam pattern in the reference frame of the sporadic-.? 

layer observed on 28 January 1974. Motion of the 
layer is 85 m s -1 at an azimuth of 190 deg. The radar 
is scanning between 210 and 270 deg. 
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A statistical analysis was made of the data taken on 28 January 1974 
between 13:45 and 13:52 AST, when electron-density irregularities were ob- 
served. A histogram of the number of observations of peak electron densities 

4 -3 

in intervals of 3 x 10 cm is shown in Figure 4.47. Each observation is 
the result of integrating 3000 samples and represents an average electron 
density in approximately 4 s of time or 2 km of horizontal distance. The 
frequency of occurrence of electron densities greater than 1.4 x 10^ cm ^ 
decreases exponentially with increasing electron density, indicating the 
possibility that higher electron densities were present within the layer, 
but not observed. 

As was mentioned, only about 4.5 ? o of the region is actually sampled by 
the incoherent-scatter measurement. If it were possible to measure elec- 
tron densities in the entire region, the number of observations in Figure 
4.47 would presumably be increased by a factor of (0.045) The extended 
exponential decrease of observations with electron density shown in Figure 
4.47 would then predict one observation at the electron density where the 
present curve is equal to 0.045. 

If it is accepted that the values of electron density shown in Figure 

4.20 are underestimated by a factor of 0.7 because of volume averaging, 

5 -3 

then an electron density of 7 * 10 cm , corresponding to the maximum 

value of foEs in the ionogram (7.5 MU z) , would appear as an electron 

5 -3 

density of 4.9 * 10 cm in the incoherent-scatter data. 

Following the above line of reasoning, a point has been included in 

5 -3 

Figure 4.47 at an electron density of 4.9 * 10 cm and a number of observa- 
tions equal to 0.045, as predicted by the ionosonde observation. There is 
good agreement between this point and the extension of a least-squares fit 
to the incoherent-scatter data. Thus, it is probable that an 



I 


i 




167 



Figure 4.47 Frequency of occurrence of peak electron densities 

between 13:45 and 13:50 AST on January 1974 (circles'), 
and prediction based on maximum ftF.s (square). 
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electron-density irregularity existed in the sporadic-# layer with a peak 
density of 7 x 10^ cm ^ , corresponding to the maximum value of foEs . 

The observation of 19 July 1975 shows the converse effect in which, be- 
tween 12:50 and 13:00 AST, the peak electron densities observed by incoherent 
scatter were much greater than expected on the basis of foEs. Two factors 
may have prevented the ionosonde from observing the echoes at 21 MHz, in- 
dicated in Figure 4.31. First, the section of the sporadic-# layer observed 
by incoherent scatter was not in the field of view of the ionosonde. 

Second, highly-ionized irregularities comprise only a small fraction of the 
sporadic-# layer. It was stated in Chapter 3 that the ionosonde would re- 
cord reflected radiation from a sporadic-# layer only if the reflection 
coefficient of the layer were greater than approximately 1 * 10 This 
means that at least 0.1% of the field of view of the ionosonde would have 
to contain a high enough electron density to reflect radiation of a given 
frequency before an echo at that frequency would be recorded on the iono- 
gram. The fact that only one irregularity was observed, although there 
were probably others in the vicinity, indicates that they are scattered 
quite thinly throughout the region. This is evidently the reason that these 
regions of high electron density have not been detected in data from sounding 
rockets. 

4 . 7 Summary 

These series of experiments have demonstrated the value and limitations 
of the incoherent-scatter technique in the study of the horizontal structure 
of sporadic-#. The observations show that considerable horizontal structure 
does exist in some intense midlatitude sporadic-# layers. The size of this 
structure ranges from hundreds of kilometers down to the spatial resolution 
of the radar beam, about 300 m at a range of 100 km, with implications of 
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even smaller structure. Small regions were observed with densities much 

greater than the maximum electron densities in the surrounding sporadic-# 

layer. The occurrence of the irregularities was correlated in time with 

the detection by the ionosonde of partially transparent echoes from the 

layers. The magnitude of the irregularities was also found in each case 

to be consistent with the values of fbEs and foEs , when the effects o* 

volume averaging and limited sampling by the radar are considered. 

0 

The observations have interesting implications for the explanation 
of the formation of sporadic-# layers. If, as appears likely, sporadic-# 
layers are formed by the vertical drift of metallic ions under the influence 
of the neutral wind system, then it is necessary to postulate either 
horizontal convergence of the ionization or an irregular source of the 
metallic ions. 
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5. DISCUSSION 

5.1 Introduction 

It was shown in Chapter 4 that, based on data from incoherent-scatter 
radar, there is small-scale structure of sufficient magnitude to account for 
the partial transparency of sporadic-# echoes seen in ionograms. The causes 
of sporadic-# irregularities will be discussed in this chapter. 

The irregular structure observed in sporadic-# layers requires either 
that the ionization present is subject to an irregular forcing mechanism, 
or that there is an irregular source of ionization. The effects of motions 
of the neutral atmosphere on sporadic-# ionization are investigated. In 
particular, it is shown that many of the irregular features observed in 
sporadic-# layers may be caused by internal gravity waves. The possible 
proximity of a sporadic-# layer to a critical level of a dominant gravity- 
wave mode is discussed as a sou c ' of intense ionization irregularities. 
Evidence is shown of the peak ioni:;ation density in a layer being limited 
by a fluid instability. G.her effects of motions of the neutral atmosphere 
which could give rise to small-scale irregularities in sporadic-# layers 
are discussed. 

A possible irregular source of ionization from ablating meteors is also 
considered. It is shown that under specific and rather rare conditions an 
ionized meteor trail can be converged by •‘.he action of the neutral wind 
into an irregularity of relatively small dimensions. 

5.2 Ion Drift Velocity and the Formation of Spcradic-E Layers 

Seme of the features of sporadic-# layers predicted by wind-shear 
theory were discusived in Chapter 2. Although, i.i general, observations of 
sporadic-# layers support the conclusions of wii.d-shear theory, not all of 
the problems have been resolved [ Layzer, 1964a, b, 1972; Whitehead, 1964, 
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1971b; Axfordy 1964]. It is, therefore, of interest to examine the results 
of this study for features either characteristic of or in contrast to the 
predictions of wind-shear theory. 

It was found in describing sporadic-# profiles by mathematical functions, 
as described in Chapter 3, that most layers had a Gaussian shape. As 
discussed in Section 2.2.4, this is consistent with wind-shear theory for 
a linear shear if the layers are diffusion limited, as would be the case 
for metallic ions. 

The concept of diffusion-limited layers may provide an explanation of 
one of the apparent discrepancies between the results of this study and the 
predictions of wind-shear theory. In Section 4.S.9 an experiment was 
described in which the ion velocities in the region of a sporadic-# layer 
were measured for an extended period. Figure 4.28 showed the ion velocities 
and also the location of a sporadic-# layer. It will be noted that no 
convergence in the vertical ion drift was measured in the vicinity of the 
sporadic-# layer. In fact the ionization seems to be moving uniformly 
downward at a speed of approximately 5 m s *. The apparent downward motion 
is likely to be due to an uncertainty in the internal bias of the radar 
[Harper et al. 3 1976]. However, correcting for this velocity bias does 
not create a convergence of ionization. 

The effect on the velocity components of an incorrect estimate of the 
internal bias of the radar can be easily found. The ion velocity vectors 
whose components are plotted in Figure 4.28 are calculated from measured 
velocities parallel to the lines of sight of the three positions of the 
radar by the relations given by equations (4.16) through (4.18). Since the 
velocity bias affects the measured velocities equally in each position, its 
effect is doubled in the vertical velocity, but cancels in the horizontal 
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velocities. The horizontal velocities are therefore not affected by the 
internal bias. 

Given the horizontal ion velocities (y and v ), the vertical component 

x y 

of ion drift (y ) can be calculated using the equations of Chapter 2. If 
z 

u is equal to zero in equation (2.17), ion drift velocity components are 
z 

given by 


v = — l —=- r (r 2 +p. 2 )u ♦ (r r +p.r )u ] 

x 2. x x x y v z J y 1 


(i + p; ) 


(5.1) 


y - ■ * K r «V>iU“* + VyW^x 

,J ) 


(5.2) 


y = =- [(r r +p.r )u + (r r -p.r )u ] 

2 £ 1+p 2^ x 2 v y J x y z z X y 1 

1 * 


(5.3) 


Equations (5.5) and (5.6) car be solved in terms of u and u and 

x y 

substituted into equation (5.7) to give 


y = — (r y -r y ) 

- p i v y X X y 


(5.4) 


The vertical ion drift velocities deduced in this fashion from the data of 

Figure 4.24 below 115 km are shown in Figure 5.1. 

The position of the layer observed in the electron-density profiles is in 

good agreement with the calculated values of u . The layer is situated near 

z 

a convergence point, and then dissipates at 14:40 AST when convergence is no 
longer indicated. However, there is a weak layer at 104 km between 15:30 
and 17:40 AST. Although the vertical ion drift does not change direction at 
this altitude, there is some suggestion, on a finer scale than is shown in 
Figure 5.1, of a slowing down of the upward ion drift velocity that may be 
responsible for the convergence of ionization. At 17:40 AST the layer is quite 
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Figure 5.1 Vertical ion drift velocities (ms ) deduced from the horizontal components 
of the ion drift velocity measured on 16 July 1974 by incoherent-scatter 
radar . 
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intense again and begins to descend. This is probably due to the strong 
convergence level which is descending into the region. 

The fact that there is no convergence in the vertical ion velocity 
measured by the incoherent-scatter radar is evidence of the predominance 
of metallic ions in sporadic-# layers. In a diffusion-limited layer, the 
convergent velocity of the ions is exactly balanced by a divergent velocity 
due to diffusion. Thus, the mean velocity, which is measured in the 
incoherent -scatter experiment, is equal to zero. 

The vertical ion drifts necessary to maintain a diffusion-limited layer 
can be found using equation (2.22). This has been done for the profiles 
shown in Figures 3.10 and 3.14. The respective shears in the vertical ion 
drift velocities are 0.2 and 1.0 m s 1 km *. Using values for from 
Wright and Fedor [1970] , the vertical drifts would result from vertical 
shears in the east-west neutral wind of 13 and 18 m s" 1 km" 1 , which 
compare closely with the values observed by Rosenberg [1968] as illustrated 
in Figure 2.2. This analysis can only be valid in regions where the 
production and loss processes of the predominant ions involved are negligibly 
small, as in the case of metallic ions. 

5.3 Altitudes of Sporadic-E Layers 

5.3.1 Combined effects of gravity waves and tides. The magnitudes of 
tides and gravity waves are roughly the same [Hines t 1963], However, the 
persistence of the tidal winds and the endom nature of the gravity-wave 
spectrum indicates that long-term features of sporadic-# layers should be 
dominated by the tidal modes. This is supported by the observations of the 
rate of descent of intense sporadic-# layers which is at most a few km per 
hour (Figure 5.2). 

It is possible that, on some occasions, gravity-wave modes of comparative 
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intensity would interfere with the tidal modes, affecting the behavior of 
the sporadic-# layer. Chimonas [1971] maintains that since the horizontal 
extent of tidal winds is much larger than that cf sporadic-# layers, and 
since these winds are very regular, a resonance of the tidal modes with 
gravity waves is necessary to produce the sporadic nature and the required 
dimensions of sporadic-# layers. 

Interference between tides and gravity waves would also have a pronounced 
effect on the vertical position and the local ionization density of a 
sporadic-# layer. For example, destructive interference caused by an active 
gravity-wave mode would prevent a layer from forming at a point where there 
would ordinarily be vertical convergence of the ionization by a shear in the 
tidal wind. If the phase of the gravity wave were descending .nore rapidly 
than that of the tidal mode, points of strong convergence could appear and 
disappear. This would be further complicated by the direction of propagation 
of the gravity wave with respect to the tidal winds. 

Figure 5.3 illustrates vertical ion drift due to the interaction of two 
waves with different phase velocities. The sinusoidal east-west wind has 
constant phase and a wavelength of 10 km while a north-south sinusoidal wind 
has a wavelength of 12.5 km and has a downward phase velocity. Layering of 
the ionization would occur at points of convergence where vanishes. In 
the sequence illustrated in Figure 5.3, the convergence of ionization at 
123 km ceases between the fifth and sixth curve. The actual structure of the 
winds is never simple, but the model is used to demonstrate what may be 
occurring in the data shown in Figure 4.29. A layer is observed to form 
at an altitude of 123 km. This layer remains at a constant altitude for a 
period of time and then begins to descend, at which time the convergence of 
ionization seems to have ceased. In the situation represented in Figure 5.3, 
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Figure 5.5 Vertical ion drift caused by two superimposed sinusoidal wind systems: one in the 

east-west plane with a vertical wavelength of 10 km and a constant phase, the other 
with a vertical wavelength of 12.5 km and a downward phase velocity. I he phase of 
the second wave changes by 2n radians in the sequence shown. 
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this effect is enhanced by the fact that the vertical transport of ionization 
is dominated by the north-south winds at the higher altitudes, and by the 
east-west at the lower altitudes. 

5.3.2 Components of tidal modes. The relative importance of the east- 
west and north-south velocity components can have an effect on the altitude 
of a sporadic-# layer when the only important winds are due to a regular, 
slowly-descending tidal mode. Chimonas [1973] showed that a layer that is 
initially located near a node in the r.orth-south wind will descend at 
essentially the phase velocity of the node until the effect of the east-west 
wind becomes appreciable. Within a relatively short vertical distance in 
the neighborhood of 130 km, the convergence due to the north-south wind shear 
fails and the layer is "released" to descend rapidly to the node in the 
east-west wind. This is a more compelling model than that of interference 
with gravity waves described in the preceding section, since it is unlikely 
that any isolated wave will dominate the gravity-wave spectrum to the 
necessary degree. However, it seems to require that sporadic-# layers descend 
steadily above 130 km, and below an east-west node near 115 km, with the 
observation of a well-defined layer being improbable between these altitudes. 
In fact, the layer in Figure 4.29 seems to have descended in a discontinuous 
manner to the 123 km level before stepping to 10a km. 

If the release of a layer from a north-south node is indeed as abrupt 
as is indicated in the work by Chimonas [1973], it could possibly account 
fo 1 " the behavior of the sporadic-# layer observed on the night of 21 July 1975 
and shown in Figure 4.37. The layer would have initially been located at a 
node in the north-south wind at 127 km. Then, possibly as a result of a 
horizontally-propagating gravity wave, the altitude of the left side of the 
scan decreased below the point where the layer was released from the 
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north-south wind shear. The ionization in this part of the layer then 
descended to the next point of convergence, which was at 107 km. At the 
same time, the ionization in the part of the layer at the right side of the 
scan was lifted by the action of the wave, and remained in the region where 
the north-south shear was effective. 

5.3.3 Standing waves. As a wave propagates upward through the 
atmosphere, it will eventually reach a point where, because of changing 
temperature and background winds, the vertical wavelength becomes very 
large and the wave is reflected [Pitteway and Hines 3 1965] . This condition 
could set up standing waves, with nodes separated by a few kilometers and 
remaining at a relatively constant altitude. The effect of standing waves 
on the ionization may explain evidence that sporadic-!? layers form in 
"preferred" altitude intervals [ Seddon 3 1962; Smith and Mechtly 3 1972]. If 
there were a fairly consistent source of atmospheric waves which were 
reflected by a similarly consistent mechanism above the region of sporadic-# 
formation, nodes would form in the standing wave pattern near certain 
altitudes, leading to the apparent preferred altitudes of the sporadic-# layers. 

Standing waves may also be responsible for the discrete levels of layer 
formation where there is no relative motion between the layers, as illustrated 
in Figures 4.15 and 4.29. The former example shows two well-defined layers 
observed on 19 January 1974 between 14:21 and 14:32 AST. The two layers 
appear to be moving with the same velocity relative to the incoherent-scatter 
radar, and have trailing edges that coincide with each other. Neither layer 
was detected during the scan following the ones illustrated in the figure. 

This is interpreted to show that they are not two merging layers controlled by 
different atmospheric waves, but are parts of the same dynamic system. 

5.3.4 Quasi-stagnation levels. A gravity wave propagating through the 
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ionosphere will collect metallic ions near its velocity nodes, causing them 
to move downward with the phase velocity of the wave [Chimonas cmd Axford , 
1968]. This phenomenon has been called the "corkscrew" effect. As a result 
of the altitude dependence of the ion-neutral collision frequency, levels 
exist where the vertical ion motion can be shown to vanish independently of 
the phase of the wave inducing the motion resulting in "quasi-stagnation 
levels" [Hooke, 1971]. 

Quasi-stagnation levels may have a direct bearing on the discontinuous 
nature of the descent of a sporadic-^ layer through the lower E region 
(see Figures 4.29 and 4.37). The spectrum of gravity waves is expected to be 
quite broad, with each Fourier component subject to its own stagnation level. 
This could lead to a series of levels at which ions following the downward 
phase of a wave are trapped. It could also mean, however, that no components 
of the wave would dominate the spectrum to a great enough extent for a 
stagnation level to affect the vertical motion of the ions. 

The levels of stagnation can be found by solving equation (2.17) for the 
vertical ion drift. It is convenient to fix the coordinate system with 
respect to the direction of propagation of the wave, such that 


u =0 

y 


(5.5) 


Thus, 


v = — — s- [(p.r +r r )u + (p. 2 +r 2 )u ] 

3 - 2 L y x z x r z J z J 

Up. * 


(5.6) 


Hines [1960] showed that for low frequency internal waves 
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where k and k are the horizontal and vertical wave numbers, respectively. 
X z 

If this is substituted into equation (5.6), it is found that the vertical 
ion velocity vanishes when 



Hooke [1971] treats the case of an east-west (geomagnetic coordinates) 
wave propagation, so that r =0, and the stagnation level occurs where 

X 



The requirement that p . is real and positive restricts the possible values 

of k and k . 
x z 

This equation can be solved if the wave numbers are known. In the 
analysis by Hooke, a value of k a /^ z = 0.1 is used, corresponding to a wave 
with a period of approximately one hour. At Arecibo, the magnetic dip 
angle is 50 deg, giving p. equal to 6.3 and 0.093 for the two solutions. 

"I* 

The corresponding heights of the stagnation layers, from a table of p. 
versus height by Wright and Fedor [1970], are approximately 106 and 155 km. 

A similar calculation for conditions at Wallops Island (magnetic dip = 70 deg) 

gives p. = 3.1 and 0.28, corresponding to altitudes of 110 and 134 km. 

z 

A major difficulty in relating this theory to the observed sporadic-E 
layers is the location of the expected stagnation levels. The period of the 
wave described above is relatively short with respect to the usual rate of 
descent of sporadic-E layers. Yet the levels of stagnation, at 106 and 155 km, 
are respectively too low and too high to be applied to the observations at 
Arecibo. Examination of equation (5.9) shows that for waves with longer 
periods, the respective stagnation levels would occur at even lower and 
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higher altitudes. 

5.4 Boundaries of Sporadic - E Layers 

An interesting feature observed in the sporadic-# layer of 17 July 1975 

(Figure 4.29) is the well-defined trailing edge of the layer. Electron-density 

contours observed during the three scans of the radar beam as the edge was 

passing through the viewing region are shown in Figure 5.4. In this figure, 

4 -3 

the electron-density contours are on a linear scale, at intervals of 10 cm 

4-3 

Contours at intervals of 5 * 10 cm are indicated by thick lines. Note 
that the vertical and horizontal distance scales are not equal. 

The abruptness of the edge of this layer is difficult to explain in 
terms of atmospheric waves. When compared with the normal horizontal scale 
sizes of the ionosphere, it would appear as a discontinuity. The column 
electron density is also greater where the layer is present, indicating that 
the edge of the layer is not a result of a change in the vertical-convergence 
mechanism. Anderson and Barth [1971] made a similar observation in an 
experiment which measured the dayglow from the ionized and neutral magnesium 
of a sporadic-# layer. The rocket carrying the instrumentation apparently 
passed below the edge of a sporadic-# layer. Within a horizontal distance 
of about 6 km the 2800A dayglow, which is an indication of the column 
density of Mg + , increased from 0.5 to 5.0 kR. 

A clue to the mechanism sustaining this large horizontal gradient is 
found in the less intense layer that is located on the right side of the 
apparent discontinuity and approximately 4 km lower in altitude than the 
intense layer. The structure of the electron density shown in Figure 5.4 
suggests that boundaries exist in the lower # region between masses of air 
with different static and dynamic properties, as in tropospheric frontal 
systems, with the wind components parallel to, and on either side of this 
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Figure 5.4 


Electron density observed on 17 July 1975 (see Figure 4.29) 
Contours are at intervals of 1 x 10 4 cm"3, with contours at 
multiples of 5 x 10 4 cm - ^ darkened. 
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"front", exhibiting different characteristics. For example, let the boundary 
at the "front" be a vertical east-west plane. Suppose wind systems on either 
side are similar, with the exception that the position of the null in the 
east-west wind is 4 km lower on the north side of the boundary than on the 
south. This is shown schematically in Figure 5.5. The drawing on the left, 
(Figure 5.5a), shows the direction of the neutral winds with respect to the 
boundary and to the geomagnetic field. The position of the nulls in the wind 
profiles coincide with the sporadic-^ layers, and are indicated by Es. 

The ion drift velocities in the north-south plane are shown in Figure 5.5b. 
It is seen that horizontal ion motions in the region of the boundary tend to 
oppose the horizontal diffusion of ionization from the edges of the layers. 

This would produce horizontal gradients, as seen in Figure 5.4. 

The required horizontal ion drift is not great. If the layer is assumed 

to be diffusion limited so that equation (2.22) applies the ion drift which 

would maintain the observed gradient can be found by the methods described in 

2 -1 

the previous section. Assuming an ambipolar diffusion coefficient of 70 m s 
[Kaiser et al. 3 1969] it is found that a horizontal ion drift of 0.02 m s * 
would be sufficient to maintain the edge of the layer for steady-state 
conditions. Bounded wind systems may also be responsible for the noticeably 
different dynamics in the two sides of the region observed on 21 July 1975 
and shown in Figure 4.32. 

There is evidence of rapid changes in horizontal velocities in data from 
incoherent scatter used in this study (see Figure 4.24) and in data from the 
University of Illinois meteor radar [ G . C. Hess t personal communication, 1976], 
However, the poor time resolution of these wind-measuring techniques precludes 
the observation of small-scale discontinuous effects of the type mentioned 
here. 




Figure 5.5 Sketch of possible configuration of winds producing 

the electron-density distribution shown in Figure 5.4. 
(a) Direction of sinusoidal east-west horizontal 
winds, (b) Ion drift velocities and locations of 
sporadic-E layers. 
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5.5 Wave-like Structure in Sporadic-E Layers 

The existence of wave- like structure in the ionosphere has been known 
and studied for many years [ Ratcliffe , 1959], Experiments using ionospheric 
sounders operating at a fixed frequency and using spaced receivers show 
apparent random velocities superimposed on a relatively constant background 
motion associated with tidal modes [Hines, 1974], The random velocity 
components observed by spaced Aonosondes are possibly the effect of dominant 
gravity-wave modes propagating through the region where the radio wave was 
reflected. The observed random "velocities" would be the superimposed phase 
velocities of a near continuum of waves, and not an actual bulk motion of 
ionization. 

Rapid changes in the observed "velocities", and even changes in 
direction, might be observed as a result of the introduction of one wave train 
or the loss of another [Hines, 1974, postscript to Paper 33], The 
irregularities may be in the form of changes in the altitude of the 
reflecting level, or actual changes in the ionization density. 

Various mechanisms have been suggested for the horizontal convergence 
of ionization by gravity waves. These include changes in the ion chemistry 
resulting from changes in neutral density and temperature [Hooke, 1969], the 
effects of vertical air motions and the geomagnetic field on ions [Wright and 
Fedor, 1969], and resonance between tidal and gravity-wave modes [Chimonas, 
1971]. 

Wave- like motions of sporadic-E layers were observed on several occasions 
during this study, ranging from periods of a few hours to a few seconds. 


Because of the nature of the data and the inability to resolve motions in 
more than one horizontal direction, it is impossible to separate variations in 

■f-irnp from cnatial vari at i on c and thprpfnrp rliffirnlt tn n«» frenuenrie*; 



and wave numbers. Where possible, wavelengths or periods are estimated and 
compared to calculations by Hines [1960] of allowed gravity-wave modes. 

Electron-density profiles in the vicinity of a sporadic-# layer are 
plotted at 5 min intervals in Figure 5.2. Evidence of short term variations 
in altitude is seen in parts of the data. However, the dominant oscillation 
seems to have a period of approximately two hours. Assuming the atmosphere 
at the sporadic-# layer to be essentially at rest with respect to the ground, 
the allowed gravity-wave modes of this period, according to the calculations 
of Hines [1960], must have horizontal wavelengths of between approximately 
100 and 1000 km, and vertical wavelengths between 10 and 50 km in order to 
have propagated from the lower atmosphere and not be dissipated by 
viscous effects. 

Another example of wave-like structure in sporadic-# height variations 
was observed on 26 July 1975 (Figure 4.43). In this case, the horizontal 
wavelength was comparable to the distance through which the radar beam was 
scanned, approximately 20 km. According to the calculations of Hines [1960] 
and assuming no net motion of the layer, these waves could propagate with 
periods of between 10 and 160 sec and vertical wavelengths of from 1 to 50 km. 

Wave- like structure is evident in the form of electron-density 
fluctuations in the data from the single-pulse experiment on 16 July 1975. 
Figure 5.6 shows contours of electron density as a function of altituc’e and 
azimuth. The altitude resolution is about 1.2 km. Contours of electron 
density greater than 1.5 * 10 3 cm -3 are accented and regions where the average 
electron density is greater than 1.75 * 10 3 cm 3 have been darkened to 
emphasize variations within the sporadic-# layer. The main component 
appears to have a wavelength of 5 km. According to calculations by Hines 
[1960], a wave of such short wavelength could have originated in the lower 
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atmosphere only if the period is ajout 10 sec and the vertical wavelength 
is about 2 m, which are inconsistent with the data. A possible explanation 
lies in a source of the gravity wave in the middle or upper atmosphere. 

5.6 Irregular Structure Associated with Critical Levels 

5.6.1 Introduction. The main purpose of this study is to describe the 
horizontal structure of sporadic-# layers. This has been done in the 
above sections for small -amplitude perturbations in terms of atmospheric 
gravity waves. However, none of the mechanisms described can satisfactorily 
account for the very localized regions of enhanced ionization shown in 
Figures 4.21 and 4.32. This section will demonstrate that a plausible 
exp 1 anation of these intense irregularities lies in their proximity to 
critical levels of internal gravity waves. 

5.6.2 Critical levels. A wave propagating through a region where there 
is a shear in the background wind is unable to propagate through a level 
where the Doppler-shifted frequency of the wave is zero in the rest frame 

of the moving atmosphere. This is known as a critical level, and may be a 
cause of instability and turbulence. 

Some of the characteristics of critical levels were derived by Hines 
[1968], beginning with first order equations, including the vertical 
variation of the scale height, and making no assumptions about the vertical 
dependence of the wave parameters. The perturbation pressure, density, and 
horizontal velocity can be eliminated in this system of equations to give a 
second order differential equation describing vertical variations of the 
vertical perturbation velocity, w z . In the limit as the Doppler-shifted 
frequency (D) of the gravity wave goes to zero, this equation becomes 

cu R 

~ 7 7 + "77 T7 “ 0 

dz a (z-zj 
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(5.10) 
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where co^ is the Brunt-Vais3ia frequency, given by 

u>/ = [y - 1 + y(dH/dz)]g/yH (5.11) 

y is the ratio of specific heats, U is the pressure scale height, and ft 
has been expanded as 

ft s a) . k x [w Qx + a (z-z^ + ...] (5.12) 
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A solution of equation (5.10) can be found in the form 


U = A (z-z ) 
z o' 


(5.13) 


where 


v = 1/2 ± i 


J B_ l 
2 “ 4 
a 


(5.14) 


2 2 

The ratio /a is known as the Richardson number, and is an indicator of 
stability of a gas subject to a velocity shear. 

It can be seen that as a gravity wave approaches a critic.'! level, the 
vertical amplitude of the wave goes to zero as 




(5.15) 


By substituting the solution for w^ 3 the perturbation amplitudes of pressure, 
density, and horizontal velocity near the critical level are found to vary 
inversely with the square root of the distance from the critical level. The 
shears in the velocity of the neutral wind also become very large near the 
critical level, varying as 


dw 

X 




(2-Z 0 ) 


-3/2 


(5.16) 



{ 

■ 


i 


191 




~ ( z - z 0 )~ h (5.17) 

Under most conditions, the effect of critical levels on the atmosphere 

is not significant. A spectrum of all gravity waves at a point in the 

atmosphere would probably not be dominated by any particular frequency. It 

has been suggested, however, that if a wave were excited for an extended time 

by a reasonably coherent source, enough of the wave could become critical at a 

large wind shear to produce a measurable effect [Hines j 1968] . 

5.6.3 Effect in sporadic-E layer formation. Whitehead [1971a] suggested 

that, due to the increased vertical wind shear, intense sporadic-# layers 

would be most likely to form at or near the height of a critical level. 

Equation (5.17) suggests that there would also be an increase in the horizontal 

wind shear. A large horizontal shear in the neutral wind at a sporadic-# 

layer would cause the horizontal convergence or divergence of the ionization. 

This mechanism may be the source of intense electron-density irregularities 

with limited horizontal extent, such as was observed on 28 January 1974 

(Figure 4.21) and 19 July 1975 (Figure 4.32). 

In the following analysis, vertical and horizontal ion drift velocities 

will be found which, under steady-state conditions, would form an irregularity 

with the peak electron density observed on 19 July 1975. In this case, the 

5 -3 

vertical shear alone must produce a peak electron density of 8 * 10 cm , 

and the superposition of the vertical and horizontal shears produce a peak 

6 — 

electron density of 6 * 10 cm in order to account for the observed 
irregularity. 

The peak ionization density in a diffusion-limited layer can be 
calculated using the continuity equation (2.22) when the vertical ion drift 
velocities are known by 
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ln{n/n o ) = J ~jj dz 


(5.18) . 


First consider the effect of vertical convergence in the absence of 

horizontal convergence. In particular, v will be assumed to be sinusoidal 

z 

and V is assumed to increase exponentially with altitude. Integrating 
equation (5.18) then gives 


n = n exp 


y /(l/tf)sin k (z-z ) - k cos k [z-z ) 
zm 2 O z 2 o' 


(i/n + k_ 2 


(5.19) 


where is the velocity amplitude of the wave, and is the unperturbed 
metallic ion density. The ratio of the electron density at the peak of 

the laver to n is then 

o 


v k 
zm z 


D td/ff 2 ) * k, 2 ] 


(5.20) 


Mass spectrometer measurements show the background metallic ion density 
to be of the order of 10^ cm 3 [A -avcisi 3 1973; Zbinden et al. 3 1975]. 

2 -1 

Assuming a vertical wavelength of 10 km in v and a value of D of 230 m s 

z 

at 115 km [Kaiser et al. 3 1969] (the altitude of the sporadic-# layer), 

with a scale height of 7 km [CIRA, 1972], the required maximum vertical ion 

5 -3 -1 

drift to maintain a peak ionization density of 8 * 10 cm is 1.0 m s 

Next consider the case of vertical and horizontal convergence acting 
together. The continuity equation (2.22) is solved Tor this two-dimensional 
case by the method of separation of variables (x and z ) . For v = 

£C 

" v xm sin [** (*-*,)] the result is 


n r v k y 1 

p zm z xm 

n = ex P I 2 T + t 

o U [(1/0 + kj] D k J 


(5.21) 
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In order to explain the irregularity in Figure 4.32, the value of 

3 

n /n must be of the order of 6 x 10 . If the vertical ion drift is assumed 

P o 

to be the value used previously and a horizontal wavelength of 50 km in 
- v is assumed, a maximum horizontal drift of only 0.05 m s *, under steady-state 

conditions, is required. 

At Arecibo, where the magnetic dip is 50 deg, east-west horizontal winds 
will produce vertical and north-south ion drifts of approximately the same 
magnitude. For a neutral wind of 100 m s * at an altitude of 115 km, each 
component of the ion drift velocity would be approximately 20 m s *. 

Horizontal convergence occurs only for the case of a shear in the east-west 
wind. The formation of irregularities by a horizontal shear in the north-south 
wind would be prohibited by the polarization electric field, in the same 
manner that sporadic-E layers are prohibited from forming at the magnetic 
equator. This is consistent with the suggestion that the irregularity shown 
in Figure 4.32 is elongated in the east-west direction. 

The ion drifts derived in the above analysis are much less than those 
observed, and the widths of the corresponding layers are much greater than 
observed. The purpose of the preceding discussion is to show that the 
horizontal ion velocities required in such a model are not unreasonable, 
especially when considered to be the .result of a gravity wave approaching 
a critical level. In his derivation of critical levels, Hines [1968] noted 
that their effects should be enhanced in a region where there is a large 
wind shear. Even if there is no dominant Fourier component, enough of the 
spectrum can go critical in a relatively small height interval to produce 
a pronounced effect on the ionization structure. 

Further work is needed in evaluating the effect of critical levels on 
sporadic-E ionization, as the observed electron-density irregularities are 
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-n 
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too highly enhanced for the approximations of the linear perturbation theory 
used here to be valid. Howe\er, on the basis of linear theory, it is 
concluded that horizontal convergence of ionization with sufficient magnitude 
to account for observed irregularities in sporadic-^ layers can occur near the 
height of a critical level. 

5.7 Irregular Structure Caused by Fluid Instabilities 

An example of a fluid instability is seen in Figure 4.29. The very 
pronounced structure seen on the top side of the layer between 04:39 and 
05:12 AST is characteristic of the Kelvin-Helmholtz instability which occurs 
in regions of large velocity shears [Corcos and Sherman , 1976; Patnaik et al. 3 
1976]. The instability is driven by oscillations in the flow on either side 
of the shear having a horizontal phase difference of tt/ 2. The boundary layer 
at the maximum shear, i.e. the location of the sporadic-E layer, then begins 
to exhibit wave-like characteristics which grow in amplitude until destroyed 
by turbulence. 

This sequence of events is evident in Figure 4.29, and appears to be 
initiated by a wind pattern associated with the descent of the upper layer. 

At the approximate time that the upper layer would have reached the altitude 
of the lower layer, irregular wave-like structure begins to develop. The 
peak ionization density also increases, indicating that the shear in the 
neutral wind has increased. A more detailed plot of the layer at this time 

is shown in Figure 5.7. The electron-density contours are drawn at intervals 

4-3 4-3 

of 10 cm , with contours at intervals of 5 * 10 cm accented. 

The amplitude of the wave continues to grow until about 05:12 AST, 
when the layer thickens suddenly, probably at the onset of turbulence. 

Electron density during this time is shown in Figure 5.8. The decrease in 
the peak electron density of the layer at this time was also detected by the 
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ionosonde, Figure 4.30, whose field of view was separated from that of the 
radar by about 30 km. 

5.8 Irregular Structure Caused by Plasma Instabilities 

A theory that in the past has been advanced as an alternative to the 
formation of sporadic E by wind shears is that of instabilities in the 
ionospheric plasma [Clermow and Johnson, 1959; Tsuda et al ., 1966; Tsuda and 
Sato , 1968], The non-blanketing sporadic E that occurs at the magnetic 
equator is caused by plasma instabilities in the electrojet [ Farley , 1963a, b; 
Reid , 1968; Prakash et al., 1971]. However, the existence of plasma 
instabilities in midlatitude sporadic E has not been established. 

Whitehead [1967, 1968, 1970b] developed the theory of a gradient 
instability, assuming midlatitude conditions, that could occur on a 
Gaussian-type layer of ionization. He found that when the equation of 
motion and the continuity equation were solved together, assuming small 
perturbations in ionization density and electric potential, irregularities 
would form near the peak of the layer and persist for tens of seconds. The 
fastest growing irregularities were found to be partially aligned with the 
magnetic field, which is in agreement with some observations of small sporadic-E 
irregularities [Gooduin and Thomas, 1963; Gooduin, 1965]. Atkinson [1973], 

showed the irregularities to have time constants of 1 to 170 s and periods of 

-2 -1 
from 10 to 17 s, corresponding to phase velocities of 5 to 500 m s 

Whitehead [1968] concluded that this instability would be unable to form 

sporadic-E layers, but could be responsible for irregularities with scale 

sizes of a few tens of kilometers within the layers. He showed that the 

maximum ratio of the perturbation ionization density to the unperturbed 

ionization density was about 0.1 or 0.2. At greater values, the perturbation 

electric field would exceed the ambient electric field, causing the 
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irregularities to self-destruct [Kato 3 1963]. 

The intense irregularities perhaps most likely to be attributed to an 
unstable plasma are shown in Figures 4.21 and 4.32. In each case the 
irregularity was observed for more than five minutes before drifting out of 
the field of observation, and each had a peak ionization density that was 
greater than the peak density in the surrounding sporadic-# layer by at least 
a factor of two in contrast to the maximum relative perturbation amplitude of 
0.2 mentioned above. Another consequence of the gradient instability discussed 
by Whitehead [1967] is that if a layer becomes less than about 1.5 km thick, it 
will be destroyed by turbulence. There was no evidence in this study of the 
breaking up of localized irregularities by turbulence. 

5.9 Irregular Structure Attributed to Meteor Trails 

5.9.1 Irregular source function. Preceding sections have considered the 
formation of irregularities in sporadic-# layers by the horizontal convergence 
of ionization by the action of the neutral wind. An alternative mechanism for 
the formation of localized enhancements of electron density is provided by the 
irregular source of metallic ions as a consequence of their origin in meteor 
trails . 

First the immediate effect of a meteor trail on the incoherent-scatter 
radar observations is described. Then a mechanism is proposed which allows the 
convergence of the ionization of a meteor trail into a localized region of high 
electron density within the sporadic-# layer. 

5.9.2 Individual meteor trails. In Chapter 2 it was noted that a 

significant number of metallic ions are left in the meteor trail when a 

meteoroid enters the atmosphere. A representative initial linear density of 

14 -1 

meteoric ions is 10 cm [Baggaley and Cummack 3 1974; Poole and Nicholson, 
1975], This is only a small contribution to the ionization already present 


in a sporadic-# layer. If such a meteor were to pass through the beam of the 

incoherent-scatter radar, it would raise the average electron density measured 

5 -3 

in a particular 600 m height gate by about 3 x 10 cm . An apparent meteor 
trail observed on 21 July 1975 is shown in Figure 4.5. The peak electron 
density is about 1.5 x 1 «j^ cm The detection of an individual meteor trail 
is a rare event and, in any case, even more rarely would occur in a sporadic-# 

1 ayer . 

5.9.3 Convergence of a meteor trail. While the ionization density left 

initially in a meteor trail is not insignificant, it is not nearly enough to 

account for the irregularity with a peak electron density possibly as great 
6 “3 

as 6 x 10 cm observed on 19 July 1975. A possible explanation lies in the 
convergence of a meteor trail by the wind shear mechanism. The trajectories 
of ions under the influence of both an external magnetic field and collisions 
with neutral particles in motion was derived in Chapter 2. It was shown that 
ion drift velocities are described by 

V = — [p 2 « + pu x r + (u«r)r] (5.22) 

1+p 

An examination of the effect of various wind-shear configurations was 
made by letting the winds act on a line of ionization as described by equation 
(5.22). The effects of both the orientation of the line of ionization and the 
orientation of the winds with respect to the geomagnetic field were determined. 

For most configurations the effect of the wind shear is to spread the 
ionization out horizontally as it converged onto a horizontal plane at the 
null point of the wind velocity. However, it is found that the effect on 
particular configurations is to converge the ions into a small, relatively 
dense region of ionization. This is illustrated in Figure 5.9, showing the 
initial and final configurations of an ionized meteor trail initially 15 deg 
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from horizontal and at an azimuth of 230 deg. The wind has an east -west 
component only, and is represented by a sine wave with an amplitude of 100 m s 
and a wavelength of 20 km. The ion-neutral collision frequency is assumed to 
decrease exponentially with height in such a way that p equals one at the 
altitude of the wind shear. 

This analysis neglects many processes such as diffusion, production, 

and recombination. It does, however, demonstrate that the wind system is 

capable of converging ionized meteor trains into irregularities such as have 

been observed in the incoherent -scatter data. If the ionization from a 20 km 

14 -1 

meteor trail with a line density of 10 cm were converged by this mechanism 

and sampled in one height gate by the radar, the average electron density at 

6 -3 

that height would be increased by 4.7- x 10 cm , which agrees very well with 
the observation of 19 July 1975. 

The time allowed for the convergence of the meteor trail in Figure 5.0 

is 20 min, although the process is essentially complete after 10 min. 

2 -1 r , 

Assuming an ambipolar diffusion coefficient of 35.6 m s [Katser et al. t 
1969] at an altitude of 102 km, the two-dimensional diffusion of a line of 
ions with an initial half-width of 1 m would result in a half-width of about 
400 m after 20 min. The effect of diffusion is lessened somewhat by the 
restriction of movement perpendicular to the geomagnetic field. The duration 
of partial reflections from sporadic E at frequencies very much greater than 
the blanketing frequency seldom exceeds the 20 minutes used in the above 
example. 

5.9.4 Correlations with meteor activity. Creation of irregularities in 
sporadic-# layers from ionization in individual meteor trails could account 
for the sporadic nature of the observed partial transparency, as these 
irregularities would form only under very specific conditions. Their 
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dependence on the action of the neutral wind shear would also give rise 
to the observed correlation between large frequency ranges of partial 
reflections and intense sporadic-E layers. 

Two factors apparently not accounted for by this mechanism are the 
diurnal variation of ftEs and the lack of correlation with reteor activity. 

Top frequencies of 13 MHz, as in the case of the 19 July 1976 observation 
(.Figure 4.28), and even greater, are observed at times in the daytime, but 
not at night. If convergence of the metallic ions left in meteor trails 
were the cause of irregularities, they would be as likely to form at night 
as in the day. 

Many researchers have attempted to correlate sporadic-E occurrence 
with meteor activity with varying degrees of success [ Fiocco 3 1965; Kotadia 
and Jani, 1967; Wright, 1967; Whitehead , 1971b]. Hedberg [1975, 1976] showed 
that a dependence of fbEs on meteor rate can be shown if a time delay on the 
order of 6 hours is included. However, in the same study, no similar 
correlation of ftEs or foEs with meteor activity was found [/?. Hedberg , 
personal communication, 1975]. Whitehead [1970a] suggested that only certain 
meteor showers may be responsible for the metallic ions present in sporadic-E 
layers. Shower meteors are travelling in essentially parallel trajectories 
and would, if properly aligned with respect to the wind shear, produce many 
irregularities. Trtskovd [1974] showed that the relationship between the 
occurrence of sporadic E and meteor activity might be better understood when 
considered together with changes in the zonal wind system. 

Many of the attempts to correlate the occurrence of sporadic E with 
meteor activity have considered only the shower meteors. However, the majority 
of the metallic ions is introducer 1 ’nto the atmosphere by randomly-oriented 
sporadic meteors. Any wind shear which converges ionization vertically will 
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also converge into localized irregularities the fraction of sporadic meteors 
that are properly aligned with respect to that wind shear. 

It is thus concluded that on some occasions, small-scale irregularities 
in electron density probably contain ionization from individual meteor trails 
which has been converged by the shear in the neutral wind. 

5 . 10 Siomary 

In this chapter, observations of sporadic-# layers by incoherent-scatter 
radar have been discussed. Several features have been noted in the data that 
support the wind-shear mechanism of layer formation. Ion drift velocity 
measurements, when interpreted in terms of diffusion-limited metallic ion 
layers, showed the altitude of a sporadic-# layer to be near a point of 
convergence in the vertical ion drift. Preferred altitudes of occurrence of 
sporadic-# layers, as well as changes in altitude, are found to be consistent 
with the effects of wind shears caused by gravity waves and tides. The large 
horizontal gradient at the edge of a layer is indicative of a horizontal shear 
in the neutral wind, and is interpreted as being located at a boundary between 
wind systems. 

The sporadic-# layers often contained pronounced small-scale structure, 
especially at times when partially transparent echoes were observed by the 
ionosonde. It was shown, based on several observations of wave-like structure, 
that internal gravity waves can be responsible for significant perturbations 
in the altitude and density of a sporadic-# layer. A sporadic-# layer at an 
altitude near a critical level of a gravity wave may be acted upon by highly 
enhanced vertical and horizontal shears in the neutral wind, creating small 
regions of high electron density. Irregular structure is also found to be 
generated by fluid instabilities in regions of large wind shears. 

Small regions of enhanced electron density may also be an effect of the 


irregular nature of the source of metallic ions. It was shown that under 
specific conditions the ions in a meteor trail will be converged by a shear 
in the neutral wind into a relatively small irregularity at the center of a 
sporadic-E layer. 

In the observations of sporadic-E layers presented here, the small-scale 
horizontal structure has been emphasized. It is concluded that these 
small-scale irregularities are consistent with wind-shear theory when the 
combined effects of gravity waves, tides, and the irregular deposition of 
metallic ions are considered. 


6. SUMMARY AND CONCLUSIONS 


6 . 1 Summary 

The objective of this study is to determine the source of the partially 
transparent echo from sporadic-# layers. The investigation is mainly con- 
cerned with two sporadic-# models that have been proposed to account for 
these reflections. In one model, the reflections are from the large elec- 
tron-density gradients at the upper and lower boundaries of sporadic-# 
layers. The other holds that the radiation is reflected by small regions 
of enhanced ionization imbedded within the sporadic-# layer. 

The basic theory of the formation of sporadic-# layers by wind shears 
is developed. Large wind shears are nearly always present in the region of 
the atmosphere where sporadic-# layers are observed, while, as the name 
suggests, the occurrence of sporadic-# layers is erratic. The role of 
metallic ions in sporadic-# formation is iscussed and a model developed of 
the diurnal, behavior of the metallic ions and atoms. It is found that a 
significant amount of metallic ionization is introduced into the ionosphere 
directly as a result of ablation from a meteoroid entering the atmosphere. 

The importance of reflections of radio waves from electron-density 
gradients associated with sporadic-# layers is discussed in Chapter 3. A 
full-wave analysis is developed whereby reflection coefficients of sporadic-# 
electron-density profiles measured by rocket-borne probes is calculated. 

It is found that the frequency range of partial transparency is small both 
for layers causing partially transparent echoes on ionograms and those where 
no partially transparent echo is detecteo. It is thus concluded that large 
ranges of partial transparency in the echoes cannot be caused by electron- 
density gradients associated with the sporadic-# layer. 
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The small-scale horizontal structure of sporadic-E layers was investi- 
gated by a series of observations at the Arecibo Observatory in Puerto Rico. 
The observations were made during January, 1974, and July, 1975. Two basic 
experiments were used, one using phase-coded pulses to record electron-den- 
sity profiles with a range resolution of 600 m, the other using multiple 
pulses to measure ion drift velocities. Highlights of each observation are 
summarized in Chapter 4. It is evident that sporadic-E layers contain 
small-scale irregularities of sufficient intensity to account for the par- 
tially transparent echoes observed by the ionosonde. Periods of intense 
irregularities and observations of partially transparent echoes are found 
to be well correlated. 

Many interesting features observed in the incoherent- scatter data 
are discussed in Chapter 5. The observed structure of the sporadic-E 
layers is first considered with respect to motions of the neutral atmos- 
phere. Various configurations consistent with wind-shear theory are dis- 
cussed. Features indicating the presence of boundaries in wind systems 
similar to tropospheric fronts are also discussed, as well as the observa- 
tion of the growth of a Kelvin-Helmholz instability and its subsequent 
destruction by turbulence. 

The irregular deposition of metallic ions by individual meteors is 
also considered as a source of irregularities . It is found that under 
specific conditions, ionization in a meteor trail will be converged by 
actions of a neutral wind shear into a small irregularity at the level 
where a sporadic-E layer will form. 

The main conclusion of Chapter 5 is that sporadic-E layers are 
significantly influenced by strong Fourier components in the gravity-wave 
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waves are concluded to be the main source of the atmospheric 
motions affecting the small-scale structure of sporadic-# 
layers. 

(v) Intense localized irregularities indicate the proximity of 

critical levels in the internal -gravity-wave spectrum. Large 
horizontal -velocity gradients that occur near critical levels 
are suggested as the mechanism by which such irregularities 
are formed. 

6.3 Suggestions for Future Work 

Observations of neutral winds should be compared with ion drift veloci- 
ties determined by the incoherent-scatter radar. This would provide a 
better check on the predictions of wind-shear theory than was possible in 
the analysis of Chapter 5. The two measurements might be made simultaneously 
using, for example, vapor trails from rockets launched from a site near the 
Arecibo Observatory. 

Electron-density measurements are also needed that coincide with ion- 
drift-velocity measurements. In the multiple-pulse experiment described in 
Chapter 4, electron-density measurements were made onlv between ion-drift- 
velocity measurements, or during about 5 min of each 25 min period. Pulse- 
compression techniques can be developed to measure a> :ocorrelation functions 
to give ion drift velocities with the height resolution necessary for studies 
of sporadic-E layers. An experiment could then he conducted in which scat- 
tered-power profiles would be recorded with good time resolution at the 
same time that autocorrelation functions are integrated to give the statisti- 
cal accuracy required for the measurement of ion drift velocities. 

Further research is needed in determining source functions and concen- 
trations of metallic atoms and ions. The number of metallic particles h. 'ng 
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spectrum. Wave-like structure is evident in many of the observations. It 
is shown that the most intense irregularities observed in this study can 
be produced by gravity-wave modes nearing critical levels as they encounter 
the sporadic-E ionization. Strong horizontal shears may be present at 
altitudes near a critical level which would converge the ionization of a 
sporadic-E layer into the observed irregularities. 

6.2 Conclusions 

The following are the main conclusions arrived at through this study: 
(i) Gradient reflections are not important in ionosonde observa- 
tions of sporadic-E. 

(ii) Small-scale horizontal structure exists in midlatitude 

sporadic-E layers. The magnitude of the irregularities in 
electron density is strongly correlated with the observation 
of transparent echoes from the sporadic-E layer. Small 
regions of enhanced electron density provide an explanation 
for the occasions when ftEs » fbEs. 

Ciii) The features of sporadic-E layers observed by incoherent -scat- 
ter radar are found to be consistent with layers composed of 
metallic ions being acted on by the motions of the neutral 
atmosphere as described by wind-shear theory. The dominant 
neutral motions are gravity waves and tidal modes. 

Civ) The irregular structure observed in midlatitude sporadic-E 
layers is mainly determined by the motions of the neutral 
atmosphere. This is especially evident in the case of night- 
time sporadic-E layers where features in the structure at 
small electron densities can be observed. Internal gravity 
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deposited in the ionosphere and the fraction of these particles that are 
deposited directly as metallic ions needs to be determined. This should be 
done for both daytime and nighttime conditions. 

More incoherent-scatter observations of sporadic-E layers are needed. 
The data presented in Chapter 4 represents a very limited sample. The 
sporadic-E layer in each observation was found to have features different 
from those in other layers observed, suggesting that there may be other 
effects that have not yet been detected. 

The small-scale structure of sporadic-E layers should be considered 
in theoretical studies of atmospheric fluid instabilities and critical 
levels. In particular, the development of a non-linear theory which addres- 
ses itself to the observed irregularities is needed. The detailed informa- 
tion on the distribution of the electron density that is possible using 
pulse^compression techniques should be exploited in the study of non-linear 
processes associated with wind shears, as well as the onset of turbulence. 
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APPENDIX I 

The Effect of the Polarization Electric Field on 
the Formation of Sporadic-^ Layers 

When ionization is moved by the action of the neutral wind, the 
different charges and collision frequencies of electrons and ions cause 
them to move in different directions. Unless the plasma is infinite and 
uniform, this gives rise to a polarization electric field opposing their 
separation. Explicit expressions for the polarization electric field are 
developed here, following the work of Whitehead [1961], These expressions 
are applied to a model of a region of large wind shears to obtain an 
estimate of the magnitude of the polarization field. It will be shown that 
the field is small at midlatitudes and can be neglected. However, at the 
magnetic equator, the polarization electric field becomes large and prevents 
the formation of layers by vertical shears in the horizontal wind. 

The effect of the polarization electric field on the motion of the 
ionization can be found by solving equation (2.17) for both electrons and 
ions. Assuming no interaction between electrons and ions, with v. and v 
as solutions of equation (2.17), the current resulting from their differing 
motions could be written as 


3 = 


k e(v r v e ) 


( 1 . 1 ) 


where n is the electron density. 

It is assumed that in a sporadic-2? layer, any horizontal current caused 
by the relative motion of electrons and ions above the peak of the layer is 
complemented by another current in the opposite direction below the peak of 
the layer. Thus, 


ite - 


n (u -v )dz = 0 
X • X 
x e 



i 


Es 


Es 


( 1 . 2 ) 
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j dz = e [ n(v -v )dz = 0 (1.3) 

>es y >e s v € y e 


where the integration is over the sporadic-# layer. It is also assumed that 
there are no vertical currents through the layer, so that 


j - en (u —v ) = 0 (1.4) 

Z 2 • Z 

i e 


In the following derivation, only the effects of horizontal winds and 
polarization electric fields, designated by E, will be considered. The 
coordinate system will be oriented so the 3 is in the y-z plane, and the 
x-axis is toward the east. Thus, equation (2.17) can be written 


v = 
x 


(1+P 2 ) 


2 EE 

p „ ♦ P r w ♦ p f * r 2 f 
3 o o . 


(1.5) 


V = 


(1+p 2 ) 


. 0 E (p 2 +r 2 ) E T T E ' 

_ ,2„2. „ x it it u z 2 

-pr w + (p +r )u - r -r— + 2 — r 2 - + 

z x y y a b p b p b 


O J 


( 1 . 6 ) 


v = 


e r r e (p 2 +r 2 ) e 
pT w + t t w + r ~ + 2 2 


(1+p 2 ) L ^ * 2 V y B 0 P B o p B 


o J 


(1.7) 


The corresponding components of the electric current are found by 
subtracting the electron velocities from the ion velocities. £ will be 
expressed in terms of its components parallel and perpendicular to the magnetic 
field, so that 


e = r e + r e, 

y y 3 1 


( 1 . 8 ) 
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E = r E - T El 
z z y J_ 


where E | | is the component parallel to and is the component in the 


y-z plane perpendicular to B. Thus, by recalling that p .>>p , the dif- 

Is & 


ference between the drift velocities of the ions and electrons can be writ- 
ten as 


p i E x *1 

V x. V x 2 P i W x + r ' z W y B p i B 

t e (1+p • ) L o o 

*ls 


(I. 10) 


p .r r e ei~| r E\ | 

y - v = — — - w + p.t w + p JL+ A - -2- (i. li) 

II . ii x i z u ID D p O 

='•» a *» flin 1 V 0 0 3 0 


’e (1+p/-) 


V - V = 

i e (1+P^ ) 


w -p.ru - p . 


E 

X 

5" 1 

r * s ll 

5 

o 

B oJ 

°e B o 


( 1 . 12 ) 


The condition that j =0 requires that v - v =0. Thus, E can 

z z i z e 


be eliminated between equation (I. 11) and (1.12), giving 

p . r . p . e 

— /i ii #*• 


k • i p . E , E\ 

i 1 t x 1 1 

V • v ,.2,. r x % y T B Y B 
s e (1+p . ) z 3 z o z o 


(1.13) 


Equations (I. 10) and (1.13) are now substituted into equations (1.2) 
through (1.4), to give 


E r E I 

f(r W + dz = ~ 

iEe z y B n Es X B n 


(1.14) 


f W E_l 

L - rr'* ■ 


°K * rr )<& 

Es y a o z 


(I. IS) 
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where 


n p. 


F = 


G = 


1 +P . 


n p ; 


i+p. 


It can be assumed that E and E i vary slowly with altitude and can be 

x 1 

treated as constants. Equations (1.14) and (I.1S) ere then solved for 
and , giving 
E 

(1.16) 


where 


rjrf'/V®* 1 * W'-arj 

o j *g L * 


~T - TTT [^4 * - T rS% - ^>] 

o r *g L - 1 


f = 


Fdz 


(1.17) 


9 = 


4 = 


Gdz 


J Er 


Fw dz, etc. 
x 


Es 


The vertical ion velocity can now be found by substituting equations 
(1.16) and (1.17) into equation (1.7). By setting equation (1.12) equal to 
zero, the parallel component of the polarization electric field is found to 
have a negligible effect on the vertical ion velocity. Thus, 

r 


v = 
z 


(l+p/urV 




p.(u - w ) + r (w 


p.(u - w )+r(u - w ) 

i 3 Xj z K y y G J 


2 

9 + 


1 - u )~| / 2 

y y F s 


(£ - w. ) - p.r ( 


x„ 


l 3 ' y p 


ti . i ' 1 

b _ 


J * 


(1.13) 
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where 


w = 7 -, etc. 

x r f 


The field strength can be estimated for the conditions in the iono- 


sphere. Suppose the ion density is constant above 80 km at a value of 
10^ cm \ and zero below 80 km. Also, let w be 100 m sec * between 110 and 

X 

120 km and zero elsewhere, w is set equal to zero everywhere. It can be 
assumed that p. decreases exponentially with altitude, and has a value of 


1 at 120 km. Thus, 


, 2-120 km. 
Pv = exp ( u ) 


(1-19) 


v ‘ ' H 

where H is the scale height of the neutral atmosphere, and is roughly equal 


to 10 km. By using these approximations, equations (1.16) and (1.17) give 


the following: E 


10 -4 V m ' 1 and ~ 10 " 3 V m *. Equation (1.12) then 


gives Em - 10 5 V m *. In midlatitudes, therefore, the polarization 


electric field strength, and its effect on ion velocities, is negligibly 
small . 

At the magnetic equator T =0. The step in the above derivation which 

2 

gives equation (1.13) involves a division by r , end cannot be used. Setting 

z 


T z equal to zero gives 


v - v 
x . x 
t e 


P- £ E \ 

17 + r - it* 


d + pi ) 


( 1 . 20 ) 


v - v 


- v 


-r s || 

P B 

e c 


p .r 
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n 2. x 

( 1+ P i ) 
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E x E L 
n “ zH 


d-21) 

( 1 . 22 ) 


Setting v r 


v = 0 , 
2 

e 
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E 

x 

°iT 
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(1.23) 
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APPENDIX II. 

APPROXIMATE DERIVATIONS OF REFLECTION COEFFICIENTS 
1 1.1 Approxination of a Layer by Multiple Thin Film 

A first approximation of the reflection coefficient is found by de- 
scribing the reflecting medium by a series of vertically-thin layers of 
differing electron density. Within each layer the electron density, and 
therefore the refractive index, is assumed to be everywhere constant. A 
derivation of the reflection coefficient of multiple dielectric films can be 
found in many optics texts [ie. Stone 3 1963]. 

Considering a boundary between two layers, Stokes theorem insures that 
the tangential E and H fields be continuous. Denoting the amplitudes of 
the incident, reflected, and transmitted fields by the subscripts i s r 3 
and t respectively, this cu . be written 


E. ♦ E = 


a. * h = h. 

t V t 


Or, expressing the spatial dependence explicitly, 

n ikn.z .. i kn.z .. -ikn^z 
n .e 1 + h e 1= He 2 

t v t 

Where n ^ and are the refractive indices of the two layers. From 
Maxwell's equations, it is seen that for the transverse components 


(11. 1) 

( 11 . 2 ) 


(II. 3) 


Thus 


H « 


dE 

3a 




(II. 4) 


‘l^t “ r J 2 t 

Equations (II. 1) and (11.5) can be combined to give the Fresnel formulas 
for the reflection (#) and transmission (T) coefficients 


(II. 5) 


I 


E n. - n- 

R = JL- _i 2 

n ~ E. ~ n, + n. 

i 1 2 


(II. 6) 


- - 2 "i 

£. n. ♦ n_ 
r 1 2 


(II. 7) 


If this layer is followed by another boundary and another layer with 
index of refraction n^, there are waves travelling in both directions on 
both sides of the initial boundary. Within the layer, the waves are denoted 


by and E . Thus equations (II. 1) and (II. 5) become 
E. + E = E + E 

i T + 


nAE. - E ) = nAE - E ) 
1 v % 2 V + - 


(11. 8) 

(11. 9) 


At the second boundary, the phase of the wave differs by kn^d from the 
wave at the first boundary, where d is the width of the layer. Thus, the 
above equations are 


E e- ihl 2 d . E e ikn 2 d . £, 
♦ - t 


(II. 10) 


( „ -ikn~d „ ikn^ „ 
n 2 (E + e 2 - E e 2 ) = n 


(II. 11) 


The fields between the boundaries can be eliminated and the result 
written in matrix form 


E. ♦ E 
t r 

n Ai ■ V 


(11.12) 


where 


cos {kn^d) — sin [kn^d) 

2 

in 2 sin(fen d) cos {kn^d) 


This process can be repeated for successive layers. If the index of 


refraction in the j layer is »i., and in the region of the incident and 

J 


i 


I 


transmitted waves, n. and n respectively, equation (11.12) can be general- 
ized to give 


E. + E 
t r 


n .[E . - E )| 
v % r'J 


= A A A A 

A ** ~«7 


where 


cos (fen .d) — sin (fen .d) 

J ^7 «7 

4m 

in^sinffen^.cf) cos (fen .d) 


The product of the /I's can be written in the form of a two-by-two 


matrix 


A. A- ... A = 
-1-2 -n 


c dj 


from which reflection and transmission coefficients can be derived as 
an .+ bn. n - a - dn^ 

r, ^ i u t 


an . * bn .n, + c + an , 
v ^ t t 


an . * bn .n, + c + cin^ 
x v t t 


( 11 . 16 ) 


1 1. 2 W.K.B. Solution 


The wave equation for E can be written 


U . Lj i ~ 

— J4- + fe n £ 
j 2 i 


(11.17) 


In a varying medium, the refractive index is a function of s, making this 
a difficult equation to solve. The derivation of the W.K.B. solut it- 1 is as 


follows [ Budden , 1961]: 
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Let 


so that 


<j>0) = ± iknz 


i${z) 


E = E e 

y yo 


This, substituted into the wave equation, (11.17), gives 


E e t<P = 0 
yo 


or 


. d 2 * 

(dt V ,2 2 

I*;?- 

[lEj * kn 


(d$\ 2 .2 2 . d 2 $ 

U) ~^2 


d 2 <t> . 


It is assumed that — ir is small, and so, to a first approximation, 
dz 


(11.18) 


(11.19) 


( 11 . 20 ) 


( 11 . 21 ) 


3* - ± kn 
dz 


This is differentiated to give 


d 2 <S> _ , dn 

O ” — » > 7 

& 2 02 


which is substituted back into (11.21) to give 


( 11 . 22 ) 


d<f> 

Iz 


.22 dn 

- k n ± vk -T- 
dz 


Then, by expanding the square root, 

\ 2 kn / 


(11.23) 


(11.24) 


The second sign in equation (11.25) is determined by the sign in 
(11.22). Therefore, if the two equations are consistent, both signs 
equation (11.25) must be the same, or, 

<2<j> ^ i dn 

di ±kn * ST 3? 

This may now be integrated to give 


(11.25) 

equation 

in 

(11.26) 
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$ = ± k ndz + ilrifn* 5 ) 

so that 

E ~-E n\ ±ikIndz 

y yo 

The reflection coefficient at a reflection point is given by 
the ratio of the downward travelling to the upward travelling wave. 


(II. 27) 


(11.28) 


|fl| = exp j-2fcj ° 


(11.29) 


where x is the imaginary part of the refractive index and z q is the altitude 
where reflection occurs. 

This approximation does not account for partial reflections from the 
ionized medium. If the frequency of the incident radiation is greater than 
the peak plasma frequency of the medium, the reflection coefficient derived 
by this method is zero. 

An expression for the reflection coefficient where total reflection 
does not occur is suggested by Yeh and Liu [1972]. They solve the wave 
equation (11.17) for E using a Green's function to obtain 

E (z) ,z 

-E = IW ex l >[ - il: ”^*1 

yo v J J o 

i r 2 fr 2 a * r 2 ' £0') 1 


1 f 2 f rZ , rZ' t (3') ■) 

* 2n(z) ex Pt~ ife n(t)dx]^ exp [ik n(T)dt]-| dz\ 

1 J c U-® >o yo J 

tz c r-¥°° >Z' E ( 2 ^) -n 

- 2n(z) ex P^ fe n(T)dT] J exp [-ik n(T)dt]-| dz'\ 

V 1 O } 0 VO J 


(11.30) 


where z ' and t are dummy height variables. 


The first and second terms on the right-hand side describe upward- 
propagating waves, while the third term describes a downward-propagating 


! 


I 


wave. As z -*• ® the second term goes to zero and the upgoing and downgoing 
waves can be compared to give the reflection coefficient. 


r jr rz E ( z ) 

■jg exp[-ik n{x)dx ]~ | dz 

-OO o yo 


(11.31) 


Thus, by substituting the expression for E given by the W.K.B. method 

y 

(11.28), an approximate expression for the reflection coefficient can be 
written in terms of known quantities. 


R = 



rz 


exp [-2 ik 


n(-r) dx]dz 


J o 


(11.32) 


At 2 = +®, the third term on the right-hand side of equation (11.30) 
vanishes, leaving only upward propagating waves. This can be compared with 
the upgoing wave at z = -<= to give the transmission coefficient. 


T = 


?»(-») 
n (+°°) 


exp [-ik 


n(x) dx] 




dn/dz exp[2i/c 




n[x)dx]dz 


(11.33) 


or 


2 ■ " 


(T)dT](l - S) 


(11.34) 


The expression for E , equation (11.28), can be substituted into the 

y 

wave equation (11.17), giving 


3 

( 1 dn'\ ,2 2 i 1 

d 2 n 

4 

^3T) - k ” ’in 

dz 2 . 


E + k 2 n 2 E = 0 

y y 


(11.35) 


This provides a condition for the region in which the W.K.B. approximation 
is valid. That is, 


3 (1 dn 
Tin 3z 


« k 2 n 2 


1_ 3 f J_ 

,7 4 1 2 3a. 


2 3 i 2 

n dz 


This condition is not violated in sporadic-# layers, except where the 
frequency of incident radiation approaches the reflection point in the 
medium. However, this fact makes it impossible to examine the reflection 
coefficient of sporadic-# layers by the W.K.B. method at frequencies near 
and below the peak plasma frequency of the layer. 

1 1.5 Solution Using Hyperyeometric Functions 

This method assumes the ionization profile to be in the shape of a 
so-called Epstein profile [ Budden , 1961]. For this case, exact algebraic 
expressions for the reflection and transmission coefficients can be de- 
rived [5rekkovskikh, I960]. 

Consider the hypergeonetric equation, 

_ (a ♦ B ♦ IK - Y dF aS F = 0 rj1 T ' 


5(1 - 0 


it ~ m - o 


Now iet 


where 


F = r(z)#, 


5 = 5 ( 3 ) 


•(a) = r o 5' v/2 (1 - 


and £(z) is an arbitrary function. Then equation (11.37) becomes 
tf 2 # 

+ = 0 
dz a 




A. 




[ 


where 


2 

... 1 d In dt - 

9W 1 ~71 ln an 

dz '• 


1 d . dl 

4 S ln & 


d i 2 r i 

Tz H [*1 + X 2 R + X 3 


(1-0 


*1 = 4 Y(Y-2) 


K ? = j [1 - Cot-B) 2 «■ YtY-2)] 


K 3 =j= [ (ot+B-y) 2 - 1] 


For equations (11.41) to describe electromagnetic waves in the 
ionosphere, it is necessary that 


. •. ,2 2 

g{z) = k n 


This is accomplished by letting 


COO = -e 


Then , 


2 - r mz mz 

g{z) - - m K. - K~ — K. — — =■ 

a 4 1 2 , , mz n 3 , , m3, 2 
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Now let m, M, and N be defined such that K^, and K 3 can be expressed as 



(11.47) 


With this substitution, (11.44) and (11.46) can be combined to give 


mz 


mz 


n = 1 - N 


(1+e ) 


4 M 


, . mz, 2 
(1+e ) 


(11.48) 


By setting either M or U equal to zero, equation (11.48) is the 
equation for either a symmetrical layer or a transition region between 
different refractive indices. This is illustrated in Figure 3.4. 

The system of equations (11.43) can be solved for a, B, and y giving 

y = 1 ± (1+4 K x ) h 

u-S = ± [-4 K 2 + 1 + y(y-2)]* (11.49) 

a+B = ±(4K +1K + y 

To make the following derivation physically meaningful, the positive 
sign in the above equations will be retained. Thus, 

a = Ml + U+4/^) 15 + (1 + 4/^)^ + [l+4(* 1 -* 2 )] 1 ' 5 } 

8 = Ml + (1+4*^) ^ + (l+4^ 3 ) i5 - [1 +4(K 1 -A: 2 ) ]**} (11.50) 

Y = 1 + ( 1 + 4^) 15 



f 


This can be simplified by defining a relative thickness of the layer, 

S =• 2 - 
m 


and substituting into equation (11.47) 


(1+4-Kj) ^ = iS 


[1 + 4(X 1 -^ 2 )] 15 = iS(l-N) h 


( 1 + 4 ^)^ = (1 -AS 2 M)^ 


Therefore 


a = Js(l + (1 -tS 2 ^ + i | [1 + (l-AO* 2 ]} 


B = %(1 + (i- 45 2 W) 1s + i | [1 - (1-tf)* 2 } 


Y = 1 + iS 


To solve for the reflected and transmitted wave fields, it is now 
necessary to find solutions to equation (11.41). This solution can be 
written in the form of a hypergeometric series, 

F(a B y 0 = 1 - Si £ + ap.+ l)B(B+T) ^ 

na,P»Y»W 1 y * + 2 ! y CY+1) * 

( 

. a(a+l) (a+2) B (B+l) (8*2) 3 

3! Y (y+D(Y+2) * 


It will be recalled that for the case in question, £ = -e . The 
hypergeometric equation (11.37) has poles at the points £ = 0, 1, °°. 
Solutions are sought where z = « (or £ = -°°) and z = (or £ = 0) . Near 

each pole there are two independent convergent hypergeometric series 
which are analytic functions. In the regions near each of these points, 
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Bvckhcvskikh [1960] gives the two convergant series as the following: 


For t. ' 0 


For 


F 1 = r(a,8,Y,0 

F, = S 1_Y F(a-yl, B-y*1. 2-y,5) 


F_ = F(a,B,a+B-Y+l, 1-0 


F 4 = (l-U Y_a ' B F(Y-a,Y-8,Y-a,6+l,l-0 


For 


F. = r B F(a,a-Y+1, a-6+l,C _1 ) 

F 6 = r 6 F(8,8 -y + 1> B-o+1,5" 1 ) 


(11.55) 

(11.56) 

(11.57) 

(11.58) 

(11.59) 

(11.60) 


The solution in the region of £ = -« can be expressed in terms of 
and bv the method of analytic continuation. This is done by solving the 

Barns' contour integral for hypergeometric functions in the two regions 
[Whittaker and Watson, 1952], 


r - H r(a*s) r(fc+s) r(-s) f ,s , 

2 ^ } _ xi r (a*e) L V 


(11.61) 


By Cauchy's Theorem, equation (11.57) can be expressed as a contour 
integral with the solution 

I = 2 (residues within the contour). (11.62) 

The poles in the gamma function occur where the argument of the function is 
equal to zero or a negative integer. The residue at a pole of the gamma 
function is 

residue (T(-nj) - n = 0,1,... (11.63) 

The contour in (11.57) is closed about the right half of the complex plane 
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if | c | <1 and about the left half plane if |t;| > 1. Taking first the case 

of id < i. 


j i ■ 


lim 


N 

l 

N -*■ ® n-0 


r(q+n) T(fc+n) (-1) 
T(c+n) n! 


(-0 


(11.64) 


or 


T _ r r(q+n)r(fc+n) n 
n! r(c+n) ^ 


n=0 


(11.65) 


This is proportional to the hypergeometric series, 


_ r(q)rQb) 
i " r(e) 


F(a,fc,c,0 


( 11 . 66 ) 


Now, the contour is closed about the left half plane (|;j > 1). This 
includes poles in r(q+s) and T(i>+s) at s = -a-n and s = -b-n respectively. 
The solution is 


^2 


lim 


N 


i (-l) 4 r(fc-q-n)r(q+n) .-a-n 

n t 0 «! T(q-q-n) 1 j 


V (-U n r(a-b-n) r(b+n) . .,-b-n 
n t 0 n! T(e-b-n) 


(11.67) 


This can be rewritten in the form of (11.59) using the relation 


Thus , 


F(s) T(l-s) = -r 


■n 


sin tts 


( 11 . 68 ) 


_ _ p (-l) n r(q-*-n) r(l-o+q+n) sin [Tr(q-q-n) ] , . -q-n 

1 2 ~ _A n n! r(l-b+q+n) sin fTT(q-q-n) 1 ^ 


n=0 


[n(t)-q-n) ] 


r (-l) n r (b+n) (1 -c+b+n) sin [n(g-p-n) ] . n -b-n f 
L n! T( 1-q+i+n) sin [ir(q-i-n) ] L- ^ J ul,oyj 


n=0 


or, by again using equation (11.68), 
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I 2 3 (-0"* F[a,\-c+a ,\-b+a 


♦ (-0 ~ b F(b,l-c+b,l-a+b,f l ) I 

The expressions and 1^ can be expressed in the form developed 
earlier by making the following substitutions: 


b = a-y+1 

( 

a = a-B+1 

t - r 1 

so that, recalling equations (11.55), (11.56), and (11.59), setting I 


equal to gives 


Fr = (-1) 


-a r(a-S-l)r(l-Y) 


r(o-y*l) (1-6) 1 


F , ♦ (- 1 ) 


y-l-a r(a-B-*-l)r(Y-l) 


r(a)r(v-B) r 2 

(11.72) 


The function F is related to the electric field through (11.58). 
Making use >f the fact that = 

y r(a) o s v J 


(11.73) 


When £ -*• • (or z -*• «) , the first term of is dominant, so that, by 
recalling equation (11.45), 


E y = r o~ l m ~ ^ exp[- m(a-B) 2 ] 


(11.74) 


When £ -► 0 (or 2 -*■ -«) , the right side of equation (11.72) is used. 
Again, only the first term in each series remains. 

U ” E - r -V s (-ijWTT-n-a 

y o 

9 -CD V 


r(a-B»l)r(l-Y) >(Y- 1)2 a r(q-B*l)r(Y-l) -^(y-I)*! 


x r i?.V v ~ .w a—'" ~ J ~ + ' v ~ r zj • 1 1 e 

Lr(a- Y +i)r(i-B) r (a) r ( y -6) 
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If, as required by (11.53), the expressions (q-B) and (y-1) are 
chosen to be purely imaginary and positive, equation (11.74) becomes an 
expression for a plane wave at z - ® travelling in the +2 direction. 

Equation (11.75) can be interpreted as two plane waves at z = -<*>, one 
travelling in the +2 direction and the other in the -2 direction. Thus, 
(11.75) can be considered to be the incident and reflected waves below the 
ionized region; and (11.74), the transmitted wave above the region. The 
reflection and transmission coefficients are given by ratios of these waves, 


_ r(q-g+i)r(i- Y ) . r(q)r( Y -B) 
r(q-Y+l)r(l-B) r(q-B+l)r(Y-l) 

_ rq-Y) r(qjr(Y-B) 
r(q- Y +i) rq-B)r(Y-i) 

_ r(q) T(y-B) 

" T(q-B+1) r(Y-l) 


(11.76) 

(11.77) 


By substituting the appropriate values of q, B, and y into (11.77), 
the reflection coefficient for a transition or a symmetrical layer are 
found. 

For a transition layer, M = 0, N t 0, and 
q = 1 + If [1 + (1 -N) k ] 

B = 1 + [1 " (11.78) 


Y = 1 + iS 


From which, 

r(-is)r{i + y~[\ * (l-^) 15 ]} r{^ [i + q -N)*]} 
r(i 5 )r{i - y- [1 - U-n) H ]} r{- H- [1 - (l-tf) 4 ]} 


(11.79) 


Thus, from equation (11.68), 
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, sinh { y- U - (l-N)' 2 ]} 

1*1 = {RB*V = ^ r— 

sinh [1 + (1 -NV]) 


(11.80) 


From (11.48) it is seen that the index of refraction a great dis- 
tance above the transition, n , is given by 

00 

n = (1-tf) 1 * (11.81) 

GD 

Whence, 


i*l . sinh ^ 

sinh [£- (1+nJ] 

For a symmetrical layer, M t 0, N = 0, and 

a = k + iS + (k - S 2 M)^ 

8 = k * (k - S 2 M) h 
Y = 1 * iS 


(11.82) 


(11.83) 


and 


r(-jg)r(»s * d * is) r(k - d + iS) 
r (is) r(»s ♦ d) rc*s - *) 

where d = {k - S 2 M )\ 

If 4 S 2 M < 1, d is real, and 


(11.84) 
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(11.85) 


If 4S M > 1, d is imaginary, and 


l*| 2 - 


cosh (tt d) 
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( 11 . 86 ) 
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